Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1309 1516 (N24W32 L=213)  3445    728        8 ( 997/1462-1587/1551) -24.1 W12-W51 N15-N33
      2  808 1347 (N28E04 L=177)   372     98        2 ( 780/1306- 844/1391)  -9.5 E05-E02 N25-N31
      3 1282 1214 (N08W19 L=201)   237     75        2 (1266/1175-1304/1257)  78.7 W18-W20 N06-N11
      4 1918 1118 (S16W66 L=247)  1916    263        2 (1873/1004-1953/1216)  85.9 W62-W70 S23-S09
      5 1818 1149 (S11W55 L=236)  1116    198        1 (1765/1092-1864/1198) -58.0 W51-W60 S14-S06
      6 1879 1005 (S21W58 L=239)   418     77        1 (1866/ 970-1894/1038)  36.9 W56-W60 S22-S18
      7 1619  751 (S27W30 L=211)  2574    635        4 (1445/ 546-1719/1030)  38.2 W13-W42 S37-S13
      8 1200  909 (S06W07 L=188)   233     76        0 (1150/ 891-1245/ 928)  -9.5 W03-W09 S07-S05
      9  508  834 (N06E35 L=146)   313     65        2 ( 487/ 805- 537/ 863)  36.9 E37-E32 N04-N08
     10  625  800 (N02E29 L=153)   381    156        0 ( 554/ 745- 688/ 853) -58.0 E31-E26 S02-N06
     11 1449  707 (S25W17 L=198)   596    176        0 (1386/ 626-1516/ 800)  26.6 W11-W23 S28-S20
     12 1130  513 (S28E08 L=173)   387    128        0 (1091/ 458-1155/ 570)  40.6 E12-E04 S30-S24
     13  750  433 (S22E34 L=147)   508    182        3 ( 716/ 371- 809/ 477)  21.8 E38-E28 S25-S20
     14  503  408 (S18E53 L=128)   515     90        1 ( 481/ 371- 526/ 452)  80.5 E54-E51 S21-S15
     15  711  304 (S30E44 L=137)   393     78        1 ( 675/ 282- 747/ 321)  -9.5 E46-E40 S31-S29
     16  925  276 (S38E32 L=149)   359    100        2 ( 886/ 260- 981/ 293)  -7.1 E35-E27 S39-S36
     17  772  245 (S36E44 L=137)   466     74        1 ( 741/ 226- 810/ 269)  -9.5 E47-E40 S37-S34
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y