<div dir="ltr"><div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:#0b5394">It shows:</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-size:large;color:rgb(11,83,148)"><span style="font-family:monospace"><a href="http://v.info">v.info</a> MDG_UNSALB_admin3_selectionSN5@PERMANENT <br> +----------------------------------------------------------------------------+<br> | Name: MDG_UNSALB_admin3_selectionSN5 |<br> | Mapset: PERMANENT |<br> | Location: southern_africa_grass |<br> | Database: /run/media/amit/InternalHD/FAO/country |<br> | Title: |<br> | Map scale: 1:1 |<br> | Name of creator: amit |<br> | Organization: |<br> | Source date: Tue Mar 22 17:47:24 2022 |<br> | Timestamp (first layer): none |<br> |----------------------------------------------------------------------------|<br> | Map format: native |<br> |----------------------------------------------------------------------------|<br> | Type of map: vector (level: 2) |<br> | |<br> | Number of points: 0 Number of centroids: 7 |<br> | Number of lines: 0 Number of boundaries: 27 |<br> | Number of areas: 7 Number of islands: 2 |<br> | |<br> | Map is 3D: No |<br> | Number of dblinks: 1 |<br> | |<br> | Projection: Africa_Albers_Equal_Area_Conic |<br> | |<br> | N: -2374864.5344628 S: -2811333.16032851 |<br> | E: 2422280.61532197 W: 2266547.34604055 |<br> | |<br> | Digitization threshold: 0 |<br> | Comment: |<br> | |<br> +----------------------------------------------------------------------------+</span></div></blockquote></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, 4 Apr 2022 at 17:20, 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">
<div style="direction:ltr">
<p><br>
</p>
<div>On 04/04/2022 14:29, Amit Ghosh wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:rgb(11,83,148)">Dear Mr. Micha,</div>
<div class="gmail_default" style="font-family:times new roman,serif;font-size:large;color:rgb(11,83,148)">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:rgb(11,83,148)">Many thanks for the
solution. It worked. However, the v.out.ogr module exports the
file without the CRS.</div>
</div>
</blockquote>
<p><br>
</p>
<p>What does <font face="monospace"><a href="http://v.info" target="_blank">v.info</a>
MDG_UNSALB_admin3_selectionSN5 </font>show ??</p>
<p><br>
</p>
<blockquote type="cite">
<div dir="ltr">
<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:rgb(11,83,148)">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" target="_blank">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>
</blockquote>
<pre cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</div>
</blockquote></div>