[Qgis-user] map projection in Myanmar

[Leo Kris Palao]

Good Day QGIS Users,

Would like to ask what map projection is suitable for Myanmar? Is it okay
to use UTM coordinate system or can you recommend other coordinate system?
Thank you.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] map projection in Myanmar

[Leo Kris Palao]

Dear Wirote,

Thank you very much. It will surely help.

-Leo


On Mon, Jan 28, 2013 at 6:38 PM, Wirote Laongmanee wrote:

> Dear Leo,
>
> In Thailand, we are using UTM Zone 47N for majority part and 48N at east
> part of Thailand. And also we are using WGS84-Latlon together, we also
> remain Indian Datum 1975 data but need to consider when work together.
>
> In Myanmar I am not sure there are specific projection or not but when we
> handle their data we selected WGS84/UTM Zone 46 N, it work well.
> May help you.
>
> Wirote
>
>
> On Mon, Jan 28, 2013 at 5:27 PM, Leo Kris Palao wrote:
>
>> Good Day QGIS Users,
>>
>> Would like to ask what map projection is suitable for Myanmar? Is it okay
>> to use UTM coordinate system or can you recommend other coordinate system?
>> Thank you.
>>
>> -Leo
>>
>> ___
>> Qgis-user mailing list
>> Qgis-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>>
>>
>
>
> --
> Wirote Laongmanee
> http://wmap.wordpress.com
> http://parn-site.blogspot.com
>
>


-- 
***

*LEO KRIS M. PALAO*
Researcher
World Agroforestry Centre (International Centre for Research in
Agroforestry)
Mobile #: +63(916)7698-367
E-mail: l.pa...@irri.org or lk.pa...@gmail.com

*Transforming Lives and Landscapes
*
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Projection

[Leo Kris Palao]

Good Day QGIS Users!

I am working on a GIS project where my area of interest (drainage basin and
river network) cover two UTM Zones. In this case, I cannot use UTM to
compute for area and distance, and to perform other analysis in which
spatial data needs to be projected. What projection can I use? Is Equal
Area Cylindrical NSIDC EASE - Grid Global (EPSG: 3410) okay to use if
working on a very large spatial dataset? Is this widely used by GIS users?
Thank you.

-Leo Kris
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] TauDEM

[Leo Kris Palao]

Dear QGIS Users,

Can I ask somebody to share me notes on how to install TauDEM in QGIS. My
TauDEM seems to be not working correctly. This is how I install TauDEM.

1. I go to fetch phyton plugin in and install Sextante.
2. I install TauDEM 5.0.6 Command Line Executables for 64 bit Windows
PC for
ArcGIs 10.
3. I install MPICH2
4. I configure TauDEM under sextante and specify MPICH2/OpenMPI bin
directory (path is C:\Program Files\MPICH2\bin) and TauDEM command line
tools folder (path is C:\Program Files(x86)\TauDEM\TauDEM5Exe)

I tried running taudem but it is not running correctly, it says *Could not
load model watersheds. model ALGORITHM:saga:catchmentarea(parallel)*
*
*
Does anyone has successfully run TauDEM in QGIS? can you share your
installation notes.

I am using Windows 7 prof - 64 bit.

Thanks.
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] TauDEM

[Leo Kris Palao]

Thanks Alex, will try your suggestion! thanks for taking the time to reply.
TauDEM is a good tool for terrain analysis and hope I can make it work.

-Leo


On Fri, Mar 15, 2013 at 3:00 PM, Alexander Bruy wrote:

> Hi Leo,
>
> On Fri, 15 Mar 2013 11:59:04 +0800
> Leo Kris Palao  wrote:
>
> > 1. I go to fetch phyton plugin in and install Sextante.
> > 2. I install TauDEM 5.0.6 Command Line Executables for 64 bit Windows
> > PC<http://hydrology.usu.edu/taudem/taudem5.0/TauDEM5Exex64.zip> for
> > ArcGIs 10.
> > 3. I install MPICH2
> > 4. I configure TauDEM under sextante and specify MPICH2/OpenMPI bin
> > directory (path is C:\Program Files\MPICH2\bin) and TauDEM command line
> > tools folder (path is C:\Program Files(x86)\TauDEM\TauDEM5Exe)
>
> As I can see, you setup all correctly and TauDEM should work with
> this settings.
>
> > I tried running taudem but it is not running correctly, it says *Could
> not
> > load model watersheds. model ALGORITHM:saga:catchmentarea(parallel)*
>
> This message not related to TauDEM, you see this message because
> SAGA provider is not enabled. Please look for another messages
> in SEXTANTE History and Log window.
>
> Also can you try to run any TauDEM tool from pure commandline and
> see if it works? For exmaple simplest one — PitRemove.
>
> --
> Alexander Bruy
>



-- 
***

*LEO KRIS M. PALAO*
Researcher
World Agroforestry Centre (International Centre for Research in
Agroforestry)
Mobile #: +63(916)7698-367
E-mail: l.pa...@irri.org or lk.pa...@gmail.com

*Transforming Lives and Landscapes
*
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] plugin for processing multiple raster files

[Leo Kris Palao]

Hi *QGIS users*,

I am working with daily rainfall raster dataset of one decade acquired from
*TRMM*. We will use this as criteria for our suitability analysis of a
farming system. I want to derive the following:

1. Mean daily rainfall *(*output is 365 raster files (representing *day
1-day 365*) with each raster file contain the *mean daily rainfall* derived
from one decade of daily rainfall data for a specific day*)*; and
 2. After computing the mean daily rainfall for one decade, I want to
compute the rainfall sum with the interval of 10 days (per dekad).

is there a plugin in QGIS that I can use that will automate this process. I
know this is possible using a script, but I do not have a background yet in
scripting. Right now I am still reading Phyton programming. Also I think it
is possible to do this in R, but right now I am still learning R.

Any help will be much appreciated. Thank you very much.

-*Leo Kris Palao*
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] plugin for processing multiple raster files

[Leo Kris Palao]

Thank you Etienne for your help, I'll try your suggestion. BTW, my data is
in netcdf format. I am having trouble reclassifying my data using raster
calculator when using hdf. I'll email you if I get stuck with the
processing =)

-Leo


On Thu, May 2, 2013 at 7:44 PM, Etienne Tourigny wrote:

> Hi,
>
> I am not sure you need GIS tools (such as QGIS) to do this kind of thing,
> which is fairly routin in the meteorological/climate research field.
>
> My advice would be to convert the TRMM data to a single netcdf file (which
> is much more suited for temporal + spatial datasets), and then process the
> data using command-line tools like CDO or NCO. In my experience this could
> be done with just a few commands.
>
> cdo: https://code.zmaw.de/projects/cdo
> nco: nco.sourceforge.net/
>
> BTW in which format did you get the TRMM data? most likely netcdf or hdf,,,
>
> let me know privately if you need a bit more help
>
> cheers
> Etienne
>
> On Thu, May 2, 2013 at 12:59 AM, Leo Kris Palao wrote:
>
>> Hi *QGIS users*,
>>
>> I am working with daily rainfall raster dataset of one decade acquired
>> from *TRMM*. We will use this as criteria for our suitability analysis
>> of a farming system. I want to derive the following:
>>
>> 1. Mean daily rainfall *(*output is 365 raster files (representing *day
>> 1-day 365*) with each raster file contain the *mean daily rainfall*derived 
>> from one decade of daily rainfall data for a specific day
>> *)*; and
>>  2. After computing the mean daily rainfall for one decade, I want to
>> compute the rainfall sum with the interval of 10 days (per dekad).
>>
>> is there a plugin in QGIS that I can use that will automate this process.
>> I know this is possible using a script, but I do not have a background yet
>> in scripting. Right now I am still reading Phyton programming. Also I think
>> it is possible to do this in R, but right now I am still learning R.
>>
>> Any help will be much appreciated. Thank you very much.
>>
>> -*Leo Kris Palao*
>>
>> ___
>> Qgis-user mailing list
>> Qgis-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>>
>>
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Customized projection

[Leo Kris Palao]

Dear QGIS Community,

Good Day!

I am currently working with a project that covers the whole country. My
problem is if I want to compute for perimeter, area, and when using terrain
analysis modules, the spatial datasets should be projected. But the country
covers two UTM zones and I want to use a much broader scale projection. I
want to use the Equal Area Cylindrical > *NSIDC EASE-Grid Global
Projection*but I want to set the central meridian to be in the middle
of my country of
interest. In this case, I can minimize the distortions in my country of
interest, right? I do not know what parameter should I change to do it.
Below is the code for the abovementioned projection.

+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +a=6371228 +b=6371228 +units=m
+no_defs

Thanks in advance.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] QGIS and QT Designer

[Leo Kris Palao]

Dear QGIS Users,

I am exploring the synergy of QGIS and QT on a tablet for field work and
field validation activities. I learn this functionality from Nathan
Woodrow. As I understand, I need QGIS-dev and QT-dev packages to facilitate
the integration of QT forms in QGIS. This require the use of OSGeo4W
installer to install the packages, hence this makes Windows-based tablets
appropriate to use.

Right now, the only way to use the QT designer in QGIS is to have a vector
point layer in which the fields are linked to the widgets in the QT form.
In this case, for QT form to work (*pop-up in the map*) a user should
digitize a point in the map canvass. Is there a way to use the GPS
location/coordinates (*like mark waypoint in GPS*) instead of digitizing my
points? For instance, I can use the current GPS location instead of
manually digitizing my points. Or it can be the other way around, where I
have an existing point layer and when you select on a specific point
feature in the map canvass the QT forms will pop-up and populate or update
the attribute information of that point. Does this makes sense?

By the way (if it is okay), can anybody recommend a good windows-based
tablet (that can support QGIS and QT, longer battery life, and with minimal
glare when used outdoors)?

Thank you.
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Qgis-user Digest, Vol 89, Issue 84

[Leo Kris Palao]

Hi Andreas,

Thanks for your detailed response and suggestions. I will certainly
consider those issues that you raised and suggestions, as well as try
NextGIS. We are really moving towards paper less fieldwork, that is why we
are exploring how tablets and QGIS can be optimized for field use. I have
limited experience in using Ubuntu/Linux OS, and for now I am more
comfortable in using QGIS on a Windows platform. Hence I preferred QGIS on
a Windows tablet. I would love to try QGIS on Ubuntu so everything will be
Opensourced.

I have tried running QGIS on our android device. I installed the QGIS
nightly in Samsung Galaxy Tab. I loaded a road map and administrative
boundary map. Everything is okay when I navigate, but when I click the
*Automatically
add point *in *Tracks* and *Automatically save added feature*, then after a
short navigation then QGIS will stop and exit.

I agree to you that QGIS on a tablet should be further improved. Anyone can
just carry a tablet device instead of the usual paper, ballpen, and GPS.

My preference on using Windows tablet is because I want to use QT forms in
QGIS as a platform for conducting field surveys. Although I have not yet
actually tried it before. But according to the blog and post of Nathan,
QGIS is working well in Windows tablet.

Hi Nathan, I really want to try QMap. Any documents or related web links
that you can share to install QMap (to have the interface like the one you
posted in your email). Right now, what I have tried is to use QT forms in
QGIS on my laptop (just to explore the functionality). I assume that
installing and using QGIS and QT is the same as installing and using it in
PCs and laptops.

Thank you,

-Leo


On Fri, Jul 19, 2013 at 4:19 PM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Re: Geostats tool (bathymetry) (pcr...@pcreso.com)
>2. QGIS and QT Designer (Leo Kris Palao)
>3. Re: QGIS and QT Designer (Andreas Neumann)
>4. Re: QGIS and QT Designer (Nathan Woodrow)
>5. Setting up an FTP server to read a json file (Christophe B)
>
>
> --
>
> Message: 1
> Date: Thu, 18 Jul 2013 12:36:16 -0700 (PDT)
> From: pcr...@pcreso.com
> To: image 
> Cc: "qgis-user@lists.osgeo.org" 
> Subject: Re: [Qgis-user] Geostats tool (bathymetry)
> Message-ID:
> <1374176176.62364.yahoomailba...@web122002.mail.ne1.yahoo.com>
> Content-Type: text/plain; charset=iso-8859-1
>
> I suggest yo look at using GMT to work with these data.
>
> http://gmt.soest.hawaii.edu
>
> turn both datasets into DEM's & compare them (subtract one from the other
> to return an anomaly grid). You can do some of this sort of work with
> GDAL/QGIS, but GMT is optimised for bathymetric data - one of the few GIS
> tools that is used more in aquatic (& geological) domains than typical
> terrestrial GIS.
>
>
>
>  2013/7/17 image 
>
>  Dear  all,
>
>
>  I have 2 bathymetric data : one raster (generated by a black
>  box from an
>
>  image satellite Kompsat2 multispectral) AND a vector data
>  (our reference)
>
>  with a field informing bathymeric values.
>
>  So we would like comparate these 2 data with a geostats
>  tool. Is there a
>
>  tool in qgis to do that (grass? sextante? taudem? saga?
>  Other tool?).
>
>  Moreover, we have a excel file with several rows with
>  several columns
>
>  informing ?: bathy values (reference+ generated data)/
>  Coordinates (image
>
>  pixel row column).
>
>
>
>  If you can tell me what availabale tools could be
>  interesting for me?
>
>
>
>  Thank you very much.
>
>
>
>  Regards.
>
>
>
>
>
>
>
>
>
>
>
>  --
>
>  View this message in context:
> http://osgeo-org.1560.x6.nabble.com/Geostats-tool-bathymetry-tp5067105.html
>
>  Sent from the Quantum GIS - User mailing list archive at
>  Nabble.com.
>
>  ___
>
>  Qgis-user mailing list
>
>  Qgis-user@lists.osgeo.org
>
>  http://lists.osgeo.org/mailman/listinfo/qgis-user
>
>
>
>
>  -Inline Attachment Follows-
>
>  ___
>

[Qgis-user] sextante GRASS r.statistics

[Leo Kris Palao]

Dear QGIS users,

Good Day!

I want to use the *GRASS r.statistics* module in *Sextante Toolbox* to know
the *majority* of *cell values* per *administrative boundary* (Provincial
level). To determine the cell majority per boundary I want to use the
*Mode*method (
The *mode* is the value that appears most often in a set of data
(wikipedia)). Below are my parameters:

*Base raster layer: *PHL_adm_1 (this is a raster administrative boundary
with x = 0.037 and y = 0.041 resolution; Data type: GDT_Float64 - Sixty
four bit floating point)
*
*
*Cover raster layer:* PET_1 (raster Potential Evaportranspiration layer
with x = 0.25 and y = 0.25 resolution, float type; Data type: GDT_Float32 -
Thirty two bit floating point )
*
*
*Method: *mode
*
*
*Cover values extracted from the category labels of the cover map: *Yes
*
*
*GRASS region: *default [Leave blank to use min covering extent]
*
*
*GRASS region cellsize:* default (value is 0.0)
*
*
*Output raster layer:* [Save to temporary file]

*Error: *Could not load layer:
C:\Users\\AppData\Local\Temp\sextante\grassstatistics.tif
Check the SEXTANTE log to look for errors

I do not know what is causing the error. Basically, I need to know the
majority of PET values per Province (PHL_adm_1).

