python - Getting Rpy2 to work with rgdal to spatially subset points -

-

so have r code works. takes bunch of points data , spatially subsets against shapefile in method of http://robinlovelace.net/r/2014/07/29/clipping-with-r.html.

#data .csv file lon lat points data_points <- spatialpoints(data) proj4string(data_points) <- crs("+proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0,0") data_ll <- sptransform(data_points, crs("+proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0,0")) melbourne <- readogr("melbourne_australia.land.coastline","melbourne_australia_land_coast") #this shapefile https://mapzen.com/data/metro-extracts subset <- data_ll[melbourne,]

plot(melbourne) points(subset)

i'm trying convert corresponding rpy2 script. far have;

import pandas pd import numpy np import rpy2.robjects ro import rpy2.robjects.numpy2ri rpy2.robjects.packages import importr rgdal = importr('rgdal') base = importr('base')  rpy2.robjects.numpy2ri.activate() data = pd.read_csv('sim.csv') data = data.values coordinates = ro.r['coordinates'] proj4string = ro.r['proj4string'] sptransform = ro.r['sptransform'] readogr = ro.r['readogr'] spatialpoints = ro.r['spatialpoints'] crs = ro.r['crs'] class_r = ro.r['class']  key = crs("+proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0,0") data_points = spatialpoints(data, proj4string = key) data_ll = sptransform(data_points, key) melbourne = readogr("melbourne_australia.land.coastline", "melbourne_australia_land_coast")  subset = data_ll[melbourne,]

which fails on last line error typeerror: 'rs4' object not subscriptable. have idea of going on?

one way convert r code function , import package.

this r code:

library(rgdal) library(sp)  #import data data <- read.csv("sim.csv", header = f)  subset_points <- function(data){     data_points <- spatialpoints(data, proj4string=crs("+proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0,0"))     data_ll <- sptransform(data_points, crs("+proj=longlat +datum=wgs84 +no_defs +ellps=wgs84 +towgs84=0,0,0"))     sink("/dev/null")       melbourne <- readogr("melbourne_australia.land.coastline", "melbourne_australia_land_coast")     sink()     subset <- data_ll[melbourne,]      final <- as.data.frame(subset)     return(final) }

and python code:

import rpy2.robjects ro import rpy2.robjects.numpy2ri rpy2.robjects.packages import importr  rpy2.robjects.packages import signaturetranslatedanonymouspackage open('subset_data.r') fh:     rcode = os.linesep.join(fh.readlines())     subset = signaturetranslatedanonymouspackage(rcode, "subset") rpy2.robjects.numpy2ri.activate()  data = pd.read_csv('sim.csv') data = data.values final = subset.subset_points(data) print(np.array(final).t)

Comments

Post a Comment

Popular posts from this blog

php - Invalid Cofiguration - yii\base\InvalidConfigException - Yii2 -

-
Image
i've installed yii2-user module/dektrium yii2-app-basic application command composer require "dektrium/yii2-user:0.9.*@dev" config/console.php return [ . . 'modules' => [ 'gii' => 'yii\gii\module', 'user' => [ 'class' => 'dektrium\user\module', ], ], . . ] config/web.php 'components' => [ . . /* 'user' => [ 'identityclass' => 'app\models\user', 'enableautologin' => true, ], */ 'modules' => [ 'user' => [ 'class' => 'dektrium\user\module', ], ], . . ] after that, run command $ php yii migrate/up --migrationpath=@vendor/dektrium/yii2-user/migrations updating database schema . but, when run http://localhost/mylawsuit/yii/web/index.php?r=u
Read more

How to show in django cms breadcrumbs full path? -

-
need advice - how show in breadcrumbs whole path article or page? i've got code: {% ance in ancestors %} <li> {% if not forloop.last %} <a href="{{ ance.get_absolute_url }}">{{ ance.get_menu_title }}</a> {% else %} <span class="active">{{ ance.get_menu_title }}</span> {% endif %} </li> {% endfor %} but shows "home" , current page. how modify show whole path in breadcrumbs? i've found in file menu_tags.py code calls ancestors. @ least there class showbreadcrumbs: class showbreadcrumb(inclusiontag): """ shows breadcrumb node has same url current request - start level: after level should breadcrumb start? 0=home - template: template used render breadcrumb """ name = 'show_breadcrumb' template = 'menu/dummy.html' options = options( argument('start_level', default=0, require
Read more

ruby on rails - npm error: tunneling socket could not be established, cause=connect ETIMEDOUT -

-
i have app uses rails framework , implements angularjs part of front end. i have pushed heroku , have heroku toolbelt installed, when try migrate db using "heroku run rake db:migrate" receive following error(s): installing core plugins heroku-cli-addons, heroku-apps, heroku-fork, heroku-git, heroku-local, heroku-run, heroku-status... error installing package. try running again gode_debug=info see more output. ! `run` not heroku command. ! perhaps meant `-h`, `2fa`, `auth`, `join`, `open`, `orgs`, `pg`, `ps` or `rake`. ! see `heroku help` list of available commands. i run command "gode_debug=info heroku run rake db:migrate" , receive error: npm err! darwin 14.5.0 npm err! argv "/users/christopher_pelnar/.heroku/node-v4.2.1-darwin-x64/bin/node" "/users/christopher_pelnar/.heroku/node-v4.2.1-darwin-x64/lib/node_modules/npm/cli.js" "install" "heroku-cli-addons" "heroku-apps" "heroku-fork
Read more