def extract(self, point_x, point_y, radius=<NUM_LIT:0>):
row, col = map_pixel(point_x, point_y, self.x_cell_size, self.y_cell_size,<EOL>self.xmin, self.ymax)<EOL>col2 = np.abs(radius/self.x_cell_size).astype(int)<EOL>row2 = np.abs(radius/self.y_cell_size).astype(int)<EOL>return GeoRaster(self.raster[max(row-row2, <NUM_LIT:0>):min(row+row2+<NUM_LIT:1>, self.shape[<NUM_LIT:0>]),max(col-col2, <NUM_LIT:0>):min(col+col2+<NUM_LIT:1>, self.shape[<NUM_LIT:1>])], self.geot,<EOL>nodata_value=self.nodata_value,projection=self.projection, datatype=self.datatype)<EOL>
geo.extract(x, y, radius=r)

Return subraster of raster geo around location (x,y) with radius r
where (x,y) and r are in the same coordinate system as geo