I have configured the GRASS successfully.

*GRASS Folder:* C:\OSGeo4W\apps\grass\grass-6.4.3RC4
*GRASS help folder: *C:\OSGeo4W\apps\grass\grass-6.4.3RC4\docs\html
*Msys folder:* C:\OSGeo4W\apps\msys

Thanks for any help.

Truly yours,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Assign value for polygon layers based on nearby value of point features

[Leo Kris Palao]

Dear QGIS Users,

Good Day!

I have a regularly spaced polygon layer for a country, and a point layer
with precipitation values. I want to assign the value of the precipitation
point layer to the nearest polygon grids.

Any advice will be very helpful.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Vector to Rast

[Leo Kris Palao]

Dear QGIS Community,

Good Day!

I tried using r.mode module of GRASS in Sextante Toolbox. It worked well
but I have to convert my vector layer to raster layer in ArcGIS to generate
an integer raster layer. Whenever I convert my vector layer to raster in
QGIS 1.8 (Win7 Prof) using the Raster > Conversion > Rasterize module it
always gives me a GDT_Float64 - Sixty four bit floating point raster (which
is not a valid input as a Base Layer in the r.mode module), even though the
vector field that I am using is an integer type. How can I do this properly?

Thanks in advance for the help.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Vector to Rast

[Leo Kris Palao]

Hi Yasser,

Thanks for your suggestion. I tried it in Raster Calculator my input is my *
.tif* raster (*GDT_Float64 - Sixty four bit floating point*). My syntax
is *round(.tif)* but it converted my 64-bit raster to GDT_Float32 -
Thirty two bit floating point raster.

Thanks,
-Leo


On Thu, Aug 29, 2013 at 2:07 AM, Yasser Said Lopez de Olmos Reyes <
biolyas...@gmail.com> wrote:

> I've never been in that situation, but you can use r.mapcalculator, after
> rasterization, with round(16bit_raster_coberture) to get it done.
>
>
> 2013/8/28 Leo Kris Palao 
>
>>  Dear QGIS Community,
>>
>> Good Day!
>>
>> I tried using r.mode module of GRASS in Sextante Toolbox. It worked well
>> but I have to convert my vector layer to raster layer in ArcGIS to generate
>> an integer raster layer. Whenever I convert my vector layer to raster in
>> QGIS 1.8 (Win7 Prof) using the Raster > Conversion > Rasterize module it
>> always gives me a GDT_Float64 - Sixty four bit floating point raster (which
>> is not a valid input as a Base Layer in the r.mode module), even though the
>> vector field that I am using is an integer type. How can I do this properly?
>>
>> Thanks in advance for the help.
>>
>> -Leo
>>
>> ___
>> Qgis-user mailing list
>> Qgis-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>>
>
>
>
> --
> Saludos,
>
> Yasser
>



-- 
*LEO KRIS M. PALAO*
GIS - Social Sciences Division
International Rice Research Institute
Website: www.irri.org
Mobile #: +63(916)7698-367
E-mail: l.pa...@irri.org or lk.pa...@gmail.com
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Converting null(no data) to 0

[Leo Kris Palao]

Dear QGIS Community,

Good Day!

We are currently processing srtm_5 minutes resolution data for asia. We
have polygons (125 features) that covers 9 raster cells each. Our task is
to compute the average cell values using the zonal statistics plugin. In
our data we have negative values and NULL values; we managed to mask those
negative values using raster calculator -- our syntax is *srtm@1 > 0) *
srtm@1*. But we are having a problem dealing with the NULL values. We tried
the raster > conversion > translate module and set the No Data value to 0.
It did convert the No Data to 0 values as indicated in the metadata. But
when we clicked on the pixels the value is -3.40282e+38. My observation
was, in the polygons where the pixel values are mix of actual elevation
values and -3.40282e+38, it creates inconsistency in the results. For
instance, we have a polygon that covers 9 cells -- the 7 cells contains
elevation values ranging from 4 to 4.59, and the other two cells have
values -3.40282e+38, the result was: *MEAN = -7.56182992530064e+22* and *SUM
= -6.80564693277058e+22*. If the NULL values have been successfully
converted to 0 the result should be MEAN = 2.90103 and SUM = 26.10925.

How can I do this properly.

Thank you in advance for any suggestion and help.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Qgis-user Digest, Vol 90, Issue 71

[Leo Kris Palao]

Hi Stefan,

Thanks for taking the time reply. It did work the Sextante > Saga > Grid
Tools > Reclassify Grid Values tool.

Thank you.

-Leo


On Fri, Aug 30, 2013 at 4:21 PM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Converting null(no data) to 0 (Leo Kris Palao)
>2. OSX - QGIS not loading Grass plugin (Graeme Merrall)
>3. How to show images of arster and vector files? (Yoshiaki Tsuzuki)
>4. Re: Converting null(no data) to 0 (skampus)
>5. Re: How to show images of arster and vector files?
>   (Alexandre Neto)
>
>
> ----------
>
> Message: 1
> Date: Fri, 30 Aug 2013 13:44:15 +0800
> From: Leo Kris Palao 
> To: QGIS Mailing List , a.mauna...@irri.org
> Subject: [Qgis-user] Converting null(no data) to 0
> Message-ID:
> <
> camktv+2rizzbed107bkcry9kq_h+dox+tzkus1b+5ob5m5k...@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear QGIS Community,
>
> Good Day!
>
> We are currently processing srtm_5 minutes resolution data for asia. We
> have polygons (125 features) that covers 9 raster cells each. Our task is
> to compute the average cell values using the zonal statistics plugin. In
> our data we have negative values and NULL values; we managed to mask those
> negative values using raster calculator -- our syntax is *srtm@1 > 0) *
> srtm@1*. But we are having a problem dealing with the NULL values. We
> tried
> the raster > conversion > translate module and set the No Data value to 0.
> It did convert the No Data to 0 values as indicated in the metadata. But
> when we clicked on the pixels the value is -3.40282e+38. My observation
> was, in the polygons where the pixel values are mix of actual elevation
> values and -3.40282e+38, it creates inconsistency in the results. For
> instance, we have a polygon that covers 9 cells -- the 7 cells contains
> elevation values ranging from 4 to 4.59, and the other two cells have
> values -3.40282e+38, the result was: *MEAN = -7.56182992530064e+22* and
> *SUM
> = -6.80564693277058e+22*. If the NULL values have been successfully
> converted to 0 the result should be MEAN = 2.90103 and SUM = 26.10925.
>
> How can I do this properly.
>
> Thank you in advance for any suggestion and help.
>
> Thanks,
> -Leo
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20130830/6fff0796/attachment-0001.html
> >
>
> --
>
> Message: 2
> Date: Fri, 30 Aug 2013 16:04:14 +1000
> From: Graeme Merrall 
> To: qgis-user@lists.osgeo.org
> Subject: [Qgis-user] OSX - QGIS not loading Grass plugin
> Message-ID:
>  gdcclriyk...@mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> I'm trying to get the Grass plugin to load in QGIS in OSX Mountain Lion.
>
> QGIS ver: 1.8.0 (kyngchaos)
> Grass ver: 6.4 (kyngchaos)
> GDAL ver: 1.10 (kyngchaos with Grass plugin)
> Python: 2.7.2
>
> The error I'm getting is:
> Failed to load
> /Applications/QGIS.app/Contents/MacOS/../PlugIns/qgis/libgrassplugin.so
> (Reason: Cannot load library
> /Applications/QGIS.app/Contents/MacOS/../PlugIns/qgis/libgrassplugin.so:
>
> (dlopen(/Applications/QGIS.app/Contents/MacOS/../PlugIns/qgis/libgrassplugin.so,
> 5): Library not loaded:
> /Applications/GRASS-6.4.app/Contents/MacOS/lib/libgrass_gis.6.4.2.dylib
>   Referenced from:
> /Applications/QGIS.app/Contents/PlugIns/qgis/libgrassplugin.so
>   Reason: image not found))
>
> Following errors include libgrassprovider and libgrassrasterprovider.
>
> >From from searching about I saw very old references to fftw install
> (tried,
> no joy) and tweaking your LD_LIBRARY_PATH.  For LD_LIBRARY_PATH I exported
> both the QGIS plugins dir and Grass lib dir and started QGIS from the
> command line - still didn't work.
>
> Can anyone point out what's going wrong?
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attac

[Qgis-user] Python Script for Multi-band zonal statistics

[Leo Kris Palao]

Dear QGIS Users,

Often we want to compute for band statistics for our stacked layers (i.e.,
bands 1 to 9). In QGIS, the available plugin for this is the Zonal
Statistics plugin but it only computes for one band at a time, and as we
suspect if our input is a stacked layer the plugin seems to get the average
values for all the stacked layers in a polygon.

So, I searched for some answers in the web, and I found a python script. I
used it and it seems to be working well, but this only works for one band,
So, what I did was to copy the script and change the band number for each
copy. However, this resulted in a very long script, and I think there is a
much better script than this. Can someone help me improve my script. I am
just new to Python. I just completed a basic Python training course.

import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(1),1)
zonalstats.calculateStatistics(None)
import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(2),2)
zonalstats.calculateStatistics(None)
import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(3),3)
zonalstats.calculateStatistics(None)
import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(4),4)
zonalstats.calculateStatistics(None)
import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(5),5)
zonalstats.calculateStatistics(None)
import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(6),6)
zonalstats.calculateStatistics(None)
import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(7),7)
zonalstats.calculateStatistics(None)

We used this tool most often to determine the average reflectance or
backscatter from our images per zone (polygon).

Thanks in advance,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] QGIS 1.8.0 Crash open Project Settings

[Leo Kris Palao]

Dear QGIS Users,

Good Day!

I have a problem with my QGIS 1.8.0. I have repeatedly install and
uninstall the program but the behavior is still the same. The QGIS crashes
when I *1.*) open Settings > Project Settings, *2.)* when I save as the
layer and select Layer CRS button or *3.)* Click the CRS Status Icon on the
lower right portion.

When I uninstall QGIS, I erase also *.qgis* folder under Users and Quantum
GIS Lisboa folder under Program Files (x86). I don't know why this crash
happens.

When I install, I also do rebooting to ensure proper installation. I have
not experienced this before. I am using Windows 7 Professional 64-bit OS.

Thanks in advance for any help!

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Sextante/Processing Toolbox

[Leo Kris Palao]

Dear *QGIS *Users,

Good Day!

Just would like to ask how to properly configure sextante/processing
toolbox GRASS and SAGA modules using the QGIS 2.0 standalone installer.

I installed QGIS 2.0 Dufour, GRASS 6.4.3 and SAGA using standalone
installers. Then in the Processing Menu > Options and Configurations I set
the following:

*GRASS Folder:* C:\Program Files (x86)\GRASS 6.4.3
*Msys Folder:* C:\Program Files (x86)\GRASS 6.4.3\Msys

*SAGA Folder:* C:\Program Files (x86)\SAGA-GIS

*I am not quite sure about the actual path name because I already
uninstalled QGIS 2.0 standalone application as well as GRASS and SAGA.*

When I ran each module of GRASS it says that my GRASS folder is not
properly configured.

I installed QGIS, GRASS, SAGA, and Orfeo Toolbox using the OSGeo4W 32-bit
installer and it is working okay.

I just want to ask how can I properly configure the processing toolbox
using QGIS, GRASS, and SAGA standalone installers. Because we are planning
to give training to our partners and I feel that the OSGeo4W will not be
practical to use as installer for the training because we might have
slow/interrupted internet connection.

Thanks in advance for any help.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Qgis-user Digest, Vol 95, Issue 37

[Leo Kris Palao]

Hi Andrew,

Good day!

Thanks for your advice, I'll copy the downloaded files in my USB. BTW, is
there any other way to configure the modules of GRASS and SAGA using
standalone installers? or the best way is to use OSGeo4W installer?

Thanks
-Leo


On Thu, Jan 23, 2014 at 4:00 AM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Sextante/Processing Toolbox (Leo Kris Palao)
>2. Re: Sextante/Processing Toolbox (Andre Joost)
>3. Problem with numerical precision of a raster layer (steffy)
>4. Re: configuraci?n otb y grass (Rafael Torres)
>5. Re: problem reprojecting some shapefiles (Giacomo Fontanelli)
>6. Re: Current version of QGIS (Jonathan Moules)
>7. Re: Problem with numerical precision of a raster layer
>   (Etienne Tourigny)
>8. Re: Current version of QGIS (J?rgen E. Fischer)
>
>
> ----------
>
> Message: 1
> Date: Wed, 22 Jan 2014 10:18:32 +0800
> From: Leo Kris Palao 
> To: QGIS Mailing List 
> Subject: [Qgis-user] Sextante/Processing Toolbox
> Message-ID:
>  oc4b0vndppnmyxkutivjnp89edzy...@mail.gmail.com>
> Content-Type: text/plain; charset="iso-8859-1"
>
> Dear *QGIS *Users,
>
> Good Day!
>
> Just would like to ask how to properly configure sextante/processing
> toolbox GRASS and SAGA modules using the QGIS 2.0 standalone installer.
>
> I installed QGIS 2.0 Dufour, GRASS 6.4.3 and SAGA using standalone
> installers. Then in the Processing Menu > Options and Configurations I set
> the following:
>
> *GRASS Folder:* C:\Program Files (x86)\GRASS 6.4.3
> *Msys Folder:* C:\Program Files (x86)\GRASS 6.4.3\Msys
>
> *SAGA Folder:* C:\Program Files (x86)\SAGA-GIS
>
> *I am not quite sure about the actual path name because I already
> uninstalled QGIS 2.0 standalone application as well as GRASS and SAGA.*
>
> When I ran each module of GRASS it says that my GRASS folder is not
> properly configured.
>
> I installed QGIS, GRASS, SAGA, and Orfeo Toolbox using the OSGeo4W 32-bit
> installer and it is working okay.
>
> I just want to ask how can I properly configure the processing toolbox
> using QGIS, GRASS, and SAGA standalone installers. Because we are planning
> to give training to our partners and I feel that the OSGeo4W will not be
> practical to use as installer for the training because we might have
> slow/interrupted internet connection.
>
> Thanks in advance for any help.
>
> -Leo
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20140122/d6028e4f/attachment-0001.html
> >
>
> ------
>
> Message: 2
> Date: Wed, 22 Jan 2014 06:11:19 +0100
> From: Andre Joost 
> To: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] Sextante/Processing Toolbox
> Message-ID: 
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Am 22.01.2014 03:18, schrieb Leo Kris Palao:
>
> > I installed QGIS, GRASS, SAGA, and Orfeo Toolbox using the OSGeo4W 32-bit
> > installer and it is working okay.
> >
> > I just want to ask how can I properly configure the processing toolbox
> > using QGIS, GRASS, and SAGA standalone installers. Because we are
> planning
> > to give training to our partners and I feel that the OSGeo4W will not be
> > practical to use as installer for the training because we might have
> > slow/interrupted internet connection.
> >
>
> Using OSGeo4W, you can save the downloaded files to disk, save that and
> "Install from disk" from the Advanced Installation menu for the other
> computers.
>
> HTH,
> Andr? Joost
>
>
>
>
> --
>
> Message: 3
> Date: Wed, 22 Jan 2014 01:07:51 -0800 (PST)
> From: steffy 
> To: qgis-user@lists.osgeo.org
> Subject: [Qgis-user] Problem with numerical precision of a raster
> layer
> Message-ID: <1390381671716-5099392.p...@n6.nabble.com>
> Content-Type: text/plain; charset=us-ascii
>
> Hi,
> I've problems to visualize a .tif raster in a singleband pseudocolor map

