<div dir="ltr"><div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:#0b5394">Dear Mr. Micha,</div><div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:#0b5394">Yes, it is EPSG: 102022,  I used the proj4 string to create the location. <br></div><div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:#0b5394">Many thanks for the solution. It worked. However, the v.out.ogr module exports the file without the CRS.</div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:#0b5394">Projection information updated<br>(Mon Apr  4 16:44:33 2022) Command finished (0 sec)                             <br>(Mon Apr  4 16:44:45 2022)                                                      <br>g.proj -p                                                                       <br>-PROJ_INFO-------------------------------------------------<br>name       : Africa_Albers_Equal_Area_Conic<br>datum      : wgs84<br>ellps      : wgs84<br>proj       : aea<br>lat_0      : 0<br>lon_0      : 25<br>lat_1      : 20<br>lat_2      : -23<br>x_0        : 0<br>y_0        : 0<br>no_defs    : defined<br>-PROJ_SRID-------------------------------------------------<br>SRID       : EPSG:102022<br>-PROJ_UNITS------------------------------------------------<br>unit       : meter<br>units      : meters<br>meters     : 1</div></blockquote><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div style="font-family:times new roman,serif;font-size:large;color:rgb(11,83,148)" class="gmail_default">v.out.ogr input=MDG_UNSALB_admin3_selectionSN5@PERMANENT output=/home/amit/Documents/mdg_unsalb_east_mdg.gpkg format=GPKG<br><span style="color:rgb(255,0,0)">proj_create: crs not found</span><br>Exporting 7 areas (may take some time)...<br>v.out.ogr complete. 7 features (Polygon type) written to <MDG_UNSALB_admin3_selectionSN5> (GPKG format).<br>(Mon Apr  4 16:45:39 2022) Command finished (0 sec)                             </div><br></div></blockquote><div><span class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:rgb(11,83,148)">Kind regards,</span></div><div><span class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:rgb(11,83,148)">Amit</span> </div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 4 Apr 2022 at 16:33, Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello:<br>
<br>
<br>
On 04/04/2022 13:38, Amit Ghosh wrote:<br>
> Hi,<br>
> I tried to export a vector file to use it in another application. The <br>
> module returns an error message that says "Unable to create OGR <br>
> spatial reference". The detailed output is here:<br>
><br>
>     v.out.ogr input=MDG_UNSALB_admin3_selectionSN5@PERMANENT<br>
>     output=/home/amit/Documents/mdg_aoi2_.gpkg<br>
>     ERROR 10: Pointer 'hSRS' is NULL in 'OSRImportFromWkt'.<br>
>     ERROR: Unable to create OGR spatial reference<br>
><br>
><br>
> Projection info<br>
><br>
>     g.proj -p<br>
>     -PROJ_INFO-------------------------------------------------<br>
><br>
<br>
I would guess that this is the problem:<br>
<br>
>     name       : unknown<br>
>     datum      : wgs84<br>
>     ellps      : wgs84<br>
>     proj       : aea<br>
>     lat_0      : 0<br>
>     lon_0      : 25<br>
>     lat_1      : 20<br>
>     lat_2      : -23<br>
>     x_0        : 0<br>
>     y_0        : 0<br>
>     no_defs    : defined<br>
>     towgs84    : 0.000,0.000,0.000<br>
>     -PROJ_UNITS------------------------------------------------<br>
>     unit       : meter<br>
>     units      : meters<br>
>     meters     : 1<br>
><br>
><br>
> I cannot figure out where I am wrong. Can you suggest an alternative?<br>
<br>
How did you initially create the GRASS Location? Is the albers equal <br>
area projection described by a standard EPSG code? It appears to be this <br>
one:<br>
<br>
<a href="https://epsg.io/102022" rel="noreferrer" target="_blank">https://epsg.io/102022</a><br>
<br>
<br>
**IF** that's the correct projection, then you should be able to do:<br>
<br>
<br>
echo <br>
"PROJCS["Africa_Albers_Equal_Area_Conic",GEOGCS["GCS_WGS_1984",DATUM["WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["False_Easting",0],PARAMETER["False_Northing",0],PARAMETER["longitude_of_center",25],PARAMETER["Standard_Parallel_1",20],PARAMETER["Standard_Parallel_2",-23],PARAMETER["latitude_of_center",0],UNIT["Meter",1],AUTHORITY["EPSG","102022"]]" <br>
| g.proj -c wkt=-<br>
<br>
<br>
<br>
That long WKT string will reset the projection definition of the <br>
Location, then you should be able to export OK.<br>
<br>
<br>
><br>
> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
<br>
-- <br>
Micha Silver<br>
Ben Gurion Univ.<br>
Sde Boker, Remote Sensing Lab<br>
cell: +972-523-665918<br>
<br>
</blockquote></div>