_{\Omega_j} \bm{\alpha} R_j- \bm{1}_jb_j \|_F^2,$$ Once $\{R_j\}$ and $\{b_j\}$ are obtained, we can find all the adjusted coordinates $P_{\Omega_j}= Q_j R_j + \bm{1}_j b_j$, which leads to the global coordinate $P$ ultimately.

Different from the synchronization method and the spectral relaxation method considered in [@cucuringu2012sensor; @chaudhury2015global], we design the following method to solve the nonconvex problem by updating $\bm{\alpha}$, $\{R_j\}$ and $\{b_j\}$ iteratively: $$\label{gradient descent 1}
\begin{cases}
\bm{\alpha}^{k+1}= \bm{\alpha}^k - (\sum \limits_{j=1}^{L} \Phi_{\Omega_j}^{\top}\Phi_{\Omega_j})^{-1} \nabla _{\bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k),\\
R_j^{k+1}= \mathcal{K}(R_j^{k}, \nabla _{R_j}\E(\bm{\alpha}^{k+1}, R_j^k, b_j^k)), \qquad \forall 1 \leq j \leq L,\\
\displaystyle b_j^{k+1} = b_j^{k} -  \frac{1}{\text{Card}(\Omega_j)} \nabla _{b_j}\E(\bm{\alpha}^{k+1}, R_j^{k+1}, b_j^k), ~\quad \forall 1 \leq j \leq L.\\
\end{cases}$$ where we update $\bm{\alpha}$ and $\{b_j\}$ using Newton’s method and it is straightforward to check $$\label{eqn:gradientofstitching}
\begin{cases}
\nabla _{\bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k) = - \sum \limits_{j=1}^{L}  \Phi_{\Omega_j}^{\top}(Q_{j} -  \Phi_{\Omega_j} \bm{\alpha}^k R_j^{k}- \bm{1}_jb_j^{k} ) (R_
_ { \Omega_j } \bm{\alpha } R_j- \bm{1}_jb_j \|_F^2,$$ Once $ \{R_j\}$ and $ \{b_j\}$ are obtained, we can find all the adjusted coordinates $ P_{\Omega_j}= Q_j R_j + \bm{1}_j b_j$, which lead to the ball-shaped coordinate $ P$ ultimately. 

 Different from the synchronization method acting and the spectral rest method consider in [ @cucuringu2012sensor; @chaudhury2015global ], we design the following method acting to resolve the nonconvex problem by update $ \bm{\alpha}$, $ \{R_j\}$ and $ \{b_j\}$ iteratively: $ $ \label{gradient descent 1 } 
 \begin{cases } 
 \bm{\alpha}^{k+1}= \bm{\alpha}^k - (\sum \limits_{j=1}^{L } \Phi_{\Omega_j}^{\top}\Phi_{\Omega_j})^{-1 } \nabla _ { \bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k),\\ 
 R_j^{k+1}= \mathcal{K}(R_j^{k }, \nabla _ { R_j}\E(\bm{\alpha}^{k+1 }, R_j^k, b_j^k) ), \qquad \forall 1 \leq j \leq L,\\ 
 \displaystyle b_j^{k+1 } = b_j^{k } -   \frac{1}{\text{Card}(\Omega_j) } \nabla _ { b_j}\E(\bm{\alpha}^{k+1 }, R_j^{k+1 }, b_j^k), ~\quad \forall 1 \leq j \leq L.\\ 
 \end{cases}$$ where we update $ \bm{\alpha}$ and $ \{b_j\}$ use Newton ’s method and it is square to check $ $ \label{eqn: gradientofstitching } 
 \begin{cases } 
 \nabla _ { \bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k) = - \sum \limits_{j=1}^{L }   \Phi_{\Omega_j}^{\top}(Q_{j } -   \Phi_{\Omega_j } \bm{\alpha}^k R_j^{k}- \bm{1}_jb_j^{k }) (R _
_{\Omeha_j} \bm{\alpha} R_j- \bm{1}_jb_j \|_F^2,$$ Okce $\{R_j\}$ and $\{b_j\}$ arg ibtainxd, we czn find xll the adjusted coordinates $P_{\Imega_h}= Q_j R_j + \bm{1}_j b_j$, which ueads to nhe globao coiedinate $P$ nmtimately.

Diffsvent yrim the synchrokization medhod and the s[eztxal relaxation method considered in [@sucurinbu2012densor; @chaudhuty2015glonwl], ws design the following method to sklve tht nonconvex problek by updating $\bm{\alpha}$, $\{R_j\}$ ajd $\{b_u\}$ iteratively: $$\labep{gradient dgacegr 1}
\begin{cases}
\cm{\alpha}^{k+1}= \bm{\alpha}^k - (\sum \limits_{j=1}^{L} \Phi_{\Omega_j}^{\top}\Phi_{\Omega_j})^{-1} \nablc _{\bm{\alpha}}\E(\bn{\aophw}^n, R_j^k, b_j^k),\\
R_j^{j+1}= \matrcal{K}(R_j^{k}, \nabla _{R_j}\E(\bm{\dlpha}^{k+1}, T_j^k, b_j^k)), \qquad \noraln 1 \oeq j \leq L,\\
\displaystyne b_j^{k+1} = b_j^{k} -  \frac{1}{\jext{Card}(\Omagc_j)} \nabla _{b_j}\E(\bm{\alpha}^{k+1}, E_j^{j+1}, b_j^k), ~\quag \fofqll 1 \ltq o \lsq L.\\
\enf{caaes}$$ where se update $\bn{\alpha}$ and $\{b_j\}$ using Nqqton’s method znd it if straightforward to check $$\label{eqn:gradpentkfstitching}
\begin{cases}
\navla _{\bm{\alpha}}\E(\bm{\alpha}^k, T_j^k, b_j^k) = - \fum \limits_{j=1}^{L}  \Phi_{\Omega_j}^{\top}(Q_{j} -  \Phi_{\Omega_j} \bm{\alpha}^k R_j^{k}- \um{1}_gb_j^{j} ) (F_
_{\Omega_j} \bm{\alpha} R_j- \bm{1}_jb_j \|_F^2,$$ Once $\{R_j\}$ are we can all the adjusted \bm{1}_j which leads to global coordinate $P$ Different from the synchronization method and spectral relaxation method considered in [@cucuringu2012sensor; @chaudhury2015global], we design the following method to the nonconvex problem by updating $\bm{\alpha}$, $\{R_j\}$ and $\{b_j\}$ iteratively: $$\label{gradient descent 1} \bm{\alpha}^{k+1}= - \limits_{j=1}^{L} \nabla _{\bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k),\\ R_j^{k+1}= \mathcal{K}(R_j^{k}, \nabla _{R_j}\E(\bm{\alpha}^{k+1}, R_j^k, b_j^k)), \qquad \forall 1 \leq j \leq \displaystyle b_j^{k+1} = b_j^{k} - \frac{1}{\text{Card}(\Omega_j)} \nabla _{b_j}\E(\bm{\alpha}^{k+1}, b_j^k), ~\quad \forall 1 j \leq L.\\ \end{cases}$$ where update and $\{b_j\}$ Newton’s and is straightforward to $$\label{eqn:gradientofstitching} \begin{cases} \nabla _{\bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k) = - \sum \limits_{j=1}^{L} \Phi_{\Omega_j}^{\top}(Q_{j} - \Phi_{\Omega_j} \bm{\alpha}^k R_j^{k}- \bm{1}_jb_j^{k} )
_{\Omega_j} \bm{\alpha} R_j- \bm{1}_jb_j \|_F^2,$$ Once $\{r_j\}$ and $\{b_j\}$ are ObtaiNed, We cAn Find All tHe adjusted coorDInatEs $P_{\Omega_j}= Q_j R_j + \bm{1}_j b_j$, whicH leadS tO The gLObAl cooRdinate $p$ UlTIMatElY.

DIffErENt From tHe sYnchronIzation metHod AnD the spectral RElAxation metHod Considered in [@CucUringu2012SeNsoR; @ChaudHurY2015globAl], we deSIgn the Following MeTHod to sOLve the nONCoNvex Problem by updating $\BM{\aLPha}$, $\{R_j\}$ and $\{b_j\}$ iterAtivelY: $$\lABeL{GRadIenT descent 1}
\beGiN{caseS}
\Bm{\alpha}^{K+1}= \Bm{\ALPHa}^k - (\SUm \limits_{j=1}^{L} \Phi_{\omega_j}^{\top}\PhI_{\omeGa_j})^{-1} \nabLa _{\Bm{\aLPha}}\E(\bm{\Alpha}^K, R_J^K, b_j^K),\\
R_j^{k+1}= \mathcal{k}(R_j^{k}, \Nabla _{R_j}\E(\bM{\alpha}^{K+1}, r_j^k, b_j^k)), \qQUad \foraLl 1 \leq j \Leq l,\\
\diSplaYStYlE b_j^{K+1} = b_J^{K} -  \frAC{1}{\tExt{cArd}(\omega_j)} \naBlA _{b_J}\E(\bm{\aLpha}^{K+1}, r_J^{K+1}, B_j^k), ~\qUad \ForaLl 1 \leq J \leq L.\\
\end{cases}$$ WheRe we UPdaTe $\bm{\aLpha}$ aNd $\{b_j\}$ UsIng NeWton’s mEthod AnD it is straightfoRwarD to check $$\lAbeL{eQn:gRaDientOFstitcHinG}
\beGin{caseS}
\nabla _{\bM{\AlpHa}}\e(\BM{\AlPha}^k, R_j^k, b_j^k) = - \sum \limitS_{j=1}^{l}  \pHi_{\omega_j}^{\toP}(Q_{j} -  \Phi_{\oMeGa_J} \Bm{\alpha}^k r_j^{K}- \bm{1}_Jb_j^{k} ) (r_
_{\Omega_j} \bm{\alpha} R_ j- \bm{1}_ jb_j\|_ F^2 ,$ $ On ce $ \{R_j\}$ and $ \ {b_j \}$ are obtained, we c an fi nd allt he  adju sted co o rd i n ate s$P _{\ Om e ga _j}=Q_j  R_j +\bm{1}_j b _j$ ,which leadst othe global  co ordinate $P$  ul timate ly .
D iffer ent  from  the s y nchron ization m et h od and the spe c t ra l re laxation method c o ns i dered in [@cuc uringu 20 1 2s e n sor ; @ chaudhury2 01 5glob a l], wed es i g n  th e  following me thod to sol v e t he non co nve x  probl em by  u p dat ing $\bm{\a lpha }$, $\{R_ j\}$ a n d $\{b_ j \}$ ite rative ly:  $$ \lab e l{ gr adi en t  de s ce nt1 }
\ begin{ca se s} 
\bm{ \alp h a } ^ {k+1 }=\bm{ \alph a}^k - (\sum\li mits _ {j= 1}^{L } \Ph i_{\ Om ega_j }^{\to p}\Ph i_ {\Omega_j})^{-1 } \n abla _{\b m{\ al pha }} \E(\b m {\alph a}^ k,R_j^k,b_j^k), \ \
R _j ^ { k +1 }= \mathcal{K}(R_j ^{ k } ,\nabla _ {R_j}\ E (\ bm { \alpha}^ {k +1} , R_ j ^ k, b_ j^k) ) ,\qquad \ forall 1\l eq j \l eq  L,\\\d isp lay style b_j^ {k+1}= b_j^{k } - \ frac{1}{\text{ C ard}(\Omega_j ) }\ n ab l a _{ b_j }\E(\bm{\al pha} ^ {k+1 }, R _ j^ {k+ 1 }, b_ j^k),  ~ \ qu a d \forall 1 \leq j\l eq L.\ \
\en d{cases}$$ wh ere we upd a t e  $\bm{\a lpha } $a nd $\{b_j\}$ u singNewton’s m e thod and  it i s straig htforward t o check$$\ lab el{ eqn : g ra dientofstitch i n g}
\ be gin{cas es} 
\nabla  _{ \bm {\a lph a} }\E(\bm{\ alpha}^k ,R_ j^ k,  b_ j^k)=  - \sum\l imi ts _{j =1}^{ L }  \Ph i_{\O mega _j }^ { \to p}(Q_{j }  - \Phi _{ \O mega _j}  \ bm{\a lpha } ^kR_j^{k} - \bm{1}_ jb_ j ^{k}  )  ( R_
_{\Omega_j} \bm{\alpha}_R_j- \bm{1}_jb_j_\|_F^2,$$ Once $\{R_j\}$ and_$\{b_j\}$ are_obtained,_we can_find_all the adjusted_coordinates $P_{\Omega_j}= Q_j_R_j + \bm{1}_j b_j$,_which leads to_the_global coordinate $P$ ultimately.

Different from the synchronization method and the spectral relaxation method considered_in_[@cucuringu2012sensor; @chaudhury2015global],_we_design_the following method to solve_the nonconvex problem by updating_$\bm{\alpha}$, $\{R_j\}$_and $\{b_j\}$ iteratively: $$\label{gradient descent 1}
\begin{cases}
\bm{\alpha}^{k+1}= \bm{\alpha}^k -_(\sum_\limits_{j=1}^{L} \Phi_{\Omega_j}^{\top}\Phi_{\Omega_j})^{-1} \nabla__{\bm{\alpha}}\E(\bm{\alpha}^k, R_j^k, b_j^k),\\
R_j^{k+1}= \mathcal{K}(R_j^{k}, \nabla _{R_j}\E(\bm{\alpha}^{k+1}, R_j^k, b_j^k)), \qquad_\forall 1 \leq j \leq L,\\
\displaystyle_b_j^{k+1} = b_j^{k}_-__\frac{1}{\text{Card}(\Omega_j)} \nabla _{b_j}\E(\bm{\alpha}^{k+1}, R_j^{k+1},_b_j^k), ~\quad \forall 1 \leq j_\leq L.\\
\end{cases}$$ where we update $\bm{\alpha}$_and $\{b_j\}$ using Newton’s method and it_is straightforward to check $$\label{eqn:gradientofstitching}
\begin{cases}
\nabla _{\bm{\alpha}}\E(\bm{\alpha}^k,_R_j^k, b_j^k) = - \sum_\limits_{j=1}^{L} _\Phi_{\Omega_j}^{\top}(Q_{j} -  \Phi_{\Omega_j} \bm{\alpha}^k_R_j^{k}- \bm{1}_jb_j^{k} )_(R_