[Qgis-user] Field calculator

[Leo Kris Palao]

Deal ALL QGIS Users,

I have a strange results when using the Field Calculator. I have an integer
value in one of my column ranging from 1 to 13. I would like to get the
percentage of it when divided by 18. So in field calculator my expression
is "(DN" / 18) * 100. When I execute this expression I got all 0 results in
all cells. My column where the results should be stored is in Real with
Width = 3 and Precision = 3. I did the same in ArcGIS but the results were
correct.

Also, I tried another method, this time I created a column first, then I
just divide the DN values to 18 but still the results are 0 in all cells in
my column.

I am using QGIS 2.2 installed using OSGeo4W installer, with Windows 7 64
bit OS.

Could anyone help me troubelshoot?

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Field calculator

[Leo Kris Palao]

It worked well. Thanks Nathan.

-Leo


On Thu, Mar 27, 2014 at 1:23 PM, Nathan Woodrow  wrote:

> Hey Leo,
>
> Try this:
>
> ("DN" / 18.00) * 100
>
> The field calculator will integer calculations unless told one of the
> values is a real.  Adding .00 at the end is enough to for that.
>
> - Nathan
>
>
> On Thu, Mar 27, 2014 at 3:17 PM, Leo Kris Palao wrote:
>
>> Deal ALL QGIS Users,
>>
>> I have a strange results when using the Field Calculator. I have an
>> integer value in one of my column ranging from 1 to 13. I would like to get
>> the percentage of it when divided by 18. So in field calculator my
>> expression is "(DN" / 18) * 100. When I execute this expression I got all 0
>> results in all cells. My column where the results should be stored is in
>> Real with Width = 3 and Precision = 3. I did the same in ArcGIS but the
>> results were correct.
>>
>> Also, I tried another method, this time I created a column first, then I
>> just divide the DN values to 18 but still the results are 0 in all cells in
>> my column.
>>
>> I am using QGIS 2.2 installed using OSGeo4W installer, with Windows 7 64
>> bit OS.
>>
>> Could anyone help me troubelshoot?
>>
>> Thanks,
>> -Leo
>>
>> ___
>> Qgis-user mailing list
>> Qgis-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>>
>
>


-- 
*LEO KRIS M. PALAO*
GIS and Remote Sensing Specialist
Social Sciences Division - GIS Unit
International Rice Research Institute
Website: www.irri.org
Mobile #: +63(999)-99-75327
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Value Tool Python Error Prompt

[Leo Kris Palao]

Hi QGIS users,

Good day!

I installed/added the Value tool plugin in QGIS 2.2. The installation was
successful. But whenever I enable the Value tool and turn on a layer in the
map legend a Python Error occurs. It doesn't crash (after the error prompt
the value tool still works), but I don't know why there is a Python error
prompt. I did not encounter this one when I am using QGIS 2.0.

Below is the Python Error Prompt:

I am using QGIS 2.2 Valmiera installed via OSGeo4W installer 32 bit. I am
using Windows 7 64 bit OS.

Thanks,
-Leo

An error has occured while executing Python code:

Traceback (most recent call last):
  File "C:/Users/lpalao/.qgis2/python/plugins\valuetool\valuewidget.py",
line 500, in invalidatePlot
if self.mplLine is not None:
AttributeError: 'ValueWidget' object has no attribute 'mplLine'

Python version:
2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]


QGIS version:
2.2.0-Valmiera Valmiera, c3a2817

Python path: ['C:/Users/lpalao/.qgis2/python/plugins\\processing',
'C:/Users/lpalao/.qgis2/python/plugins\\LecoS',
'C:/OSGeo4W/apps/qgis/./python', u'C:/Users/lpalao/.qgis2/python',
u'C:/Users/lpalao/.qgis2/python/plugins',
'C:/OSGeo4W/apps/qgis/./python/plugins',
'C:\\OSGeo4W\\apps\\orfeotoolbox\\wrapping\\lib',
'C:\\OSGeo4W\\apps\\orfeotoolbox\\wrapping\\Python',
'C:\\OSGeo4W\\apps\\orfeotoolbox\\python', 'C:\\OSGeo4W\\bin',
'C:\\OSGeo4W\\bin\\python27.zip', 'C:\\OSGeo4W\\apps\\Python27\\DLLs',
'C:\\OSGeo4W\\apps\\Python27\\lib',
'C:\\OSGeo4W\\apps\\Python27\\lib\\plat-win',
'C:\\OSGeo4W\\apps\\Python27\\lib\\lib-tk', 'C:\\OSGeo4W\\apps\\Python27',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\PIL',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32\\lib',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\Pythonwin',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win32.egg',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode',
'C:\\Users\\lpalao\\.qgis2\\python\\plugins\\mmqgis/forms',
'C:\\OSGeo4W\\apps\\qgis\\python\\plugins\\fTools\\tools']

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Tablet vs GPS accuracy

[Leo Kris Palao]

Hi QGIS Users,

I would like to inquire if there are some user that already assessed the
positioning (location) accuracy of tablet vs commercial gps. We would like
to use qgis android for fieldwork. We would also like to try qmap.

Thanks,
Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Tablet vs GPS accuracy

[Leo Kris Palao]

Hi Guys,

Thanks all for your valuable insights about my inquiry. I really appreciate
your comments. I can use your comments as reference on how we will also
evaluate GIS applications in Tablets.

Our application mainly falls on agriculture. We want to use the tablets for
surveying in Rural areas. And we want to use QGIS android so we can load
our shapefile and raster maps for validation. On the other hand, we will
also want to use tablets to take location of field corners of farmers field
to get the plot size (this is where accuracy is important).

I observe that when there is a wifi the location significantly increases in
tablets. But of course, wifi is not available in rural areas. So, if we
will use tablet we will just rely on available GPS signal, weather, and the
device.

Thank so much, and again I really appreciate the comments,
-Leo


On Fri, Apr 4, 2014 at 2:20 AM, Michael.Dodd wrote:

> I agree with the points Richard has made here.  I mentioned in an earlier
> message on this topic that I'd done a lot of testing with smartphones vs
> consumer gps  see
> http://www.open.ac.uk/blogs/otih/?p=466 for details on some of the tests.
>
> Also tested consumer grade gps and a range of other surveying techniques
> with results and discussion shown here:
> http://oro.open.ac.uk/30066/
> 
> From: Richard McDonnell [richard.mcdonn...@opw.ie]
> Sent: Thursday, April 03, 2014 4:09 PM
> To: Leo Kris Palao
> Cc: QGIS Mailing List
> Subject: Re: [Qgis-user] Tablet vs GPS accuracy
>
> Hi,
> It all boils down to the level of accuracy your job requires, Commercial
> or Survey Grade GNSS (Global Navigation Satellite System) now utilises NRTK
> (Network Real Time Kinematic) which can produce sub centimetre accuracies.
> We have found that the accuracies is often too good, with our surveys
> conflicting with Maps. We have seen +-40mm accuracies in the field, where a
> stake in the ground is surveyed and the kit can take us right back to the
> point days later, to within 40mm. it is worth noting that Survey grade GPS
> has a stated accuracy of sub meter without NRTK.
> The next step down in accuracies would be GPS enabled Field Computers,
> again a dedicated kit but hand held form factor, the Trimble Geo-Explorer 7
> is an example of one of these. They now also have the ability to utilise
> NRTK, bringing there accuracies to within 0.5m (or better).
> The issue with phones and tablets is that, unlike dedicated systems, you
> have less control over how the coordinates are recorded. They don't take
> into account things like multipath (GPS signal bouncing off buildings,
> Ground etc.)  The number of satellites (the more the better) Most survey
> kits wont record a point with less than 5 satellites.
> There is a way to get an idea of the accuracy for your proposed device.
> Take several readings over a day (7 or more) of 2 fixed points approx 50m
> apart,  every day, for a week. Then compare the results. you will have a
> spread of points, none of which are likely to be perfect, but will give you
> an idea of the accuracies of the device, the spread of the points being
> your expected error. A shorter test over a day would yield similar results
> but over a week, you would have a more comprehensive result. I have seen
> some tests, where accuracies with an i-Phone were stated as +-3m.
> One other thing you should all note when using a phone or similar device
> to record your location, Use the native Latitude and Longitude to record
> the position, don't use the device to transform the coordinates to some
> other CRS. Use a dedicated CRS Transformation tool on a desktop/Laptop
> computer to convert to your desired CRS. Handhelds are not suitable for
> that purpose.
> I hope this helps!!
> Regards,
>
> Richard.
>
> On 03/04/2014 12:30, Leo Kris Palao wrote:
>
> Hi QGIS Users,
>
> I would like to inquire if there are some user that already assessed the
> positioning (location) accuracy of tablet vs commercial gps. We would like
> to use qgis android for fieldwork. We would also like to try qmap.
>
> Thanks,
> Leo
>
>
>
> ___
> Qgis-user mailing list
> Qgis-user@lists.osgeo.org<mailto:Qgis-user@lists.osgeo.org>
> http://lists.osgeo.org/mailman/listinfo/qgis-user
>
>
>
> **
>
> Email Disclaimer: http://www.opw.ie/en/disclaimer/
>
> **
>
>
>
> -- The Open University is incorporated by Royal Charter (RC 000391), an
> exempt charity in England & Wales and a charity registered in Scotland (SC
> 038302).
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Cell Statistics in QGIS

[Leo Kris Palao]

Hi ALL QGIS users,

Good day!

I want to calculate the sum of the raster layers, but I want to ignore the
no data values in the calculation. This is similar to cell statistics
(spatial analyst) enable the *Ignore NoData in calculations (optional)* in
ArcGIS.

So far I am doing all my process in QGIS (batch processing is great in the
processing toolbox). This is the only step that I am using ArcGIS because I
cannot find the right tool in QGIS.

Thanks in advance for any help.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Cell Statistics in QGIS

[Leo Kris Palao]

Hi Micha and Victor,

Thanks for your reply. I will try R.Series module of GRASS. I think this is
what I am looking for...

-Leo
On Apr 14, 2014 5:57 PM, "Victor Olaya"  wrote:

> Using the Processing raster calculator to sum up layers should
> correctly handle no-data values, and leave a no-data value in those
> cells for which one of such values appears in the sum.
>
> 2014-04-14 11:29 GMT+02:00 Micha Silver :
> >
> > On 14/04/2014 06:19, Leo Kris Palao wrote:
> >
> > Hi ALL QGIS users,
> >
> > Good day!
> >
> > I want to calculate the sum of the raster layers, but I want to ignore
> the
> > no data values in the calculation. This is similar to cell statistics
> > (spatial analyst) enable the Ignore NoData in calculations (optional) in
> > ArcGIS.
> >
> >
> > You can choose the GRASS module "r.series". Regarding no data cells, they
> > are handled as follows (from the manual page):
> >
> > Without -n flag, the complete list of inputs for each cell (including
> NULLs)
> > is passed to the aggregate function. Individual aggregates can handle
> data
> > as they choose. Mostly, they just compute the aggregate over the non-NULL
> > values, producing a NULL result only if all inputs are NULL.
> >
> > So if you use "method=sum" you'll get an output raster which sums all non
> > null input cells.
> >
> >
> >
> > So far I am doing all my process in QGIS (batch processing is great in
> the
> > processing toolbox). This is the only step that I am using ArcGIS
> because I
> > cannot find the right tool in QGIS.
> >
> > Thanks in advance for any help.
> >
> > -Leo
> >
> > This mail was received via Mail-SeCure System.
> >
> >
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
> > This mail was received via Mail-SeCure System.
> >
> >
> >
> >
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Raster Uniques Values

[Leo Kris Palao]

Hi QGIS users,

Good Day!

I would just like to know the unique values of my raster layer, can you
help me find the right tool?

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Field Calculator Get lon values of points

[Leo Kris Palao]

Dear QGIS users,

Good day!

I want to get the longitude and latitude coordinate values of my points. I
used Geometry > $x and Geometry > $y to derive longitude and latitude of my
points. It give me correct values but when I exit editing mode the values
were rounded to whole number. For instance:

my Lon values is 95.94592 but when I exit editing mode this value is
converted to 95.0

Thanks in advance for any help.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Field Calculator Get lon values of points

[Leo Kris Palao]

Hi Alex,

Thanks for you comment. I have already solved the issue. When I first
created the field I set the width to 3 and precision to 5 - this is where
my values were rounded to whole number. But when I tried it the second
time, I leave the width to its default value of 12 and adjusted the
precision to 5 my value were not rounded (i got the correct values). The
width supposed to be the number of characters for the integer and precision
is the decimal places.

Thanks,
-Leo


On Tue, Apr 29, 2014 at 11:29 AM, Alex Mandel wrote:

> On 04/28/2014 06:54 PM, Leo Kris Palao wrote:
> > Dear QGIS users,
> >
> > Good day!
> >
> > I want to get the longitude and latitude coordinate values of my points.
> I
> > used Geometry > $x and Geometry > $y to derive longitude and latitude of
> my
> > points. It give me correct values but when I exit editing mode the values
> > were rounded to whole number. For instance:
> >
> > my Lon values is 95.94592 but when I exit editing mode this value is
> > converted to 95.0
> >
> > Thanks in advance for any help.
> >
> > -Leo
>
> New field or existing field?
> Hopefully with a new field you picked Text, and changed the maximum
> width. That is assuming you wanted to get 38.345 95.945 etc
> If you wanted to fields, then use Decimal Number(real), and do each one
> at a time. Now when you pick real, you have set the Precision. Precision
> is the number of values after the decimal (to the right). Typical GIS is 8.
>
> So to store 1 coordinate:
> Decimal Field (real)
> Width 12
> Precision 8
>
> Now if it's an existing field, just make a new field because it sounds
> like the old field isn't set right.
>
> If that doesn't help, can you tell use what version of QGIS and what
> operating system(including version).
>
> Thanks,
> Alex
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Field Calculator Get lon values of points

[Leo Kris Palao]

Thanks Alex, it cleared some of my questions. Before I thought that the
width is only for the integer and do not consider the number of decimal
places.

-Leo


On Wed, Apr 30, 2014 at 12:22 AM, Alex Mandel wrote:

