Path: csiph.com!v102.xanadu-bbs.net!xanadu-bbs.net!feeder.erje.net!eu.feeder.erje.net!xlned.com!feeder5.xlned.com!newsfeed.xs4all.nl!newsfeed1.news.xs4all.nl!xs4all!post.news.xs4all.nl!not-for-mail
To: python-list@python.org
From: Mark Lawrence <breamoreboy@yahoo.co.uk>
Subject: Re: Drawing shaded area depending on distance with latitude and altitude coordinate
Date: Mon, 06 Jan 2014 20:16:53 +0000
References: <7f58342b-9b75-4a2c-b98e-cc43a4bb05a6@googlegroups.com>
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
User-Agent: Mozilla/5.0 (Windows NT 6.1; rv:24.0) Gecko/20100101 Thunderbird/24.2.0
In-Reply-To: <7f58342b-9b75-4a2c-b98e-cc43a4bb05a6@googlegroups.com>
Precedence: list
Newsgroups: comp.lang.python
Message-ID: <mailman.5052.1389039434.18130.python-list@python.org>
Lines: 69
NNTP-Posting-Host: 2001:888:2000:d::a6
Xref: csiph.com comp.lang.python:63332

On 06/01/2014 20:08, Isaac Won wrote:
> I have tried to make a plot of points with longitude and latitude coordinate, and draw shaded area with distance from one point. So, I thought that I could uae contourf function from matplotlibrary. My code is:
>      import haversine
>      import numpy as np
>      import matplotlib.pyplot as plt
>      with open(filin, 'r') as f:
>      arrays = [map(float, line.split()) for line in f]
>          newa = [[x[1],-x[2]] for x in arrays]
>
>      lat = np.zeros(275)
>      lon = np.zeros(275)
>      for c in range(0,275):
>          lat[c] = newa[c][0]
>          lon[c] = newa[c][1]
>
>      with open(filin, 'r') as f:
>          arrays = [map(float, line.split()) for line in f]
>      newa = [[x[1],-x[2]] for x in arrays]
>
>      lat = np.zeros(275)
>      lon = np.zeros(275)
>      for c in range(0,275):
>          lat[c] = newa[c][0]
>          lon[c] = newa[c][1]
>
>
>      dis = np.zeros(275)
>
>      for c in range(0,275):
>          dis[c] = haversine.distance(newa[0],[lat[c],lon[c]])
>
>      dis1 = [[]]*1
>
>      for c in range(0,275):
>          dis1[0].append(dis[c])
>
>
>      cs = plt.contourf(lon,lat,dis1)
>      cb = plt.colorbar(cs)
>
>      plt.plot(-lon[0],lat[0],'ro')
>      plt.plot(-lon[275],lat[275],'ko')
>      plt.plot(-lon[1:275],lat[1:275],'bo')
>      plt.xlabel('Longitude(West)')
>      plt.ylabel('Latitude(North)')
>      plt.gca().invert_xaxis()
>      plt.show()
>
> My idea in this code was that I could made a shaded contour by distance from a certain point which was noted as newa[0] in the code. I calculated distances between newa[0] and other points by haversine module which calculate distances with longitudes and latitudes of two points. However, whenever I ran this code, I got the error related to X, Y or Z in contourf such as:
>      TypeError: Length of x must be number of columns in z, and length of y must be number of rows.
>
> IF I use meshgrid for X and Y, I also get:
>      TypeError: Inputs x and y must be 1D or 2D.
>
> I just need to draw shaded contour with distance from one point on the top of the plot of each point.
>
> If you give any idea or hint, I will really apprecite. Thank you, Isaac
>

Sorry I can't help directly but can point you here 
https://lists.sourceforge.net/lists/listinfo/matplotlib-users or perhaps 
stackoverflow.

-- 
My fellow Pythonistas, ask not what our language can do for you, ask 
what you can do for our language.

Mark Lawrence