add ability to handle extrapolations for both hemispheres

[peterkuriyama]

No description provided.

[James-Thorson]

Thanks Peter!

Why do you have this line:

Data_Extrap_s$N_km <- Data_Extrap_s$N_km - 10000

...?

[James-Thorson]

Also have you tried using grid_in_UTM=FALSE instead?

[peterkuriyama]

I believe it's because converting lat/long in the southern hemisphere results in northings that are in the northern hemisphere? @HaikunXu can you elaborate or modify this explanation?

I have tried using grid_in_utm=FALSE. The problem is that the output list has zone=NA, which then results in an error in the Convert_LL_to_EastNorth_Fn call downstream. I think The CRS.new object in Line 29 of the function has to have a zone (e.g. 12 or 13) to do the calculation.