> Close, I believe width is total number of characters. And yes precision
> is just the number of characters to the right of the decimal.
> so 12 w 5 p means 7char.5char
> Which is why I suggested 12 w 8 p
> -120.12345678
> 38.12345678
>
> I don't recall if a negative counts as a char, I don't think it does on
> a real. 8 after the decimal is a common amount of precision, the last
> couple of digits take you down to meter/submeter. I'd have to look up
> exactly where it gets you I don't recall right now.
>
> Thanks,
> Alex
>
> On 04/28/2014 10:23 PM, Leo Kris Palao wrote:
> > Hi Alex,
> >
> > Thanks for you comment. I have already solved the issue. When I first
> > created the field I set the width to 3 and precision to 5 - this is where
> > my values were rounded to whole number. But when I tried it the second
> > time, I leave the width to its default value of 12 and adjusted the
> > precision to 5 my value were not rounded (i got the correct values). The
> > width supposed to be the number of characters for the integer and
> precision
> > is the decimal places.
> >
> > Thanks,
> > -Leo
> >
> >
> > On Tue, Apr 29, 2014 at 11:29 AM, Alex Mandel <
> tech_...@wildintellect.com>wrote:
> >
> >> On 04/28/2014 06:54 PM, Leo Kris Palao wrote:
> >>> Dear QGIS users,
> >>>
> >>> Good day!
> >>>
> >>> I want to get the longitude and latitude coordinate values of my
> points.
> >> I
> >>> used Geometry > $x and Geometry > $y to derive longitude and latitude
> of
> >> my
> >>> points. It give me correct values but when I exit editing mode the
> values
> >>> were rounded to whole number. For instance:
> >>>
> >>> my Lon values is 95.94592 but when I exit editing mode this value is
> >>> converted to 95.0
> >>>
> >>> Thanks in advance for any help.
> >>>
> >>> -Leo
> >>
> >> New field or existing field?
> >> Hopefully with a new field you picked Text, and changed the maximum
> >> width. That is assuming you wanted to get 38.345 95.945 etc
> >> If you wanted to fields, then use Decimal Number(real), and do each one
> >> at a time. Now when you pick real, you have set the Precision. Precision
> >> is the number of values after the decimal (to the right). Typical GIS
> is 8.
> >>
> >> So to store 1 coordinate:
> >> Decimal Field (real)
> >> Width 12
> >> Precision 8
> >>
> >> Now if it's an existing field, just make a new field because it sounds
> >> like the old field isn't set right.
> >>
> >> If that doesn't help, can you tell use what version of QGIS and what
> >> operating system(including version).
> >>
> >> Thanks,
> >> Alex
> >>
> >
> >
> >
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
> >
>
>


-- 
*LEO KRIS MARIANO PALAO*
Social Sciences Division
International Rice Research Institute
Website: www.irri.org
Mobile #: 0932-847-6644
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Projection String

[Leo Kris Palao]

Hi ALL,

Good Day!

I have two files: 1) a raster data; 2) a vector fileI would like to
overlay my datasets in QGIS but I cannot seem to get the correct Projection
String.

I extracted the Projection Information of the Raster file in *Raster
*> *Miscellaneous
*> *Information*. And the details are as follows:

Driver: ENVI/ENVI .hdr Labelled
Files: Drive:/Folder/Folder/Folder/Folder/data_h30v07
   Drive:/Folder/Folder/Folder/Folder\filename.hdr
