I am using New York's PLUTO lot data in QGIS for a project.
Right now, the shp file automatically imports with a custom CRS, allowing it to visually line up with a separate layer of subway entrances (with that subway layer using WGS84).
That all works fine visually, with both layers overlaying each other on my map.
However, when I attempt to do some analysis in python on my PLUTO features via
print feature.geometry().centroid().asPoint() I find that instead of a lat/lon pair, I get a weird set of x/y coordinates with values in the hundreds of thousands.
I assume that this is related to some sort of "state plane" or "local plane" coordinate reference system.
What I would rather see is a set of coordinates closer to lat 40 and lon -73.
Is there any way to convert my PLUTO layer so it will act more like a WGS84 layer?
Here is the CRS that gets displayed in the "Layer Properties" dialogue box:
Selected CRS (USER:100000, * Generated CRS (+proj=lcc +lat_1=40.66666666666666 +lat_2=41.03333333333333 +lat_0=40.16666666666666 +lon_0=-74 +x_0=300000 +y_0=0 +ellps=GRS80 +towsgs48=... (at which point the screen cuts off the extended text).
Simply using the "save as" dropdown on the layers palette to resave the layer as a WGS84 shapefile does not solve this problem. The new shapefile generated by that process still imports with a custom CRS, and analysis of that generated shapefile still returns coordinates like
x:997248.532881, y:221815.926017. So this does not seem to be a duplicate of the flagged-as-duplicated problem.