I am planning to spend a year abroad for academic purposes in the upcoming year. Currently, I am in the process of identifying an appropriate laboratory or supervisor, but I am finding it challenging to locate the right connections. My research interests are primarily in the fields of 3D reconstruction and computer vision. I would be immensely grateful if anyone could offer some guidance or suggestions in this regard. Should you be willing to assist, please do not hesitate to reach out to me via my personal email: 374774222@qq.com. Thank you in advance for your consideration.
This is VastGaussian: Vast 3D Gaussians for Large Scene Reconstruction
unofficial implementation, since this is my first time to recreate the complete code from scratch, the code may have some errors, and the code writing may seem a bit naive compared to some experts. Lack of engineering skills. But I got my foot in the door. I couldn't find any implementation of VastGaussian on the web, so I gave it a try.
-
Refactor repository Important ! -
Camera-position-based region division is implemented -
Position-based data selection is implemented -
Visibility-based camera selection is implemented -
Coverage-based point selection is implemented -
Decoupled Appearance Modeling is implemented -
Fixed high RAM and VRAM usage issues, it is now possible to train the first few iteration on a 12GB GPU and 32GB of RAM -
Fix Decoupled Appearance Modeling that didn't work well -
Seamless Merging is implemented -
For non-standard scenes by manual Manhattan alignment -
Parallel training ofThanks livion a lot.$m\times n$ regions on$n$ GPUs is implemented after dividing the point cloud -
Refine Data Partition -
Refine Seamless Merging -
update color correct and LPIPS (using alex net) while evaluation
- In order not to change the original directory structure of 3DGS, I added a new
scene/vastgs
module to store VastGaussian. Part of the code I called the existing functions in thescene
folder. Also to fix theimport
error, I moved the Scene class into the datasets.py folder
- The naming of the files is consistent with the method mentioned in the paper for easy reading
data_partition.py
corresponding to theProgressive Data Partitioning
in the paper.
scene/vastgs/appearance_network.py
corresponding to theDecoupled Appearance Modeling
in the paper. I refer to the implementation of gaussian-opacity-fields
graham_scan.py
convex hull calculation is used to project the partition cube onto the camera plane and calculate the intersection of the projected region and the image region when implementing Visibility based camera selection.
seamless_merging.py
corresponding to theSeamless Merging
in the paper.
- I have added a new file
train_vast.py
to modify the process of training VastGaussian, if you want to train the original 3DGS, please usetrain.py
. - The paper mentioned
Manhattan world alignment, so that the Y-axis of the world coordinate is perpendicular to the ground plane
, I asked the experts to know that this thing can be adjusted manually usingthreejs
: https://threejs.org/editor/ or the softwarecloudcompare
, after manually adjusting the scene you get the --position and --rotation parameters, just take them as command line arguments and train.
- Now you can adjust your point cloud so that the ground is perpendicular to the y-axis and the boundaries are as parallel as possible to the x- and z-axis with the help of the options on the left, but of course you can also use the editing area on the right to directly enter the corresponding values. - Then you can get the appropriate parameters in the right edit area.
- After importing your initial point cloud via File-->Import you can check if your initial point cloud needs to be Manhattan aligned, if it looks like this
- Open the cloudcompare software and import the sparse point cloud into the software.
- Use the
cross section
tool in the toolbar to reduce the scope of the point cloud to only the areas you are interested in (for easier alignment), or you can leave them out. Then you can use the toolbar on the left of the red arrow to adjust your viewing Angle (there are 6 viewing angles), and finally drag the arrow pointed by the green arrow to adjust the area you are interested in.
- After you have adjusted the point cloud, you can export it as a new point cloud, noting that there is no transformation of coordinates involved. Then close the box on the right. Select the exported point cloud and deselect the initial point cloud.
- Use the
Translate/Rotate
tool on the toolbar to adjust the pose of the point cloud. ClickRotation
to select the axis around which you want to rotate. If you want to adjust both the rotation matrix and the transfer vector, you can tickTy Ty Tz
Also select the toolbar on the left to adjust the viewing Angle.
- The Manhatton alignment mentioned in the paper can be realized by manually adjusting the pose of the point cloud so that the boundary frame x and z axis of the point cloud are parallel. And you can get the transformation matrix relative to the initial point cloud after this adjustment in the software console. Let's call it
A1
- If you only adjust once, then A1 is the final transformation matrix (
A=A1
), if you adjust the pose of the point cloud several times in A row, assuming 3 adjustments, and get the transformation matrixA1 A2 A3
, then the final transformation matrix isA= A3*A2*A1
- Enter the resulting transformation matrix into the command line.
- In the process of implementation, I used a small range of data provided by 3DGS for testing. Larger data can not run on the native computer, and a large range of data requires at least 32G video memory according to the instructions of the paper.
- In the implementation process, some operations in the paper, the author is not very clear about the details, so some implementation is based on my guess and understanding to complete, so my implementation may have some bugs, and some implementation may be a little stupid in the eyes of the expert, if you find problems in the use of the process, please contact me in time, progress together.
You can refer to the parameter Settings in the files train_vast.sh
and eval_vast.sh
- The data format is the same as 3DGS, and the training command is basically the same as 3DGS. I didn't make too many personalized changes, you can refer to the following command (see
arguments/parameters.py
for more parameters): if you want to perform manhattan alignment:
Using threejs
for Manhattan alignment
python train_vast.py -s datasets/xxx --exp_name xxx --manhattan --platform tj --pos "xx xx xx" --rot "xx xx xx"
Using cloudcompare
for Manhattan alignment
# The 9 elements of the rotation matrix should be filled in rot
python train_vast.py -s datasets/xxx --exp_name xxx --manhattan --platform cc --pos "xx xx xx" --rot "xx xx xx xx xx xx xx xx xx"
Train without Manhattan alignment:
This may cause an error because the Manhattan alignment is not performed, which may cause an error when the data partition is executed
python train_vast.py -s datasets/xxx --exp_name test
I get the preprocessed data from https://vastgaussian.github.io/, and implement Manhattan alignment, you can use my pos and rot params.
# train rubble
python train_vast.py -s ../datasets/Mill19/rubble \
--exp_name rubble \
--manhattan \
--eval \
--llffhold 83 \
--resolution 4 \
--pos "25.607364654541 0.000000000000 -12.012700080872" \
--rot "0.923032462597 0.000000000000 0.384722054005 0.000000000000 1.000000000000 0.000000000000 -0.384722054005 0.000000000000 0.923032462597" \
--m_region 3 \
--n_region 3 \
--iterations 60_000
# train building
python train_vast.py -s ../datasets/Mill19/building \
--exp_name building \
--manhattan \
--eval \
--llffhold 83 \
--resolution 4 \
--pos "-62.527942657471 0.000000000000 -15.786898612976" \
--rot "0.932374119759 0.000000000000 0.361494839191 0.000000000000 1.000000000000 0.000000000000 -0.361494839191 0.000000000000 0.932374119759" \
--m_region 3 \
--n_region 3 \
--iterations 60_000
I added new parameters in arguments/__init__.py
New Parameters for train_vast.py
Experiment name
store_true
, Whether to perform Manhattan alignment
Platform for Manhattan alignment, choose in "cc" and "tj" which means "cloudcompare" and "threejs"
Translation vector
rotate matrix
default=None, transformational matrix
the number of regions in the x direction
the number of regions in the z direction
The rate of boundary expansion
Airspace-aware visibility rate
Evaluation method: During training, every 83 images were used as a test image (llffhold=83
). For rubble, 21 images were selected and added to the test set.
# render rubble
python render.py -s ../datasets/Mill19/rubble \
--exp_name rubble \
--eval \
--manhattan \
--resolution 4 \
--pos "25.607364654541 0.000000000000 -12.012700080872" \
--rot "0.923032462597 0.000000000000 0.384722054005 0.000000000000 1.000000000000 0.000000000000 -0.384722054005 0.000000000000 0.923032462597" \
--load_iteration 60_000
# eval rubble
python metrics.py -m output/rubble
# render building
python render.py -s ../datasets/Mill19/building \
--exp_name building \
--manhattan \
--eval \
--llffhold 83 \
--resolution 4 \
--pos "-62.527942657471 0.000000000000 -15.786898612976" \
--rot "0.932374119759 0.000000000000 0.361494839191 0.000000000000 1.000000000000 0.000000000000 -0.361494839191 0.000000000000 0.932374119759" \
--load_iteration 60_000
# eval building
python metrics.py -m output/building
-
Urbanscene3D
: https://github.com/Linxius/UrbanScene3D -
Mill-19
: https://opendatalab.com/OpenDataLab/Mill_19/tree/main/raw
https://vastgaussian.github.io/ have uploaded the pre-processed data for
Urbanscene3D
andMill-19
- test data for this implementation: https://repo-sam.inria.fr/fungraph/3d-gaussian-splatting/datasets/input/tandt_db.zip
Happily, we now have several contributors working on the project, and we welcome more contributors to join us to improve the project. Thank you all for your work.