Size is 4800, 4800
Coordinate System is:
PROJCS["Sample Sinusoidal",
GEOGCS["GCS_WGS_1984",
DATUM["WGS_1984",
SPHEROID["WGS_84",6378137.0,298.257223563]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Sinusoidal"],
PARAMETER["False_Easting",0.0],
PARAMETER["False_Northing",0.0],
PARAMETER["longitude_of_center",30.0],
UNIT["Meter",1]]
Origin = (13343406.2363100,2223901.0392980)
Pixel Size = (231.65635825990,-231.65635825990)

The projection string that I created for the above projection is:

*+proj=sinu +lon_0=30 +x_0=0 +y_0=0 +a=6378137.0 +b=298.257223563 +units=m
+no_defs*

I created a custom projection in QGIS, but still it won't overlay. Is there
something wrong with my string?

The *+b *is the semi-minor axis right? If so, the values 298 is I think
weird.

Thanks in advance for any replies,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection String

[Leo Kris Palao]

Hi Even and QGIS users,

Thanks for the reply.

I open the metadata of my raster file and copy the projection string:

+proj=sinu +lon_0=30 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0
+units=m +no_defs

I used the projection string in the metadata to modify my custom
projection. But still it won't overlay...I open the .hdr file of my raster
file, and here is the metadata:

ENVI
description = {
  Applied Mask Result [Mon Oct 21 11:27:39 2013]}
samples = 4800
lines   = 4800
bands   = 368
header offset = 0
file type = ENVI Standard
data type = 1
interleave = bip
sensor type = Unknown
byte order = 0
map info = {Sample Sinusoidal, 1., 1., 13343406.2363, 2223901.0393,
2.3165635826e+002, 2.3165635826e+002, WGS-84, units=Meters}
coordinate system string = {PROJCS["Sample
Sinusoidal",GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Sinusoidal"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",30.0],UNIT["Meter",1.0]]}
wavelength units = Unknown
band names = {

What could be the possible parameter(s) that I did not include in my
projection string?

Thanks again,
-Leo

On Thu, May 8, 2014 at 6:03 PM, Even Rouault
wrote:

>
> > The *+b *is the semi-minor axis right? If so, the values 298 is I think
> > weird.
>
> 298.xxx is a typical value for the reverse flattening, not the semi-minor
> axis.
> So you should use +rf instead of +b
> But as those values are in fact the definition of the WGS84 ellipsoid, you
> could more simply replace the "+a= +rf=" by +ellps=WGS84
>
> --
> Geospatial professional services
> http://even.rouault.free.fr/services.html
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] r.recode in processing toolbox

[Leo Kris Palao]

Dear QGIS users,

Good Day!

I am recoding the values in my raster data. So I use the r.recode module in
the processing toolbox. I am using a recode rule from a text file that says:

I want to replace the value 3 with the new value 1448.8

3:3:1448.8:1448.8

It run completely Okay. No problem with the processing. I just find it
weird that when I select a cell using the Identify Feature it gave me a
value of 1448.800048828125. How can I retain the number of decimal places?

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Link CSV file to Raster (TIFF)

[Leo Kris Palao]

Dear QGIS users,

I would like to ask if it is possible to link the values of a column in a
csv file to a raster TIFF format, then use it to create symbology?

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Share Python Script for Point Analyses Tool

[Leo Kris Palao]

Hi QGIS Users,

I am processing multi-temporal EVI data derived from MODIS, and the number
of datasets could be around 70 images. I have points and want to extract
pixel values from those points. I know that the Point Analyses Tool can do
it. However, it is cumbersome to load all the MODIS data because it takes
too long. Is there a way to extract the pixel values from points without
actually loading the datasets? Maybe a script (GDAL or Python) can do this.
Can somebody help me or point me to the right blog or website where I can
use as reference.

Thanks in advance for the help.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Share Python Script for Point Analyses Tool

[Leo Kris Palao]

Thanks Saber,

I think that is the tool that I need. I will explore how it works.

Thanks very much.
-Leo


On Wed, Jun 4, 2014 at 3:32 PM, Saber Razmjooei <
saber.razmjo...@lutraconsulting.co.uk> wrote:

> Have you tried this:
>
> http://gis-lab.info/qa/extract-values-rasters-eng.html
>
>
>
> Cheers,
>
> Saber
>
>
>
>
>
> *From:* qgis-user-boun...@lists.osgeo.org [mailto:
> qgis-user-boun...@lists.osgeo.org] *On Behalf Of *Leo Kris Palao
> *Sent:* 04 June 2014 06:50
> *To:* QGIS Mailing List
> *Subject:* [Qgis-user] Share Python Script for Point Analyses Tool
>
>
>
> Hi QGIS Users,
>
> I am processing multi-temporal EVI data derived from MODIS, and the number
> of datasets could be around 70 images. I have points and want to extract
> pixel values from those points. I know that the Point Analyses Tool can do
> it. However, it is cumbersome to load all the MODIS data because it takes
> too long. Is there a way to extract the pixel values from points without
> actually loading the datasets? Maybe a script (GDAL or Python) can do this.
> Can somebody help me or point me to the right blog or website where I can
> use as reference.
>
> Thanks in advance for the help.
>
> -Leo
>
>
> --
>
> This email and any files transmitted with it are confidential and intended
> solely for the use of the individual or entity to whom they are addressed.
> If you have received this email in error please notify the system manager.
> This message contains confidential information and is intended only for the
> individual named. If you are not the named addressee you should not
> disseminate, distribute or copy this e-mail. Please notify the sender
> immediately by e-mail if you have received this e-mail by mistake and
> delete this e-mail from your system. If you are not the intended recipient
> you are notified that disclosing, copying, distributing or taking any
> action in reliance on the contents of this information is strictly
> prohibited.
>
> Whilst reasonable care has been taken to avoid virus transmission, no
> responsibility for viruses is taken and it is your responsibility to carry
> out such checks as you feel appropriate.
>
> If this email contains a quote or offer to sell products, carry out work
> or perform services then our standard terms and conditions
> <http://www.lutraconsulting.co.uk/downloads/Lutra%20Consulting%20Standard%20Terms%20and%20Conditions.pdf>
> shall apply unless explicitly stated otherwise.
>
> Saber Razmjooei and Peter Wells trading as Lutra Consulting.
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Share Python Script for Point Analyses Tool

[Leo Kris Palao]

Hi Saber,

Have you tried to use the tool before that you suggested me
http://gis-lab.info/qa/extract-values-rasters-eng.html . I tried to run it
but I get an error message - I think I got the input-output argument right.
Can you help me debug it?

I am running it in five MODIS-EVI dataset first for testing. I am using
QGIS 2.2 with GDAL and GDAL-Python installed using OSGeo4W installer. and I
am using Windows 7 64 bit professional OS. I have the extract_values.py
script, rasters, and shapefile in one folder.

Thanks,
-Leo

GDAL 1.10.1, released 2013/08/26

C:\>cd F:\MODIS_Historical\First_season\yr2012\h27v07\Test

C:\>F:

F:\MODIS_Historical\First_season\yr2012\h27v07\Test>extract_values.py -c
GT_and_
farmsurvey_modis.shp -d F:\MODIS_Historical\First_season\yr2012\h27v07\Test
Traceback (most recent call last):
  File "F:\MODIS_Historical\First_season\yr2012\h27v07\Test\
extract_values.py",
line 34, in 
import gdal, ogr, osr
  File "C:\OSGeo4W\apps\Python27\lib\site-packages\gdal.py", line 2, in


from osgeo.gdal import deprecation_warn
  File "C:\OSGeo4W\apps\Python27\lib\site-packages\osgeo\__init__.py", line
21,
in 
_gdal = swig_import_helper()
  File "C:\OSGeo4W\apps\Python27\lib\site-packages\osgeo\__init__.py", line
17,
in swig_import_helper
_mod = imp.load_module('_gdal', fp, pathname, description)
ImportError: DLL load failed: %1 is not a valid Win32 application.

F:\MODIS_Historical\First_season\yr2012\h27v07\Test>


On Wed, Jun 4, 2014 at 3:32 PM, Saber Razmjooei <
saber.razmjo...@lutraconsulting.co.uk> wrote:

> Have you tried this:
>
> http://gis-lab.info/qa/extract-values-rasters-eng.html
>
>
>
> Cheers,
>
> Saber
>
>
>
>
>
> *From:* qgis-user-boun...@lists.osgeo.org [mailto:
> qgis-user-boun...@lists.osgeo.org] *On Behalf Of *Leo Kris Palao
> *Sent:* 04 June 2014 06:50
> *To:* QGIS Mailing List
> *Subject:* [Qgis-user] Share Python Script for Point Analyses Tool
>
>
>
> Hi QGIS Users,
>
> I am processing multi-temporal EVI data derived from MODIS, and the number
> of datasets could be around 70 images. I have points and want to extract
> pixel values from those points. I know that the Point Analyses Tool can do
> it. However, it is cumbersome to load all the MODIS data because it takes
> too long. Is there a way to extract the pixel values from points without
> actually loading the datasets? Maybe a script (GDAL or Python) can do this.
> Can somebody help me or point me to the right blog or website where I can
> use as reference.
>
> Thanks in advance for the help.
>
> -Leo
>
>
> --
>
> This email and any files transmitted with it are confidential and intended
> solely for the use of the individual or entity to whom they are addressed.
> If you have received this email in error please notify the system manager.
> This message contains confidential information and is intended only for the
> individual named. If you are not the named addressee you should not
> disseminate, distribute or copy this e-mail. Please notify the sender
> immediately by e-mail if you have received this e-mail by mistake and
> delete this e-mail from your system. If you are not the intended recipient
> you are notified that disclosing, copying, distributing or taking any
> action in reliance on the contents of this information is strictly
> prohibited.
>
> Whilst reasonable care has been taken to avoid virus transmission, no
> responsibility for viruses is taken and it is your responsibility to carry
> out such checks as you feel appropriate.
>
> If this email contains a quote or offer to sell products, carry out work
> or perform services then our standard terms and conditions
> <http://www.lutraconsulting.co.uk/downloads/Lutra%20Consulting%20Standard%20Terms%20and%20Conditions.pdf>
> shall apply unless explicitly stated otherwise.
>
> Saber Razmjooei and Peter Wells trading as Lutra Consulting.
>
>
>


-- 
*LEO KRIS MARIANO PALAO*
Social Sciences Division
International Rice Research Institute
Website: www.irri.org
Mobile #: 0932-847-6644
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Share Python Script for Point Analyses Tool

[Leo Kris Palao]

Hi Saber,
The code does not want long shapefile filename - so I rename my shapefile.
I also modifed the argument not to take the whole directory but specific
files. Since my raster data are in ENVI format, it throws an error when the
code reads the .hdr files.

This is my new argument that I executed in OSGeo4W shell:

E:\Test_modis>python extract_values.py -c pts.shp 2013041 2013049 2013057
2013065 2013073

But it still throw an error message. When I look at the attribute table of
the pts.shp it created new fieldnames based on raster filename but all the
values are NULL.

Thanks,
-Leo

this is the error message that I got:

E:\Test_modis>python extract_values.py -c pts.shp 2013041 2013049 2013057
201306
5 2013073
Found GDAL version: 1.10.1

ERROR 4: `pts.shp' not recognised as a supported file format.

Unable to open shapefile -c

E:\Test_modis>python extract_values.py pts.shp 2013041 2013049 2013057
2013065 2
013073
Found GDAL version: 1.10.1

[   ] \
0.1%Tr
aceback (most recent call last):
  File "extract_values.py", line 411, in 
inFeat.SetField( shortName, value )
  File "C:\OSGeo4W\apps\Python27\lib\site-packages\osgeo\ogr.py", line
2572, in
SetField
return _ogr.Feature_SetField(self, *args)
NotImplementedError: Wrong number of arguments for overloaded function
'Feature_
SetField'.
  Possible C/C++ prototypes are:
SetField(OGRFeatureShadow *,int,char const *)
SetField(OGRFeatureShadow *,char const *,char const *)
SetField(OGRFeatureShadow *,int,int)
SetField(OGRFeatureShadow *,char const *,int)
SetField(OGRFeatureShadow *,int,double)
SetField(OGRFeatureShadow *,char const *,double)
SetField(OGRFeatureShadow *,int,int,int,int,int,int,int,int)
SetField(OGRFeatureShadow *,char const *,int,int,int,int,int,int,int)


On Thu, Jun 5, 2014 at 9:19 AM, Leo Kris Palao  wrote:

> Hi Saber,
>
> Have you tried to use the tool before that you suggested me
> http://gis-lab.info/qa/extract-values-rasters-eng.html . I tried to run
> it but I get an error message - I think I got the input-output argument
> right. Can you help me debug it?
>
> I am running it in five MODIS-EVI dataset first for testing. I am using
> QGIS 2.2 with GDAL and GDAL-Python installed using OSGeo4W installer. and I
> am using Windows 7 64 bit professional OS. I have the extract_values.py
> script, rasters, and shapefile in one folder.
>
> Thanks,
> -Leo
>
> GDAL 1.10.1, released 2013/08/26
>
> C:\>cd F:\MODIS_Historical\First_season\yr2012\h27v07\Test
>
> C:\>F:
>
> F:\MODIS_Historical\First_season\yr2012\h27v07\Test>extract_values.py -c
> GT_and_
> farmsurvey_modis.shp -d F:\MODIS_Historical\First_
> season\yr2012\h27v07\Test
> Traceback (most recent call last):
>   File "F:\MODIS_Historical\First_season\yr2012\h27v07\Test\
> extract_values.py",
> line 34, in 
> import gdal, ogr, osr
>   File "C:\OSGeo4W\apps\Python27\lib\site-packages\gdal.py", line 2, in
> 
>
> from osgeo.gdal import deprecation_warn
>   File "C:\OSGeo4W\apps\Python27\lib\site-packages\osgeo\__init__.py",
> line 21,
> in 
> _gdal = swig_import_helper()
>   File "C:\OSGeo4W\apps\Python27\lib\site-packages\osgeo\__init__.py",
> line 17,
> in swig_import_helper
> _mod = imp.load_module('_gdal', fp, pathname, description)
> ImportError: DLL load failed: %1 is not a valid Win32 application.
>
> F:\MODIS_Historical\First_season\yr2012\h27v07\Test>
>
>
> On Wed, Jun 4, 2014 at 3:32 PM, Saber Razmjooei <
> saber.razmjo...@lutraconsulting.co.uk> wrote:
>
>> Have you tried this:
>>
>> http://gis-lab.info/qa/extract-values-rasters-eng.html
>>
>>
>>
>> Cheers,
>>
>> Saber
>>
>>
>>
>>
>>
>> *From:* qgis-user-boun...@lists.osgeo.org [mailto:
>> qgis-user-boun...@lists.osgeo.org] *On Behalf Of *Leo Kris Palao
>> *Sent:* 04 June 2014 06:50
>> *To:* QGIS Mailing List
>> *Subject:* [Qgis-user] Share Python Script for Point Analyses Tool
>>
>>
>>
>> Hi QGIS Users,
>>
>> I am processing multi-temporal EVI data derived from MODIS, and the
>> number of datasets could be around 70 images. I have points and want to
>> extract pixel values from those points. I know that the Point Analyses Tool
>> can do it. However, it is cumbersome to load all the MODIS data because it
>> takes too long. Is there a way to extract the pixel values from points
>> without actually loading the datasets? Maybe a script (GDAL or Python) can
>> do this. Can somebody help me or point me to the right blog or website
>> where I can use as reference.
>&g

[Qgis-user] QgsZonalStatistics

[Leo Kris Palao]

Dear QGIS Users,

I am trying to run a script in python console to calculate zonal statistics
but it always throws an error. This is the script that I am running:

import qgis.analysis
vectorlayer=r'E:\filepath\vector.shp'
rasterfile = r'E:\filepath\raster.tif'
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(1), 1)
zonalstats.calculateStatistics(None)

The error says:

TypeError: arguments dit not match any overloaded call:
QgsZonalStatistics(QgsVectorLayer, QString, QString attributePrefix="", int
rasterBand=1): argument 1 has unexpected type 'str'
QgsZonalStatistics(QgsZonalStatistics): argument 1 has unexpected type 'str'

However, if I load my raster and vector file in QGIS and run the script
below, it is working properly:

import qgis.analysis
vectorlayer=qgis.utils.iface.mapCanvas().layer(0)
rasterfile = qgis.utils.iface.mapCanvas().layer(1).source()
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(1), 1)
zonalstats.calculateStatistics(None)

Thanks in advance for any help.
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] atcorr implementation in QGIS

[Leo Kris Palao]

Hi Users,

I want to run i.atcorr module in the processing toolbox but with no
success. The process throws the error/warning below.

I am processing Landsat 8 data. I've tried different options in the atcorr
interface - with DEM and no DEM; output is floating and not floating
number. But still gets the error/warning message below.

Other details:
Os - Windows 7 64 bit
version: QGIS 2.2 32 bit installed using OSGeo Installer

Thanks,
-Leo

Oooops! The following output layers could not be open

   - Name for output raster map:
   E:/Landsat_8_trial/LC81330482013106LGN01/Atcor/b1.tif

The above files could not be opened, which probably indicates that they
were not correctly produced by the executed algorithm

Checking the log information might help you see why those layers were not
created as expected

This algorithm requires GRASS to be run. A test to check if GRASS is
correctly installed and configured in your system has been performed, with
the following result:

GRASS seems to be correctly installed and configured
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] ZonalStats memory issue

[Leo Kris Palao]

Hi QGIS Users,

Is there memory limitation when using Zonal Stats? When I run it in QGIS
interface or Python console after a progress of 1% it throws this error
message then crash/close.

Error msge:
writing of minidump to
C:\Users\NAME\AppData\Local\Temp\qgis-20140709-174257-37188-37436-c3a2817.dmp
failed (80070008)

Thanks in advance for help/clarification.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Vector symbology from attribute table

[Leo Kris Palao]

Hi QGIS users,

Is there a way to automatically define the symbology of vector layer using
a column in the attribute table that defines the RGB?

I want to do this to have a uniform color coding for each land use types.
Some land use types might have more or less classes than for each
administrative boundary. I am thinking if I can encode the RGB values in my
attribute table and later use it to define my symbology.

Thanks for any help and advice,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Vector symbology from attribute table

[Leo Kris Palao]

Hi Zoltan,

I tried it and it was working, it defines well the symbology of my vector
layer. However, there is a little bit of problem. In my legend items it
only create one color (black). the color scheme of my legend doesn't seem
to follow the color combinations that I set in my attribute table. it would
be better if the color scheme of my attribute table will be reflected also
in my legend items.

This is how I encoded the color codes in my attribute table:
class  color (R, G, B)
1 197, 197, 197

Thanks,
-Leo


On Mon, Aug 4, 2014 at 5:31 PM, Siki Zoltan  wrote:

> Dear Leo,
>
> you can do it in the properties dialog of the layer.
> For example if you have a polygon layer, and the fill color is stored in
> the attribute table of the layer than
> 1. open properties dialog of the layer
> 2. select style tab (left)
> 3. click on "Simple fill" under symbol layers (if you can't see Simple
> fill, click on the triangle before Fill to extend the tree)
> 4. click on data defined properties... button
> 5. you can give an attribute or an expression of attributes to define the
> fill color
>
> I hope it helps you,
> Zoltan
>
>
> On Mon, 4 Aug 2014, Leo Kris Palao wrote:
>
>  Hi QGIS users,
>>
>> Is there a way to automatically define the symbology of vector layer using
>> a column in the attribute table that defines the RGB?
>>
>> I want to do this to have a uniform color coding for each land use types.
>> Some land use types might have more or less classes than for each
>> administrative boundary. I am thinking if I can encode the RGB values in
>> my
>> attribute table and later use it to define my symbology.
>>
>> Thanks for any help and advice,
>> -Leo
>>
>>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] QGIS on android and QT Forms

[Leo Kris Palao]

Hi QGIS users,

Good day!

Would like to ask if how can I implement QT custom forms in Android -
preferably in Samsung Galaxy Tab? we would like to use it for our
fieldwork. Is it possible to build the QT forms in PC and just copy the .ui
file in the android tablet?

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Python Path

[Leo Kris Palao]

Hi QGIS users,

Good Day!

Would like to ask how can I execute my python script in Python Shell? It
keeps on asking me that the module that I am using does not exist.

I converted the OSGeo4W shell to a python shell using the python command.
C:\> python

Then I execute my file >>>
execfile('folderpath\\multiband_zonalstat_pythonshell.py')

*Then it throws an error:*
*Traceback (most recent call last):*
*  File "", line 1, in *
*  File "folderpath\multiband_zonalstat_pythonshell.py", line 1, in *
*import qgis.analysis*

This is the script that I am trying to run in python shell:

import qgis.analysis
vectorlayer = QgsVectorLayer(r"folderpath/YEU.shp", "Ye-U", "ogr")
rasterfile = "rasterpath/cdz_raw_image.img"
for num in range(1):
zonalstats =
qgis.analysis.QgsZonalStatistics(vectorlayer,rasterfile,str(num),num)
 zonalstats.calculateStatistics(None)

I tried checking the system path (please see below), and I think it can
locate my OSGeo python installation path.

I also tried editing my system variables path and include the OSGeo python,
and OSGeo QGIS path but still it throws the same error.

I am trying to run my script in a python shell to test if the processing
time will be much shorter compare to running it in the python console in
QGIS.

On the other hand, I also tried to include the following lines in the
beginning of my script with no success:
*import sys*
*sys.path.extend([r'C:\OSGeo4W\apps\qgis\python',
r'C:\OSGeo4W\apps\qgis\bin'])*

it throws this error:

*Traceback (most recent call last):*
*  File "", line 1, in *
*  File "E:\Python_scripts\multiband_zonalstat_pythonshell.py", line 3, in
*
*import qgis.analysis*
*  File "C:\OSGeo4W\apps\qgis\python\qgis\__init__.py", line 35, in
*
*from qgis.core import QgsFeature, QgsGeometry*
*ImportError: DLL load failed: The specified module could not be found.*

I am not sure why it can't find qgis.analysis module or is there something
wrong with my script.

Thanks for the help and clarification,
-Leo

GDAL 1.10.1, released 2013/08/26

C:\>python
Python 2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]
on win
32
Type "help", "copyright", "credits" or "license" for more information.
>>> import sys
>>> sys.path
['', 'C:\\OSGeo4W\\apps\\orfeotoolbox\\wrapping\\lib',
'C:\\OSGeo4W\\apps\\orfeo
toolbox\\wrapping\\Python', 'C:\\OSGeo4W\\apps\\orfeotoolbox\\python',
'C:\\', '
C:\\OSGeo4W\\bin\\python27.zip', 'C:\\OSGeo4W\\apps\\Python27\\DLLs',
'C:\\OSGeo
4W\\apps\\Python27\\lib', 'C:\\OSGeo4W\\apps\\Python27\\lib\\plat-win',
'C:\\OSG
eo4W\\apps\\Python27\\lib\\lib-tk', 'C:\\OSGeo4W\\bin',
'C:\\OSGeo4W\\apps\\Pyth
on27', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages',
'C:\\OSGeo4W\\apps\\Py
thon27\\lib\\site-packages\\PIL',
'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packag
es\\win32', 'C:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32\\lib',
'C:\\
OSGeo4W\\apps\\Python27\\lib\\site-packages\\Pythonwin',
'C:\\OSGeo4W\\apps\\Pyt
hon27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win32.egg',
'C:\\OSGeo4W\\apps\\
Python27\\lib\\site-packages\\wx-2.8-msw-unicode']
>>> execfile('E:\\folderpath\\multiband_zonalstat_pytonshell.py')
Traceback (most recent call last):
  File "", line 1, in 
IOError: [Errno 2] No such file or directory: 'folderpath\\multiband_zon
alstat_pytonshell.py'
>>> execfile('fodlerpath\\multiband_zonalstat_pythonshell.py')
Traceback (most recent call last):
  File "", line 1, in 
  File "folderpath\multiband_zonalstat_pythonshell.py", line 1, in 
import qgis.analysis
ImportError: No module named qgis.analysis
>>>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] polygon centroids

[Leo Kris Palao]

Hi QGIS users,

I don't know if this question has been asked before.

Can you help me create a syntax in field calculator to add in the attribute
table the centroids coordinates for each or selected polygon features?

I am using Qgis 2.2 "Valmiera" version in Windows.

Thanks in advance for help,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] polygon centroids

[Leo Kris Palao]

Hi Alexander,

Thanks for your reply. In my understanding, this only applies if I have
point coordinates, right? If I have polygons I think I need another syntax?

I can create centroids in Vector > Geometry > Polygon Centroids tool then
apply $x and $y to extract lat lon coordinates, but in some cases I want to
skip this process if I am dealing with more than 300 thousand features in
my polygon layer.

Thanks,
-Leo


On Thu, Aug 28, 2014 at 3:01 PM, Alexandre Neto 
wrote:

> Hello Leo,
>
> That's quite straightforward actually. In the field calculator (or in any
> other instance of expression builder) there are 2 geometry functions to
> achieve what you need,  $x and $y.
>
> You can either put them in separate fields or create a string to populate
> a single one. Something like this,
>
> $x || ' , ' || $y
>
> Would produce something like this:
>
> '1234.56 , 9876.43'
>
> Notice that the double pipe character || works as a string concatenator.
> And that the result is a string therefore not useful for any arithmetic
> operations after.
>
> Hope this helped.
>
> Best regards,
>
> Alexandre Neto
> Em 27/08/2014 22:37, "Leo Kris Palao"  escreveu:
>
>> Hi QGIS users,
>>
>> I don't know if this question has been asked before.
>>
>> Can you help me create a syntax in field calculator to add in the
>> attribute table the centroids coordinates for each or selected polygon
>> features?
>>
>> I am using Qgis 2.2 "Valmiera" version in Windows.
>>
>> Thanks in advance for help,
>> -Leo
>>
>> ___
>> Qgis-user mailing list
>> Qgis-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>>
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Qgis-user Digest, Vol 103, Issue 13

[Leo Kris Palao]

Hi Paolo,

Is there a plan to make earthbeat plugin available in qgis 2.2 "Valmiera".
I tried to download it in the repository and copy it in the Plugins folder
in .qgis2\python\plugins. But in the Plugins window I cannot activate it.
The message say: *This plugin is incompatibe with this version of QGIS.
Plugin designed for QGIS 1.7.4 - 1.99*

Thanks,
-Leo


On Tue, Sep 9, 2014 at 3:00 AM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Re: scale not working in Qgis (Amedeo Fadini)
>2. Re: Error in grid in Print composer in Qgis Master (Nyall Dawson)
>3. unable to open Attribute Table (Jane Kogelschatz)
>4. Re: scale not working in Qgis (Mike Flannigan)
>5. Re: Spectral Profile (Paolo Cavallini)
>6. Re: scale not working in Qgis (Amedeo Fadini)
>7. Plugin to easilly filter several layers with a field  values
>   (Alexandre Neto)
>8. selecting features containing a polygon (Giacomo Fontanelli)
>9. Re: GRASS r.flow command error (john polo)
>
>
> --
>
> Message: 1
> Date: Sun, 7 Sep 2014 22:04:31 +0200
> From: Amedeo Fadini 
> To: Mike Flannigan 
> Cc: Qgis Users List - International 
> Subject: Re: [Qgis-user] scale not working in Qgis
> Message-ID:
> <
> cakgz3gc1qvwo0h-ob-eq3xb_c_3cjvksmo33_szpummxl4r...@mail.gmail.com>
> Content-Type: text/plain; charset=UTF-8
>
> Hi
>
> On Sat, Sep 6, 2014 at 10:09 PM, Mike Flannigan  wrote:
> >
> > Howdy,
> >
> > I assume you can zoom out (using the mouse wheel or any
> > other zoom method) and the scale changes??
> >
>
> the zoom changes but not the scale indicatore
>
> af
>
>
> --
>
> Message: 2
> Date: Mon, 8 Sep 2014 07:04:29 +1000
> From: Nyall Dawson 
> To: Gerardo Jim?nez 
> Cc: qgis-user 
> Subject: Re: [Qgis-user] Error in grid in Print composer in Qgis
> Master
> Message-ID:
>  ekhd1d1jzactxx1qnczhceh8lmdwyt-5p...@mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> On 06/09/2014 8:30 pm, "Gerardo Jim?nez"  wrote:
> >
> > Hi
> >  I am using Qgis 5.0 86987b5. I am using some layers in EPSG 26714 (UTM
> NAD27 14N). I tried to put a grid in EPSG 4267 (lat-long NAD27) in print
> composer and label it. I choosed the degree minute second with suffix
> option but  the labels do not place well the hemisphere indicator. All my
> longitudes are labeled "S" and my latitudes are labeled "E". I do not know
> if this is a bug or if there is an option to change this behavior.
> >
>
> Gerardo,
>
> Do you have a project you can share that demonstrates this? I'm trying to
> track it down and a sample project would help a lot. I don't need any data
> files, just the qgs project file.
>
> Regards,
> Nyall
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20140908/28dd1868/attachment-0001.html
> >
>
> --
>
> Message: 3
> Date: Sun, 7 Sep 2014 17:21:30 -0400
> From: "Jane Kogelschatz" 
> To: 
> Subject: [Qgis-user] unable to open Attribute Table
> Message-ID: <001701cfcae1$b1cce000$1566a000$@gmail.com>
> Content-Type: text/plain; charset="us-ascii"
>
> Hello,
>
> I am unable to open Attribute Table in QGIS Browser or Desktop 2.4.0
>
> Browser error:  "qgis-browser-bin.exe has stopped working"
>
> Just feature id appears
>
> Thank you for any suggestions.
>
> Computer is free of viruses and registry is clean.
>
>
>
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20140907/cd7f135f/attachment-0001.html
> >
>
> --
>
> Message: 4
> Date: Sun, 07 Sep 2014 16:46:24 -0500
> From: Mike Flannigan 
> To: Amedeo Fadini , Qgis Users List - International
> 
> Subject: Re: [Qgis-user] scale not working in Qgis
> Message-ID: <540cd230.1080...@att.net>
> Content-Type: text/plain; charset=UTF-8; format=flowed
>
>
> That is a strange one for sure.  Can you pick the small
> black triangle to the right of the scale and select a
> drop-down scale, like 1:500,000?
>
>
> Mike
>
>
> On 9/7/2014 3:04 PM, Amedeo Fadini wrote:
> > Hi
> >
> > the zoom changes but not the scale indicatore
> >
> > af
> >
>
>
>
> --
>
> Message: 5
> Date: Mon, 08 Sep 2014 09:10:19 +0200
> From: Paolo Cavallini 
> To: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] Spectr

[Qgis-user] Clip polygon by raster

[Leo Kris Palao]

Dear QGIS users,

I would like to ask how to clip a polygon using a raster layer?

basically, I have a raster layer and I want to use it as a clip layer to my
vector layer.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] GDAL utilities or GDAL Python book/ebook

[Leo Kris Palao]

Hi Guys,

Just would like to ask if you can suggest a good website or book on the
following topics:

1) Batch scripting using GDAL utilities

2) Book (ebook or hardcopy) on mapping using GDAL-Python with remote
sensing, and raster analysis topics.

I am not a hardcore Linux OS user. I would like it if the examples are in
Windows machine.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Generate Random Points in QGIS Processing Toolbox

[Leo Kris Palao]

Hi ALL,

Good day!

I am trying to build a model from the Graphical Modeler but cannot find the
particular tool to generate random points like the one in vector > research
tools > *random points* but cannot get/find the equivalent function in the
processing toolbox.

In one of the steps in my model, I would like to generate a specified
number of points per ID in my shapefile layer.

Can this tool be used in the graphical modeler?

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Generate Random Points in QGIS Processing Toolbox

[Leo Kris Palao]

Hi Victor,

under the Processing toolbox > QGIS geoalgorithms > Vector create tools, I
only have Create grid, and Points layer from table. I am using QGIS 2.6.1.

Thanks,
-Leo

On Wed, Feb 4, 2015 at 6:18 PM, Victor Olaya  wrote:

> You should have at least the following algorithms (all under the QGIS
> algorithms/vector creation tools group):
>
> Random points in extent
> Random points in layer bounds
> random points in polygons(fixed)
> random points in polygons(variable)
>
>
> Hope this helps
>
> 2015-02-04 8:16 GMT+01:00 Alexander Bruy :
> > Hi Leo,
> >
> > which version of QGIS and Processing do you use? Tools for
> > generating random points were added some time ago.
> > Did you tried to use filter field in Processing toolbox to find
> > algorithms?
> >
> > 2015-02-04 3:17 GMT+02:00 Leo Kris Palao :
> >> Hi ALL,
> >>
> >> Good day!
> >>
> >> I am trying to build a model from the Graphical Modeler but cannot find
> the
> >> particular tool to generate random points like the one in vector >
> research
> >> tools > random points but cannot get/find the equivalent function in the
> >> processing toolbox.
> >>
> >> In one of the steps in my model, I would like to generate a specified
> number
> >> of points per ID in my shapefile layer.
> >>
> >> Can this tool be used in the graphical modeler?
> >>
> >> Thanks,
> >> -Leo
> >>
> >> ___
> >> Qgis-user mailing list
> >> Qgis-user@lists.osgeo.org
> >> http://lists.osgeo.org/mailman/listinfo/qgis-user
> >
> >
> >
> > --
> > Alexander Bruy
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Date and time format

[Leo Kris Palao]

Hi QGIS users,

I have a column Date -MM-DD and Time HH:MM:SS in QString format. I
would like to convert this to Date Time format -MM-DD HH:MM:SS.

We would like to explore using the date time values into the interpolation
model.

Any help will be much appreciated.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] "Date" and "Time" Column to "Date Time" column

[Leo Kris Palao]

Hi ALL,

I have two columns, namely "Date" and "Time" in QString format. How can I
convert these columns to Date Time (-mm-dd hh:mm:ss) format.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] batch conversion of raster file to GeoTIFF format

[Leo Kris Palao]

Hi QGIS users,

I am having a problem converting my *539 raster* files to *GeoTIFF format*.
I tried using the GDAL *Translate* option in the *Processing Toolbox* but
it reads only the **.hdr* file of my raster not the file that holds my
raster data. When I execute the batch process file using *file.hdr* as my
input it returns:

Processing of the following files ended with error:
> ERROR 1: The selected file is an ENVI header file, but to
> open ENVI datasets, the data file should be selected
> instead of the .hdr file. Please try again selecting
> the data file corresponding to the header file:
> J:/MODIS_EVI_smooth/mod_2002345_EVI_smth.hdr

>> J:/MODIS_EVI_smooth/mod_2002345_EVI_smth.hdr
>

I cannot find the appropriate file type option when selecting which files
to load as my inputs. Even thought I select* files(*.)* type option, still
I cannot find my files in my directory.

I have also tried to load the *file.hdr* and later remove the **.hdr*
extension before executing the batch process but with no success.

This is the sample file that I am working on:

mod_2002345_EVI_smthFile...45,000KB
> mod_2002345_EVI_smth.hdr.. HDR.1KB

mod_2002345_EVI_smth.sml..SML File3KB


I have successfully converting one raster file using Translate module
under *Raster
> Conversion > Translate* but I can only process one at a time. But I don't
want to do this for my remaining raster data. Please also note the I have
tried the batch mode but with no success.

I copied the GDAL script (below) and I wonder if I can run a batch
processing of this script in OSGeo Shell. Just my guess.

gdal_translate -a_srs "+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181
> +b=6371007.181 +units=m +no_defs" -a_nodata -3000 -of GTiff
> J:/MODIS_EVI_smooth/mod_2002345_EVI_smth
> J:/MODIS_EVI_smooth_TIFF/mod_2002345_EVI_smth.tif


Can you help me with this.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] QGIS Export Map

[Leo Kris Palao]

Hi QGIS users,

I have a bit of a problem. I am using Open Layers Google Physical Layer as
background of my other layers. When I layout my maps in Print Composer, the
layers are well aligned. But when I export it to JPEG or PNG the output
product is misaligned. I allow on-the-fly option when the layer from Open
Layers were added. So, I thought this might be causing the issue. But then
I transformed all of my layers to EPSG 3857 but still problem persist.

Can you help me with this?

I am using QGIS 2.6.1 Brighton installed via OSGeo.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS Export Map

[Leo Kris Palao]

Hi Nyall,

Thanks for the quick answer.

-Leo

On Mon, Mar 9, 2015 at 11:02 AM, Nyall Dawson 
wrote:

> On 9 March 2015 at 13:56, Leo Kris Palao  wrote:
> > Hi QGIS users,
> >
> > I have a bit of a problem. I am using Open Layers Google Physical Layer
> as
> > background of my other layers. When I layout my maps in Print Composer,
> the
> > layers are well aligned. But when I export it to JPEG or PNG the output
> > product is misaligned. I allow on-the-fly option when the layer from Open
> > Layers were added. So, I thought this might be causing the issue. But
> then I
> > transformed all of my layers to EPSG 3857 but still problem persist.
> >
> > Can you help me with this?
>
> Short answer is that the openlayers plugin doesn't support the map
> composer. There's various hacks which reportedly work for some people
> and not others, but I wouldn't recommend wasting time trying to
> explore these.
>
> Nyall
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] Zonal Statistics Computation

[Leo Kris Palao]

Hi QGIS Users,

I would like to ask if anybody have experienced getting different results
using *A)* Zonal Statistics Plugin and *B)* Zonal Statistics in Graphical
Modeler?

Here is the description of our *Dataset*:

*Vector*
5 features (2 features have area of 1.74 and 1.89 ha; 3 features have less
than 1 ha. area). These are farmer plots.

*Raster*
pixel size is 180m X 180m = 3.24ha.

*Processing*

Running Zonal Statistics in Graphical Modeler:
Only computes realistic statistics for 2 features (area having 1.74 and
1.89 ha) but gives weird computation for features having less than 1 ha.
For features less than 1 ha, it gives values 0. and 1.

Running Zonal Statistics in Raster > Zonal Statistics > Zonal Statistics:
It gives realistic values for all features.

This happens in QGIS "Brighton" and "Wien".

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Qgis-user Digest, Vol 109, Issue 23

[Leo Kris Palao]

HI Stefan,

Yes, I read the issue in the link that you gave me. In one thread, there is
an issue of both zonal statistics plugin and zonal statistics in processing
toolbox are not working. But in my case the ZonalStats Plugin is working
and the other is not. I browsed my data again and found that some polygons
are not in the middle of the raster. However, it seems that ZonalStats
plugin considers the polygons even though other the polygons are not in the
center of the raster. Is there a difference between the algorithm between
the two ZonalStats option. By the way, I am using the ZonalStatistics
under *QGIS
 geoalgorithms* > *Raster tools*.

Is there something that I missed from the threads that you shared with me?

Thanks for the help.

-Leo

On Fri, Mar 13, 2015 at 3:01 AM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Re: Customised attributes form (emmexx)
>2. Re: Customised attributes form (Matthias Kuhn)
>3. Re: Zonal Statistics Computation (Blumentrath, Stefan)
>4. Attribute copy from point to polygon (Joris Hintjens)
>5. Re: Customised attributes form (emmexx)
>6. Re: Attribute copy from point to polygon (Raymond Nijssen)
>7. Cannot activate grass (Lu?s de Sousa)
>8. Re: Cannot activate grass (Carlos Cerd?n)
>9. Re: Problem loading ArcGIS georeferenced map in QGIS 2.8
>   (Nicole Stoffels)
>   10. Re: GRASS dependency problem on OS X (William Kyngesburye)
>   11.  Cannot activate grass (Lu?s de Sousa)
>   12. Re: Cannot activate grass (Carlos Cerd?n)
>   13. TOA Reflectance (Mariana Buzzi)
>
>
> --
>
> Message: 1
> Date: Thu, 12 Mar 2015 11:21:07 +0100
> From: emmexx 
> To: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] Customised attributes form
> Message-ID: <55016893.3070...@tiscalinet.it>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Il 03/12/2015 10:46 AM, emmexx scrisse:
> > I want to use a custom interface - I need to organize attributes better
> > than the qgis form can - but I don't understand if, and how, I can mix a
> > custom interface, custom widgets and custo validation. Is that possible?
>
> I overlooked the Drag and drop designer option. Is that the right path
> to accomplish what I want?
>
> thanks
> maxx
>
>
> --
>
> Message: 2
> Date: Thu, 12 Mar 2015 11:32:02 +0100
> From: Matthias Kuhn 
> To: emmexx 
> Cc: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] Customised attributes form
> Message-ID: <55016b22.1050...@gmx.ch>
> Content-Type: text/plain; charset=UTF-8
>
> The d'n'd designer is about grouping, positioning and layouting
> different fields on a form.
>
> If I understand you correctly, you have 1 single attribute that contains
> a list of values.
> For that, the linked example is the way to go.
>
> In "setValue" you take 1 value (which is a list of values) and highlight
> the appropriate ones on your tree
>
> In "value" you inspect the (user modified) tree and generate a value
> from it to save back to the attribute.
>
> This tree widget can then be positioned with the d'n'd designer (or a
> QtDesigner form if you like) together with other attributes this layer has.
>
> Matthias
>
> On 03/12/2015 11:21 AM, emmexx wrote:
> > Il 03/12/2015 10:46 AM, emmexx scrisse:
> >> I want to use a custom interface - I need to organize attributes better
> >> than the qgis form can - but I don't understand if, and how, I can mix a
> >> custom interface, custom widgets and custo validation. Is that possible?
> >
> > I overlooked the Drag and drop designer option. Is that the right path
> > to accomplish what I want?
> >
> > thanks
> > maxx
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
>
>
>
> ----------
>
> Message: 3
> Date: Thu, 12 Mar 2015 10:42:55 +
> From: "Blumentrath, Stefan" 
> To: Leo Kris P

[Qgis-user] QGIS tracking/fieldwork application

[Leo Kris Palao]

Hi QGIS Users,

Good day!

I watched the live tracking video of QGIS in
https://www.youtube.com/watch?v=0WevRW4tbzs&t=253
(Live Tracking and sampling with Quantum GIS, a tablet PC and a Bluetooth
datalogger)

I am interested to the same. I am particularly interested on how I can
customize my forms such that when I digitize a point, I can use drop down
menu options, combo box, etc. If you could point me in the right direction
(blogs, websites), will be very helpful.

Thanks,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS tracking/fieldwork application

[Leo Kris Palao]

Thanks Bo,

I will have a look on the Roam QGIS. Can you also suggest some blogs or
websites using Android tablet?

Thanks,
-Leo

On Tue, Mar 17, 2015 at 3:01 AM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Re: Line symbols with edge along the line and without edge at
>   the ends? (Andreas Neumann)
>2. Re: QGIS tracking/fieldwork application (Bo Victor Thomsen)
>3. Oracle connection Timeout (Hillers Heather)
>4. Re: Oracle connection Timeout (J?rgen E. Fischer)
>5. Can't access Oracle sequences (PIERRE Sylvain)
>6. Add New Records to Relation or Non-Spatial Table Using "Line
>   Entry" (vs "Form Entry") (Jeffrey Durrence)
>7. Building an input form with choice-dependent dropdown lists
>   (Jeffrey Durrence)
>
>
> --
>
> Message: 1
> Date: Mon, 16 Mar 2015 11:18:11 +0100
> From: Andreas Neumann 
> To: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] Line symbols with edge along the line and
> without edge at the ends?
> Message-ID: <5506ade3.4050...@carto.net>
> Content-Type: text/plain; charset="windows-1252"; Format="flowed"
>
> Hi,
>
> This is possible. Make sure that you enable symbol levels (button
> "Advanced -> Symbol levels at the lower right of the symbology dialogue)
> and set the line endings to "flat" (not square or round).
>
> Andreas
>
> On 16.03.2015 10:47, Elstermann, Mike wrote:
> >
> > Line symbols with edge along the line and without edge at the ends?
> >
> > Left NOT OK, Right OK
> >
> > THX, mikeE.
> >
> >
> >
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
>
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20150316/5991d658/attachment-0001.html
> >
> -- next part --
> A non-text attachment was scrubbed...
> Name: not available
> Type: image/jpeg
> Size: 5997 bytes
> Desc: not available
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20150316/5991d658/attachment-0001.jpe
> >
>
> --
>
> Message: 2
> Date: Mon, 16 Mar 2015 11:25:57 +0100
> From: Bo Victor Thomsen 
> To: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] QGIS tracking/fieldwork application
> Message-ID: <5506afb5.9030...@gmail.com>
> Content-Type: text/plain; charset="windows-1252"; Format="flowed"
>
> Try googling "Roam QGIS". Roam is a new "skin" aka user interface (and
> more) for QGIS, enabling it to function better on a windows tablet.
>
> Regards
> Bo Victor Thomsen
> AestasGIS
> Denmark
>
>
>
> Den 16-03-2015 kl. 08:54 skrev Leo Kris Palao:
> > Hi QGIS Users,
> >
> > Good day!
> >
> > I watched the live tracking video of QGIS in
> > https://www.youtube.com/watch?v=0WevRW4tbzs&t=253
> >
> >
> >   (Live Tracking and sampling with Quantum GIS, a tablet PC and a
> >   Bluetooth datalogger)
> >
> >
> > I am interested to the same. I am particularly interested on how I can
> > customize my forms such that when I digitize a point, I can use drop
> > down menu options, combo box, etc. If you could point me in the right
> > direction (blogs, websites), will be very helpful.
> >
> > Thanks,
> > -Leo
> >
> >
> > ___
> > Qgis-user mailing list
> > Qgis-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/qgis-user
>
> -- next part --
> An HTML attachment was scrubbed...
> URL: <
> http://lists.osgeo.org/pipermail/qgis-user/attachments/20150316/c4e1093b/attachment-0001.html
> >
>
> --
>
> Message: 3
> Date: Mon, 16 Mar 2015 14:12:23 +0100
> From: Hi

Re: [Qgis-user] Zonal Statistics Computation

[Leo Kris Palao]

Hi Stefan,

Yes, I read the issue in the link that you gave me. In one thread, there is
an issue of both zonal statistics plugin and zonal statistics in processing
toolbox are not working. But in my case the ZonalStats Plugin is working
and the other is not. I browsed my data again and found that some polygons
are not in the middle of the raster. However, it seems that ZonalStats
plugin considers the polygons even though other the polygons are not in the
center of the raster. Is there a difference between the algorithm between
the two ZonalStats option. By the way, I am using the ZonalStatistics
under *QGIS
 geoalgorithms* > *Raster tools*.

Is there something that I missed from the threads that you shared with me?

Thanks for the help.

-Leo

On Fri, Mar 13, 2015 at 8:38 AM, Leo Kris Palao  wrote:

> HI Stefan,
>
> Yes, I read the issue in the link that you gave me. In one thread, there
> is an issue of both zonal statistics plugin and zonal statistics in
> processing toolbox are not working. But in my case the ZonalStats Plugin is
> working and the other is not. I browsed my data again and found that some
> polygons are not in the middle of the raster. However, it seems that
> ZonalStats plugin considers the polygons even though other the polygons are
> not in the center of the raster. Is there a difference between the
> algorithm between the two ZonalStats option. By the way, I am using the
> ZonalStatistics under *QGIS  geoalgorithms* > *Raster tools*.
>
> Is there something that I missed from the threads that you shared with me?
>
> Thanks for the help.
>
> -Leo
>
> On Fri, Mar 13, 2015 at 3:01 AM, 
> wrote:
>
>> Send Qgis-user mailing list submissions to
>> qgis-user@lists.osgeo.org
>>
>> To subscribe or unsubscribe via the World Wide Web, visit
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>> or, via email, send a message with subject or body 'help' to
>> qgis-user-requ...@lists.osgeo.org
>>
>> You can reach the person managing the list at
>> qgis-user-ow...@lists.osgeo.org
>>
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of Qgis-user digest..."
>>
>>
>> Today's Topics:
>>
>>1. Re: Customised attributes form (emmexx)
>>2. Re: Customised attributes form (Matthias Kuhn)
>>3. Re: Zonal Statistics Computation (Blumentrath, Stefan)
>>4. Attribute copy from point to polygon (Joris Hintjens)
>>5. Re: Customised attributes form (emmexx)
>>6. Re: Attribute copy from point to polygon (Raymond Nijssen)
>>7. Cannot activate grass (Lu?s de Sousa)
>>8. Re: Cannot activate grass (Carlos Cerd?n)
>>9. Re: Problem loading ArcGIS georeferenced map in QGIS 2.8
>>   (Nicole Stoffels)
>>   10. Re: GRASS dependency problem on OS X (William Kyngesburye)
>>   11.  Cannot activate grass (Lu?s de Sousa)
>>   12. Re: Cannot activate grass (Carlos Cerd?n)
>>   13. TOA Reflectance (Mariana Buzzi)
>>
>>
>> --
>>
>> Message: 1
>> Date: Thu, 12 Mar 2015 11:21:07 +0100
>> From: emmexx 
>> To: qgis-user@lists.osgeo.org
>> Subject: Re: [Qgis-user] Customised attributes form
>> Message-ID: <55016893.3070...@tiscalinet.it>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>
>> Il 03/12/2015 10:46 AM, emmexx scrisse:
>> > I want to use a custom interface - I need to organize attributes better
>> > than the qgis form can - but I don't understand if, and how, I can mix a
>> > custom interface, custom widgets and custo validation. Is that possible?
>>
>> I overlooked the Drag and drop designer option. Is that the right path
>> to accomplish what I want?
>>
>> thanks
>> maxx
>>
>>
>> --
>>
>> Message: 2
>> Date: Thu, 12 Mar 2015 11:32:02 +0100
>> From: Matthias Kuhn 
>> To: emmexx 
>> Cc: qgis-user@lists.osgeo.org
>> Subject: Re: [Qgis-user] Customised attributes form
>> Message-ID: <55016b22.1050...@gmx.ch>
>> Content-Type: text/plain; charset=UTF-8
>>
>> The d'n'd designer is about grouping, positioning and layouting
>> different fields on a form.
>>
>> If I understand you correctly, you have 1 single attribute that contains
>> a list of values.
>> For that, the linked example is the way to go.
>>
>> In "setValue" you take 1 value (which is a list of values) and highlight
>> the appropriate ones on your tree
>>
>

[Qgis-user] Snap Raster

[Leo Kris Palao]

Hi QGIS Users,

Would like to inquire if there is an option in GDAL or in QGIS GUI to snap
cells of the other raster to another raster layer.

Reason for this is that, I have a raster projected in WGS 84 then I would
like to project it to MODIS Sinusoidal projection. But I want the output
reprojected layer to adjust and match the cell alignment of the
specified/reference raster.

Similar to Snap Raster option in ArcGIS Environment Settings.

Thanks in advance for the help.

-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS: great at sea...(teaching & research in the field)

[Leo Kris Palao]

Hi QGIS Users,

I have also a very nice experience with QGIS, that I would like to write a
simple case study for sharing. Last November, we used QGIS live tracking
capability to collect precise ground-truth data using medium resolution
imagery for complex land use pattern mapping. We simply overlay MODIS grids
(vector format) and landsat imagery in QGIS. By overlaying the grids in
QGIS we know our relative position on the ground with respect to each MODIS
pixels. Then we orient our eyes based on the extent and direction of the
pixels. We can accurately categorize the land use pattern for each pixel.
And this helps us a lot, since we might be collecting/evaluating ground
information for wrong pixels. We then use Add feature function to
automatically grab the current point. Thanks QGIS!

-Leo

On Mon, Mar 23, 2015 at 3:01 AM,  wrote:

> Send Qgis-user mailing list submissions to
> qgis-user@lists.osgeo.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://lists.osgeo.org/mailman/listinfo/qgis-user
> or, via email, send a message with subject or body 'help' to
> qgis-user-requ...@lists.osgeo.org
>
> You can reach the person managing the list at
> qgis-user-ow...@lists.osgeo.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Qgis-user digest..."
>
>
> Today's Topics:
>
>1. Re: QGIS: great at sea...(teaching & research in the field)
>   (Piotr Kania)
>
>
> --
>
> Message: 1
> Date: Sun, 22 Mar 2015 12:12:42 +0100
> From: Piotr Kania 
> To: qgis-user@lists.osgeo.org
> Subject: Re: [Qgis-user] QGIS: great at sea...(teaching & research in
> the field)
> Message-ID: <550ea3aa.4090...@op.pl>
> Content-Type: text/plain; charset=windows-1252; format=flowed
>
> My two cents - company I work in pays huge money for Arcgis, but Arcmap
> doesn't work with wmts. That's why I encourage everyone to use QGIS and
> they are surprised how good it is - for free they can use valuable
> public wmts data, what is impossible in Arcgis 10.2.2, 10.3
>
> W dniu 2015-03-21 o 10:18, Paolo Cavallini pisze:
> > Il 20/03/2015 23:24, Falk Huettmann ha scritto:
> >
> >> In summary, I can only say the best how this worked out on the ship and
> >> with students!
> > ...
> >> So in conclusion, I can recommend it highly for (ocean) field schools
> >> and field lab settings.
> >> Over time, I would be delighted to hear from other people and their
> >> experiences re. QGIS use and teaching etc, too.
> > Glad to see you enjoy our work.
> > Have you thought of joining the community and helping us? There are lot
> > of areas where some help would be important for the progress of QGIS.
> > Also, please consider writing a short case study about your experience.
> > All the best, and thanks.
> >
>
>
>
> --
>
> ___
> Qgis-user mailing list
> Qgis-user@lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/qgis-user
>
> End of Qgis-user Digest, Vol 109, Issue 45
> **
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] QGIS 2.8 runtime error

[Leo Kris Palao]

Hi QGIS Users,

I have this problem when opening QGIS. It reports:

Runtime error:
> Program C:\OSGeo4W\bin\qgis-bin.exe
> R6034
> An application has made an attempt to load C runtime library incorrectly.
> Please contact the application's support team for more information.


After clicking "OK" it loads QGIS. I have not tried to process anything so
I have no idea if this will affect some of the functions in QGIS.

I am using Windows 7 64 bit. Installed QGIS via OSGeo 32bit.

note: I already restarted my laptop.

Thanks in advance for any help,
-Leo
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS 2.8 runtime error

[Leo Kris Palao]

HI Ross and QGIS Users,

I tried the solution from the forum but still I am receiving runtime error.
I run process explorer and found that msvcr71.dll and msvcr100.dll in my
C:\OSGeo4W/bin\ directory. I tried to 1) remove both; 2) remove
msvcr71.dll, and 3) msvcr100.dll but still I am getting the runtime error
prompt when opening QGIS.

When I installed QGIS 2.8 I erased first all the qgis files in my C:\
directory before running the OSGeo installer. I do not know if this has
something to do with this issue. Because before I installed QGIS 2.6,
installed also using OSGeo installer, but I have no issue like this.

Thanks,
-Leo

On Tue, Mar 31, 2015 at 3:56 PM, McDonaldR  wrote:

>  Hi Leo
>
>
>
> I had the same error a while back and I posted this to the UK QGIS user
> group:
>
>
>
> Anyone getting the following error on the latest update to QGIS via
> OSGeo4W? The problem is a copy of msvcr90.dll in the
> C:\OSGeo4W\apps\qgis\bin folder which clashes with the copy in
> C:\Windows\winsxs\ folder.  Removing the DLL from the QGIS folder allows
> the program to start without a hitch.  If you have QGIS DEV installed there
> is no issue.
>
> This post provided the answer and method:
> http://stackoverflow.com/questions/14552348/runtime-error-r6034-in-embedded-python-application/14680947#14680947
>
> This happened after using the OSGeo4W installer to upgrade to QGIS 2.6.1-2
>
>
>
> Whether this is just my install environment or affects everyone, I don’t
> know.
>
>
>
> Hope it helps
>
>
>
> Ross
>
>
>
> *From:* qgis-user-boun...@lists.osgeo.org [mailto:
> qgis-user-boun...@lists.osgeo.org] *On Behalf Of *Leo Kris Palao
> *Sent:* 31 March 2015 03:05
> *To:* QGIS Mailing List
> *Subject:* [Qgis-user] QGIS 2.8 runtime error
>
>
>
> Hi QGIS Users,
>
>
>
> I have this problem when opening QGIS. It reports:
>
>
>
> Runtime error:
> Program C:\OSGeo4W\bin\qgis-bin.exe
> R6034
> An application has made an attempt to load C runtime library incorrectly.
> Please contact the application's support team for more information.
>
>
>
> After clicking "OK" it loads QGIS. I have not tried to process anything so
> I have no idea if this will affect some of the functions in QGIS.
>
>
>
> I am using Windows 7 64 bit. Installed QGIS via OSGeo 32bit.
>
>
>
> note: I already restarted my laptop.
>
>
>
> Thanks in advance for any help,
>
> -Leo
>
>
> This message is strictly confidential. If you have received this in error,
> please inform the sender and remove it from your system. If received in
> error you may not copy, print, forward or use it or any attachment in any
> way. This message is not capable of creating a legal contract or a binding
> representation and does not represent the views of Angus Council. Emails
> may be monitored for security and network management reasons. Messages
> containing inappropriate content may be intercepted. Angus Council does not
> accept any liability for any harm that may be caused to the recipient
> system or data on it by this message or any attachment.
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] QGIS 2.8 runtime error

[Leo Kris Palao]

Hi Ross and QGIS users,

I updated to QGIS 2.8.2 (thanks QGIS) but still I got the runtime error
R6034. I know this the issue is beyond QGIS, but can you give me some light
on how the procedure below works (if this works for most users with the
same issue):

Basically, I have the process explorer already and found out that there are
msvcr??.dll that is not in my windows system folder. msvcr71.dll and
msvcr100.dll located in path:\OSGeo4W\bin\ directory.

Step 6: Check the path before the application run. -- So I open QGIS and
the message is Runtime Error! Program: path:\OSGeo4W\bin\qgis-bin.exe. Is
this the path step 6 is referring to?

The solution is to remove the offending entry from the path using a batch
file. -- this is not quite clear to me how I will implement this one.

Thanks,
-Leo


   1.

   Download Process Explorer here
   <http://technet.microsoft.com/en-ca/sysinternals/bb896653.aspx>.
   2.

   Start your application and reproduce runtime error R6034.
   3.

   Start Process Explorer. In the "View" menu go to "Lower Pane View" and
   choose "DLLs".
   4.

   In the top pane, locate your application and click on it. The bottom
   pane should show a list of DLLS loaded for your application.
   5.

   Locate "msvcr??.dll" in the list. There should be several. Look for the
   one that is not in the "winsxs" folder, and make a note of it.
   6.

   Now, check the path just before your application runs. If it includes
   the folder you noted in step 5, you've probably found the culprit.

How to fix the problem? You'll have to remove the offending entry from the
path before running your program. In my case, I don't need anything else in
the path, so I wrote a simple batch file that looks like this:

path=
myprogram.exe

That's it. The batch file simply clears the path before my program runs, so
that the conflicting runtime DLL is not found.

Hope this helps!

On Wed, Apr 8, 2015 at 9:51 AM, Leo Kris Palao  wrote:

> HI Ross and QGIS Users,
>
> I tried the solution from the forum but still I am receiving runtime
> error. I run process explorer and found that msvcr71.dll and msvcr100.dll
> in my C:\OSGeo4W/bin\ directory. I tried to 1) remove both; 2) remove
> msvcr71.dll, and 3) msvcr100.dll but still I am getting the runtime error
> prompt when opening QGIS.
>
> When I installed QGIS 2.8 I erased first all the qgis files in my C:\
> directory before running the OSGeo installer. I do not know if this has
> something to do with this issue. Because before I installed QGIS 2.6,
> installed also using OSGeo installer, but I have no issue like this.
>
> Thanks,
> -Leo
>
> On Tue, Mar 31, 2015 at 3:56 PM, McDonaldR  wrote:
>
>>  Hi Leo
>>
>>
>>
>> I had the same error a while back and I posted this to the UK QGIS user
>> group:
>>
>>
>>
>> Anyone getting the following error on the latest update to QGIS via
>> OSGeo4W? The problem is a copy of msvcr90.dll in the
>> C:\OSGeo4W\apps\qgis\bin folder which clashes with the copy in
>> C:\Windows\winsxs\ folder.  Removing the DLL from the QGIS folder allows
>> the program to start without a hitch.  If you have QGIS DEV installed there
>> is no issue.
>>
>> This post provided the answer and method:
>> http://stackoverflow.com/questions/14552348/runtime-error-r6034-in-embedded-python-application/14680947#14680947
>>
>> This happened after using the OSGeo4W installer to upgrade to QGIS 2.6.1-2
>>
>>
>>
>> Whether this is just my install environment or affects everyone, I don’t
>> know.
>>
>>
>>
>> Hope it helps
>>
>>
>>
>> Ross
>>
>>
>>
>> *From:* qgis-user-boun...@lists.osgeo.org [mailto:
>> qgis-user-boun...@lists.osgeo.org] *On Behalf Of *Leo Kris Palao
>> *Sent:* 31 March 2015 03:05
>> *To:* QGIS Mailing List
>> *Subject:* [Qgis-user] QGIS 2.8 runtime error
>>
>>
>>
>> Hi QGIS Users,
>>
>>
>>
>> I have this problem when opening QGIS. It reports:
>>
>>
>>
>> Runtime error:
>> Program C:\OSGeo4W\bin\qgis-bin.exe
>> R6034
>> An application has made an attempt to load C runtime library incorrectly.
>> Please contact the application's support team for more information.
>>
>>
>>
>> After clicking "OK" it loads QGIS. I have not tried to process anything
>> so I have no idea if this will affect some of the functions in QGIS.
>>
>>
>>
>> I am using Windows 7 64 bit. Installed QGIS via OSGeo 32bit.
>>
>>
>>
>> note: I already restarted my laptop.
>>
>>
>>
>> Thanks in advance for any help,
&

[Qgis-user] Relocator Plugin

[Leo Kris Palao]

Hi QGIS Users,

I am trying to use the relocator plugin but it throws the error below.

I have 46 rasters open in QGIS in which I made symbology. I want to share
the project with my colleague and hoping to use the Relocator plugin but it
crashes.

Help please.

I am using QGIS 2.8.2, Windows 7, installed via OSGeo4W

Thanks,
-Leo

An error has occured while executing Python code:
> Traceback (most recent call last):
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\relocator\relocator_dialog.py", line
> 83, in relocator
> relocator_exec(path_r, zipit)
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\relocator\relocator_exec.py", line
> 72, in relocator_exec
> print file1 + " is a match!"
> IOError: [Errno 9] Bad file descriptor
> Python version:
> 2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]
>
> QGIS version:
> 2.8.2-Wien Wien, 1b929ef
> Python path: ['E:/OSGeo4W/apps/qgis/./python/plugins\\processing',
> 'E:/OSGeo4W/apps/qgis/./python', u'C:/Users/lpalao/.qgis2/python',
> u'C:/Users/lpalao/.qgis2/python/plugins',
> 'E:/OSGeo4W/apps/qgis/./python/plugins',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\setuptools-0.9.8-py2.7.egg',
> 'E:\\OSGeo4W\\apps\\orfeotoolbox\\wrapping\\lib',
> 'E:\\OSGeo4W\\apps\\orfeotoolbox\\wrapping\\Python',
> 'E:\\OSGeo4W\\apps\\orfeotoolbox\\python', 'C:\\OSGeo4W\\apps\\Python27',
> 'C:\\OSGeo4W\\apps\\qgis\\bin', 'C:\\OSGeo4W\\apps\\qgis\\python',
> 'C:\\Python27\\Python2.7.9', 'C:\\Python27\\Python2.7.9\\Scripts',
> 'E:\\OSGeo4W\\bin\\python27.zip', 'E:\\OSGeo4W\\apps\\Python27\\DLLs',
> 'E:\\OSGeo4W\\apps\\Python27\\lib',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\plat-win',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\lib-tk', 'E:\\OSGeo4W\\bin',
> 'E:\\OSGeo4W\\apps\\Python27',
> 'e:\\osgeo4w\\apps\\python27\\lib\\site-packages\\jinja2-2.7.2-py2.7.egg',
> 'e:\\osgeo4w\\apps\\python27\\lib\\site-packages\\markupsafe-0.23-py2.7-win32.egg',
> 'e:\\osgeo4w\\apps\\python27\\lib\\site-packages\\python_dateutil-2.2-py2.7.egg',
> 'e:\\osgeo4w\\apps\\python27\\lib\\site-packages\\pytz-2014.2-py2.7.egg',
> 'e:\\osgeo4w\\apps\\python27\\lib\\site-packages\\six-1.6.1-py2.7.egg',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\PIL',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\win32\\lib',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\Pythonwin',
> 'E:\\OSGeo4W\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode',
> u'C:/Users/lpalao/.qgis2//python',
> 'E:\\OSGeo4W\\apps\\qgis\\python\\plugins\\fTools\\tools', 'E:/Spatial
> Data/GIS_projects']
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

[Qgis-user] OpenLayers issue QGIS 2.10

[Leo Kris Palao]

Hi Guys,

I want to use OpenLayers plugin in my QGIS 2.10 Pisa but getting an error
message. This error started when I upgraded my QGIS to version 2.8 Wien,
and until now, in my QGIS 2.10 still have this. Can I have some advise how
to solve this issue?

I am using Windows 7 OS.

Thanks in advance,
-Leo

An error has occured while executing Python code:

Traceback (most recent call last):
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\weblayers\weblayer.py",
> line 97, in addLayer
> self._addLayerCallback(self)
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\openlayers_plugin.py",
> line 130, in addLayer
> layer = OpenlayersLayer(self.iface, self._olLayerTypeRegistry)
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py",
> line 330, in __init__
> self.olWebPage = OLWebPage(self)
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py",
> line 48, in __init__
> proxy = getProxy()
>   File
> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\tools_network.py",
> line 42, in getProxy
> proxy.setPort(settings.value("/proxyPort", 10, type=int))
> TypeError: unable to convert a QVariant of type 10 to a QMetaType of type 2
> Python version:
> 2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]
>
> QGIS version:
> 2.10.1-Pisa Pisa, d20c5b7
> Python path:
> ['C:/PROGRA~1/QGISPI~1/apps/qgis/./python/plugins\\processing',
> 'C:/PROGRA~1/QGISPI~1/apps/qgis/./python',
> u'C:/Users/lpalao/.qgis2/python', u'C:/Users/lpalao/.qgis2/python/plugins',
> 'C:/PROGRA~1/QGISPI~1/apps/qgis/./python/plugins',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\matplotlib-1.3.1-py2.7-win-amd64.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\nose-1.3.3-py2.7.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\tornado-4.0.1-py2.7-win-amd64.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\backports.ssl_match_hostname-3.4.0.2-py2.7.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\certifi-14.05.14-py2.7.egg',
> 'C:\\OSGeo4W\\apps\\Python27', 'C:\\OSGeo4W\\apps\\qgis\\bin',
> 'C:\\OSGeo4W\\apps\\qgis\\python', 'C:\\Python27\\Python2.7.9',
> 'C:\\Python27\\Python2.7.9\\Scripts',
> 'C:\\PROGRA~1\\QGISPI~1\\bin\\python27.zip',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\DLLs',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\plat-win',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\lib-tk',
> 'C:\\PROGRA~1\\QGISPI~1\\bin', 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\PIL',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\jinja2-2.7.2-py2.7.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\markupsafe-0.23-py2.7-win-amd64.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\pytz-2012j-py2.7.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\win32',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\win32\\lib',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\Pythonwin',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win-amd64.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\six-1.3.0-py2.7.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\xlrd-0.9.2-py2.7.egg',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\xlwt-0.7.5-py2.7.egg',
> u'C:/Users/lpalao/.qgis2//python',
> 'C:\\Users\\lpalao\\.qgis2\\python\\plugins\\mmqgis/forms',
> 'C:\\PROGRA~1\\QGISPI~1\\apps\\qgis\\python\\plugins\\fTools\\tools']
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] OpenLayers issue QGIS 2.10

[Leo Kris Palao]

Hi Enrique,

Thanks. It worked great.

-Leo

On Thu, Aug 27, 2015 at 4:49 PM, Enrique  wrote:

> It seems a bug of the plugin. I have found this workaround (QGIS 2.10.1
> Pisa, Win7) when you have Internet Access through proxy:
>
> 1. Settings>Options>Network>Use proxy for web access>Select httpproxy (or
> Sockets, i.e.)> Set any port value (i.e. 9)
> 2. Select again DefaultProxy>OK
>
> Go!
>
> This workaround is only required for the first load of OSM layers.
>
> Enrique.
>
>
> 2015-08-27 4:02 GMT+02:00 Leo Kris Palao :
>
>> Hi Guys,
>>
>> I want to use OpenLayers plugin in my QGIS 2.10 Pisa but getting an error
>> message. This error started when I upgraded my QGIS to version 2.8 Wien,
>> and until now, in my QGIS 2.10 still have this. Can I have some advise how
>> to solve this issue?
>>
>> I am using Windows 7 OS.
>>
>> Thanks in advance,
>> -Leo
>>
>> An error has occured while executing Python code:
>>
>> Traceback (most recent call last):
>>>   File
>>> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\weblayers\weblayer.py",
>>> line 97, in addLayer
>>> self._addLayerCallback(self)
>>>   File
>>> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\openlayers_plugin.py",
>>> line 130, in addLayer
>>> layer = OpenlayersLayer(self.iface, self._olLayerTypeRegistry)
>>>   File
>>> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py",
>>> line 330, in __init__
>>> self.olWebPage = OLWebPage(self)
>>>   File
>>> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py",
>>> line 48, in __init__
>>> proxy = getProxy()
>>>   File
>>> "C:/Users/lpalao/.qgis2/python/plugins\openlayers_plugin\tools_network.py",
>>> line 42, in getProxy
>>> proxy.setPort(settings.value("/proxyPort", 10, type=int))
>>> TypeError: unable to convert a QVariant of type 10 to a QMetaType of
>>> type 2
>>> Python version:
>>> 2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]
>>>
>>> QGIS version:
>>> 2.10.1-Pisa Pisa, d20c5b7
>>> Python path:
>>> ['C:/PROGRA~1/QGISPI~1/apps/qgis/./python/plugins\\processing',
>>> 'C:/PROGRA~1/QGISPI~1/apps/qgis/./python',
>>> u'C:/Users/lpalao/.qgis2/python', u'C:/Users/lpalao/.qgis2/python/plugins',
>>> 'C:/PROGRA~1/QGISPI~1/apps/qgis/./python/plugins',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\matplotlib-1.3.1-py2.7-win-amd64.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\nose-1.3.3-py2.7.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\tornado-4.0.1-py2.7-win-amd64.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\backports.ssl_match_hostname-3.4.0.2-py2.7.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\certifi-14.05.14-py2.7.egg',
>>> 'C:\\OSGeo4W\\apps\\Python27', 'C:\\OSGeo4W\\apps\\qgis\\bin',
>>> 'C:\\OSGeo4W\\apps\\qgis\\python', 'C:\\Python27\\Python2.7.9',
>>> 'C:\\Python27\\Python2.7.9\\Scripts',
>>> 'C:\\PROGRA~1\\QGISPI~1\\bin\\python27.zip',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\DLLs',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\plat-win',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\lib-tk',
>>> 'C:\\PROGRA~1\\QGISPI~1\\bin', 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\PIL',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\jinja2-2.7.2-py2.7.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\markupsafe-0.23-py2.7-win-amd64.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\pytz-2012j-py2.7.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\win32',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\win32\\lib',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\Pythonwin',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win-amd64.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\six-1.3.0-py2.7.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\xlrd-0.9.2-py2.7.egg',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\Python27\\lib\\site-packages\\xlwt-0.7.5-py2.7.egg',
>>> u'C:/Users/lpalao/.qgis2//python',
>>> 'C:\\Users\\lpalao\\.qgis2\\python\\plugins\\mmqgis/forms',
>>> 'C:\\PROGRA~1\\QGISPI~1\\apps\\qgis\\python\\plugins\\fTools\\tools']
>>
>>
>> ___
>> Qgis-user mailing list
>> Qgis-user@lists.osgeo.org
>> http://lists.osgeo.org/mailman/listinfo/qgis-user
>>
>
>
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user