<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
<title></title>
</head>
<body bgcolor="#ffffff" text="#000000">
Dear Hamish and Yves,<br>
<br>
Results OK on one hand and screwed up on the other:<br>
<br>
<blockquote type="cite">
<pre wrap="">In the corineurope location try
g.region rast=fc
and then go back into your franceLIIe mapset and try again.</pre>
</blockquote>
has not been conclusive, BUT<br>
<blockquote type="cite">
<pre wrap="">another trick is in your source location (corine) zoom around your map
and run v.in.region to make a box, then in Lambert location use v.proj
to import the box, then "g.region vect=box" to get the zoom ready for
the raster map.</pre>
</blockquote>
Worked excellent...<br>
<br>
So definitely one can go a step further with this improved definition
of GRASS regions (thanks for the trick, will be most helpful,
Hamish), <br>
<br>
HOWEVER,<br>
<br>
The raster obtained failed to fall into the right coordinate range...
which may well fit with Yves' remark:<br>
<blockquote type="cite">
<pre wrap="">This not a lambert 2 etendu region. 1 000 000 diffents to the west (I may be wrong)</pre>
</blockquote>
I am a bit loosing my mind on that. The shapefiles given by BD carto
have a coordinate range and projections given to be "Lambert II etendu"
(eg Besancon is 2256000 880000) and corine landcover has Lambert
Azimutal
Equal Area. I have parametrized projections from indications given on
Corine and IGN notices respectively. To sum up, within the Corine Land
Cover I have got this
region for the Franche-Comte window:<br>
<br>
GRASS 6.0.2 (corineurope):~ > g.region -p<br>
projection: 99 (Lambert Azimuthal Equal Area)<br>
zone: 0<br>
datum: ** unknown (default: WGS84) **<br>
ellipsoid: grs80<br>
north: 2781245<br>
south: 2571341<br>
west: 3952172<br>
east: 4110683<br>
nsres: 100.00190567<br>
ewres: 100.00694006<br>
rows: 2099<br>
cols: 1585<br>
<br>
An acceptable window in for France in Lambert II etendu
would be<br>
north: 2677442<br>
south: 1620430<br>
west: 47679<br>
east: 1197823<br>
<br>
When I project from what I believe to be Corine
Land Cover to what I believe to be Lambert II etendu, I get this region
for the Franche-Comte region:<br>
<br>
GRASS 6.0.2 (franceLambertIIe):~ > g.region -p<br>
projection: 99 (Lambert Conformal Conic)<br>
zone: 0<br>
datum: ** unknown (default: WGS84) **<br>
ellipsoid: clark80<br>
north: 5781377.23822138<br>
south: 5575874.90520341<br>
west: 5197476.01104252<br>
east: 5391683.44065733<br>
nsres: 100.00113529<br>
ewres: 100.00382575<br>
rows: 2055<br>
cols: 1942<br>
<br>
Which is really crazy...<br>
<br>
Any suggestion welcome,<br>
<br>
Best regards,<br>
<br>
Patrick<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
Hamish a écrit :
<blockquote cite="mid20070326223050.16dfc17c.hamish_nospam@yahoo.com"
type="cite">
<pre wrap="">Patrick Giraudoux wrote:
</pre>
<blockquote type="cite">
<pre wrap="">GRASS 6.0.2 (franceLIIe):~ > r.proj input=fc location=corineurope
mapset=giraudoux output=fc
</pre>
</blockquote>
<pre wrap=""><!---->
In the corineurope location try
g.region rast=fc
and then go back into your franceLIIe mapset and try again.
another trick is in your source location (corine) zoom around your map
and run v.in.region to make a box, then in Lambert location use v.proj
to import the box, then "g.region vect=box" to get the zoom ready for
the raster map.
good luck,
Hamish
</pre>
</blockquote>
<br>
</body>
</html>