Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

create ascii file from tileplot variable from the command line and obtain projection information for GIS from shapefile with same name #332

Open
lizadams opened this issue Jan 9, 2024 · 0 comments
Assignees

Comments

@lizadams
Copy link
Contributor

lizadams commented Jan 9, 2024

Describe the bug
Projection information is not saved when saving variable information from a tile plot using ascii format from the command line.
It is available at the top of the file when saving from the gui.

To Reproduce
Steps to reproduce the behavior:

  1. Create a tile plot ./verdi.command -f $cwd/data/model/CCTM47_aac_spr06.12k.CONC.2006113.O3.ncf -s "O3[1]" -g tile
  2. Click on File > Export As Image/GIS
  3. Select name in Save As text box as ozone
  4. Select ASCII as the File Format
  5. Exit

Expected behavior
cat ozone.asc
When saving the ascii file from the gui, the top of the file contains the following information.

ncols 98
nrows 110
xllcorner -1.41600e+06
yllcorner -696000
cellsize 12000.0
NODATA_value -9999.0
0.0673905 0.0670518 0.0678343 0.0714933 0.0721743 0.0723521 0.0723628 0.0716080 0.0712018 0.0713046 0.0715053 0.0706815 0.0707447 0.0698938 0.0679396 0.0678840 0.0702346 0.0719731 0.0710497 0.0705922 0.0707587 0.0701048 0.0695556 0.0680646 0.0682396 0.0683228 0.0677354 0.0677105 0.0672405 0.0680479 0.0671431 0.0665232 0.0665867 0.0652284 0.0641645 0.0609215 0.0660184 0.0729519 0.0731665 0.0731936 0.0740846 0.0730665 0.0731765 0.0739741

When you save using the following command line option
./verdi.command -f $cwd/data/model/CCTM47_aac_spr06.12k.CONC.2006113.O3.ncf -s "O3[1]" -g tile -save2ascii $cwd/ozone.cli.asc -quit

Then the file output is of a different format:
./verdi.command -f $cwd/data/model/CCTM47_aac_spr06.12k.CONC.2006113.O3.ncf -s "O3[1]" -g tile -save2ascii $cwd/ozone.cli.asc

cat ozone.cli.asc

"Tile Probe: O3[1] (1, 1, 1 - 98, 1 - 110)"
Y,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98
1,0.05153941363096237,0.05198071897029877,0.05232636258006096,0.052306897938251495,0.05240609124302864,0.05262250080704689,0.05268995463848114,0.05316111817955971,0.05362606793642044,0.0538458451628685,0.05428161099553108,0.05476885288953781,0.05482390150427818,0.055143069475889206,0.056226808577775955,0.05736694484949112,0.05755554512143135,0.0583982989192009,0.05914485454559326,0.059180308133363724,0.05912591889500618,0.05896548926830292,0.058850549161434174,0.05749578773975372,0.05633653700351715,0.057069189846515656,0.05715515464544296,0.05722464621067047,0.05804986134171486,0.05760649964213371,0.057398147881031036,0.05892327055335045,0.059149064123630524,0.05908510088920593,0.060123711824417114,0.06012614816427231,0.06034191697835922,0.06086967512965202,0.06102697178721428,0.061074480414390564,0.06082504987716675,0.06093691661953926,0.060828980058431625,0.0628865584731102,0.06395664811134338,0.06394921988248825,0.06481121480464935,0.06536665558815002,0.06351026147603989,0.06292349100112915,0.06313914805650711,0.06322331726551056,0.06291322410106659,0.06284542381763458,0.06290999799966812,0.06393613666296005,0.06387770920991898,0.06408362835645676,0.06613340973854065,0.06606398522853851,0.06596691161394119,0.0679367333650589,0.06741441786289215,0.06764499098062515,0.0693579614162445,0.06972607225179672,0.06976162642240524,0.07185278832912445,0.07270196080207825,0.0723111554980278,0.070078045129776,0.0700555220246315,0.0699421614408493,0.06914334744215012,0.0689915120601654,0.06891454756259918,0.0678599551320076,0.06828218698501587,0.06843023747205734,0.06877056509256363,0.06857455521821976,0.06888900697231293,0.06884121894836426,0.06887107342481613,0.06875975430011749,0.0676107257604599,0.06759228557348251,0.06732814013957977,0.06318158656358719,0.062214527279138565,0.06228265538811684,0.06323616951704025,0.06338950991630554,0.06337607651948929,0.064142607152462,0.06412360817193985,0.06378137320280075,0.06221222132444382

  1. We need to fix the bug so that the ascii format exported by the GUI and the CLI are the same.

If you export the file to shapefile using GUI you get the projection information. IF you export to ascii from the GUI with the same name as shapefile that ascii get the projection information from the shapefile and its opens in right position using any GIS software. When saving to ascii it is a raster format. When saving to a shapefile it is a vector format which shows what the projection is.

  1. We need to update the documentation to give instructions on how to save both the shapefile and the ascii file to allow users to use the ascii in a GIS system.

This makes save to ascii a two step process, to save the projection information to shapefile and then save to ascii using the same file name. The goal is to load ascii and get the correct projection information to use in any GIS system.

  1. Or, when we save to an asc file, we need to also save the projection file.

For this example, the ozone.prj file contains:

PROJCS["unknown", GEOGCS["GCS_unknown", DATUM["D_unknown", SPHEROID["unknown", 6370000.0, 0.0]], PRIMEM["Greenwich", 0.0], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH]], PROJECTION["Lambert_Conformal_Conic_2SP"], PARAMETER["central_meridian", -97.0], PARAMETER["latitude_of_origin", 40.0], PARAMETER["standard_parallel_1", 45.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["standard_parallel_2", 33.0], UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH]]

@lizadams lizadams self-assigned this Jan 25, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant