The PAMAP dataset comes divided between two different projections: Pennsylvania North State Plane and Pennsylvania South State Plane.
I’d like a great big DEM of the whole state of Pennsylvania. gdal_merge.py
can help me make this happen, but I’m going to need to reproject at least one half of the dataset first.
Is there a standard way of choosing how to do such a reprojection? Should I map everything to the North or South Plane or, perhaps, construct an intermediate projection and map everything to that?