The function "houghpoint()" will apply the Hough transform to a collection of points, interpolating the offset for each angle bin. This function also allows each point to have an associated weight. The "houghimg()" wrapper can be used to apply the transform to images. If you want a different origin for the transform, consider modifying "houghimg()" to use different coordinates.
Example usage: >> HT = houghpoint(XYcoords,anglebins,offsetbins);
Use ">> help houghpoint" and ">> help houghimg" for more usage details.
Additionally, these function can produce a matrix that can be used to apply the Hough transform to an image using a matrix-vector multiply instead of a function. This is slower for a single image (especially images with many zeros), but for many inputs (of the same size for images, or same point locations for point clouds) it is faster to compute the matrix form once and apply it to the images of interest.