GitHub user samlamont created a discussion: Can I register a custom CRS to a user-defined EPSG code to use with RS_SetSRID?
Hello! I am reading netcdf files from s3 and saving locally as geotiff using Sedona and PySpark. If I do not explicitly set the SRID with `RS_SetSRID` I get an error saying it must be defined when writing to disk. However it seems `RS_SetSRID` only accepts an integer EPSG code and the netcdf data are in a non-standard CRS with no EPSG code. I do know the proj4/WKT strings of the CRS. So is there a way I can register the WKT string of the netcdf data to a user-defined EPSG code (ex. 99999) to use with the `RS_SetSRID` in Sedona? The netcdf data is the National Water Model v3.0 retrospective rainfall data, ex: https://noaa-nwm-retrospective-3-0-pds.s3.amazonaws.com/index.html#CONUS/netcdf/FORCING/1979/ I'm using PySpark 4.0, Sedona 1.8.0, and the GeoTools extension (org.datasyslab:geotools-wrapper:1.8.0-33.1) Ex. read: ```python raster_sdf = ( spark .read .format("binaryFile") .load(s3_glob_pattern) .selectExpr("RS_FromNetCDF(content, 'RAINRATE', 'x', 'y') as raster", "path as filepath") .selectExpr("RS_SetSRID(raster, 99999) as raster", "filepath") # <-- can I set 99999 to a custom WKT string? .selectExpr("RS_AsGeoTiff(raster, 'LZW', 1) as raster", "filepath") ) ``` Then writing to geotiff like: ```python raster_sdf.write.format("raster").mode("overwrite").partitionBy("value_time").save(f"{output_dir.as_posix()}") ``` Any tips/guidance would be greatly appreciated, thanks! GitHub link: https://github.com/apache/sedona/discussions/2604 ---- This is an automatically sent email for [email protected]. To unsubscribe, please send an email to: [email protected]
