regridder to_dataset method returns undersized n

[Huite]

I'm having a target grid with ny * nx values, but the resulting __regrid_n is smaller.

This is likely due to the fact that the scipy matrix is initialized without an explicit shape argument. In that case, it'll infer the shape using the highest index.