Skip to content

Commit

Permalink
Revert "Added divr and div to destination maps"
Browse files Browse the repository at this point in the history
This reverts commit 15a9ca6.
  • Loading branch information
Crispin committed Feb 2, 2021
1 parent 15a9ca6 commit dbc9e42
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 23 deletions.
20 changes: 6 additions & 14 deletions sDNA/sdna_vs2008/calculation.cpp
Expand Up @@ -604,10 +604,9 @@ void SDNAIntegralCalculation::process_origin(SDNAPolyline *origin,int r,shared_p
//this is very approximate, but simplifies the appearance of sdna output measures:
//if we don't estimate origin diversion ratio then that breaks the concept of normalizing by weight in radius
//as the weight does include the origin itself
double origin_crow_flight = origin->estimate_average_inner_crow_flight_distance_ignoring_oneway(proportion_of_origin_within_radius) ;
total_crow_flies(origin,r) += origin_crow_flight * destweight_of_origin_within_radius;
double origin_diversion_ratio = origin->estimate_average_inner_diversion_ratio_ignoring_oneway() ;
diversion_ratio(origin,r) += origin_diversion_ratio * destweight_of_origin_within_radius;
total_crow_flies(origin,r) += origin->estimate_average_inner_crow_flight_distance_ignoring_oneway(proportion_of_origin_within_radius)
* destweight_of_origin_within_radius;
diversion_ratio(origin,r) += origin->estimate_average_inner_diversion_ratio_ignoring_oneway() * destweight_of_origin_within_radius;

//each point experiences intra-origin traffic of 1/6*(origweight*destweight) in each direction
betweenness.increment_both(origin,r,origin_self_betweenness_weight/6.);
Expand All @@ -632,10 +631,7 @@ void SDNAIntegralCalculation::process_origin(SDNAPolyline *origin,int r,shared_p
(float)origweight_of_origin_within_radius,
(float)two_stage_betweenness_origin_weight,
(float)(cost_of_origin_within_radius / 3.),
(float)(length_of_origin_within_radius / 3.),
(float)origin_crow_flight,
(float)origin_diversion_ratio,
(float)((origin_diversion_ratio-1)*(length_of_origin_within_radius / 3.))),
(float)(length_of_origin_within_radius / 3.)),
1));
shared_ptr<sDNAGeometryPointsByEdgeLength> origasdestination(new sDNAGeometryPointsByEdgeLength());
orig_as_destinationdata->add_part(origin_geom);
Expand Down Expand Up @@ -813,8 +809,7 @@ void SDNAIntegralCalculation::process_geodesic(DestinationEdgeProcessingTask &de
}
total_geodesic_length_weighted(origin_link,r) += euclidean_cost_of_geodesic * dest_weight;

double diversion_ratio_this_destination = euclidean_cost_of_geodesic / crow_flies_distance;
diversion_ratio(origin_link,r) += diversion_ratio_this_destination * dest_weight;
diversion_ratio(origin_link,r) += euclidean_cost_of_geodesic / crow_flies_distance * dest_weight;
#ifdef _SDNADEBUG
cout << "euclidean geodesic from link " << origin_link->arcid << "-" << destination_link->arcid << " is " << euclidean_cost_of_geodesic << endl;
cout << "crow flies distance is " << crow_flies_distance << endl;
Expand Down Expand Up @@ -851,10 +846,7 @@ void SDNAIntegralCalculation::process_geodesic(DestinationEdgeProcessingTask &de
(float)origin_weight,
(float)two_stage_betweenness_weight,
(float)get_geodesic_analytical_cost(dest,anal_best_costs_reaching_edge),
(float)euclidean_cost_of_geodesic,
(float)crow_flies_distance,
(float)diversion_ratio_this_destination,
(float)(euclidean_cost_of_geodesic-crow_flies_distance)),
(float)euclidean_cost_of_geodesic),
1));
shared_ptr<sDNAGeometryPointsByEdgeLength> destination(new sDNAGeometryPointsByEdgeLength());
destination->add_edge_length_from_start_to_end(dest.routing_edge,dest.length_of_edge_inside_radius);
Expand Down
11 changes: 2 additions & 9 deletions sDNA/sdna_vs2008/calculation.h
Expand Up @@ -706,20 +706,16 @@ class SDNAIntegralCalculation : public Calculation, public SDNAPolylineDataSourc
FIELD(fFLOAT,"Weight","Weight");
return result;
}
static vector<SDNAVariant> get_destination_data(long orig,long dest,double radius,float origweight,float tpweight,float anal_cost,float radial_cost,
float crowflight,float divr,float divdiff)
static vector<SDNAVariant> get_destination_data(long orig,long dest,double radius,float origweight,float tpweight,float anal_cost,float radial_cost)
{
vector<SDNAVariant> result(10);
vector<SDNAVariant> result(7);
result[0]=SDNAVariant(orig);
result[1]=SDNAVariant(dest);
result[2]=SDNAVariant((float)radius);
result[3]=SDNAVariant(origweight);
result[4]=SDNAVariant(tpweight);
result[5]=SDNAVariant(anal_cost);
result[6]=SDNAVariant(radial_cost);
result[7]=SDNAVariant(crowflight);
result[8]=SDNAVariant(divr);
result[9]=SDNAVariant(divdiff);
return result;
}
static vector<FieldMetaData> get_destination_field_metadata()
Expand All @@ -732,9 +728,6 @@ class SDNAIntegralCalculation : public Calculation, public SDNAPolylineDataSourc
FIELD(fFLOAT,"TPWeight from origin","TPWeightO");
FIELD(fFLOAT,"Analytic Metric","AnalMet");
FIELD(fFLOAT,"Length","Length");
FIELD(fFLOAT,"Crow flight","CrowFlt");
FIELD(fFLOAT,"Diversion Ratio","DivR");
FIELD(fFLOAT,"Diversion","Div");
return result;
}

Expand Down

0 comments on commit dbc9e42

Please sign in to comment.