Hi,
I have this calculation here to calculate the distance between 2 points on a globe with their latitudes and longitudes.
[arclen, az] = distance([40.333, -10.036 ], [40.333, -9.46]);km = deg2km(arclen);display(km)
This displays:
>> stationdistkm = 48.8236
However, I have multiple latitudes and longitudes like so:
loc1 = [40.333 -10.036;40.333 -9.46;40.333 -9.767;40.333 -12.219;41.383 -13.888;42.581 -15.461;43.78 -17.032;45.05 -18.505;46.544 -19.672;48.039 -20.848;49.529 -22.017;50.278 -22.603;53.019 -24.752;55.506 -26.71;57.004 -27.879;58.207 -29.725;58.843 -31.267;59.102 -33.828;59.363 -36.397;59.623 -38.954;59.773 -41.297;59.902 -43.015;59.823 -42.399;59.799 -42.004;59.753 -45.112;59.434 -45.666;59.068 -46.083;56.916 -47.422;55.842 -48.093;53.692 -49.433;53 -51.1;]
I need to perform the calculation between each row, then each distance needs to be a sum of itself and previous distances. It should look like this, assuming only 4 columns are used and that stations are calculated to be roughly 50km apart.
>> stationdistkm = 48.8236 96.3574 156.3246 201.4127
How would I create this loop with my given matrix?
Charlie
Best Answer