Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  349 1100 (S16E45 L=207)  1156    243        1 ( 282/1008- 417/1198)  66.8 E49-E41 S22-S08
      2  249 1133 (S16E55 L=197)   446     93        1 ( 221/1097- 272/1169)  59.0 E57-E53 S18-S13
      3 1861 1053 (N21W57 L=309)   458    144        0 (1822/1033-1883/1094)  51.3 W52-W59 N19-N24
      4  431  981 (S21E36 L=216)   282     63        1 ( 407/ 956- 458/ 997)  45.0 E37-E33 S22-S19
      5  180  818 (S37E58 L=194)  1374    314        1 ( 130/ 779- 228/ 882)  36.0 E64-E52 S39-S31
      6  401  783 (S34E36 L=216)   834    156        1 ( 336/ 743- 463/ 822)  -5.2 E41-E30 S35-S32
      7  799  764 (S25E07 L=246)  1739    474        1 ( 558/ 684- 998/ 815)   9.5 E23-W07 S29-S21
      8 1360  669 (S15W30 L=282)  1746    401        3 (1158/ 635-1575/ 721)  17.2 W15-W44 S18-S07
      9  554  545 (S46E20 L=232)   704     85        1 ( 516/ 506- 592/ 586) -15.9 E24-E16 S48-S42
     10  692  526 (S43E08 L=244)   256     85        1 ( 648/ 513- 724/ 546)  33.7 E11-E05 S45-S40
     11 1244  382 (S35W36 L=289)   285     40        0 (1214/ 365-1267/ 399)  56.3 W34-W37 S36-S33
     12 1404  280 (S36W54 L=307)   510     69        0 (1368/ 268-1438/ 300)  53.1 W51-W56 S37-S33
  

    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