Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  747 1491 (N41E09 L=352)  1341    290        0 ( 651/1437- 876/1573) -41.2 E17-W00 N35-N49
      2 1158 1349 (N23W15 L= 16)   263    129        1 (1142/1291-1181/1399)  81.9 W14-W17 N19-N26
      3  320 1213 (N30E48 L=312)  2736    421        9 ( 161/1084- 449/1337) -22.0 E78-E36 N20-N39
      4 1738  939 (S14W45 L= 46)   547     88        0 (1712/ 887-1763/ 985)  33.7 W42-W48 S17-S11
      5  544  950 (N11E31 L=330)   316     80        0 ( 509/ 919- 585/ 969)   0.0 E34-E28 N09-N11
      6  237  780 (N04E62 L=299)   700     92        2 ( 207/ 760- 269/ 807)  12.5 E66-E57 N03-N06
      7  828  718 (S07E18 L=343)   265    100        1 ( 809/ 670- 860/ 763)  45.0 E19-E14 S10-S05
      8  741  639 (S11E25 L=336)   625    256        0 ( 641/ 584- 819/ 689)   4.4 E32-E19 S13-S07
      9  460  644 (S06E45 L=316)   298    134        0 ( 420/ 621- 501/ 670)  -8.1 E48-E41 S07-S04
     10 1691  527 (S42W42 L= 43)   744     77        3 (1660/ 492-1718/ 559)   9.5 W38-W46 S44-S39
     11 1248  426 (S36W02 L=  3)   488     93        0 (1206/ 406-1292/ 448)   0.0 E01-W05 S37-S34
     12  922  346 (S34E24 L=337)   607    132        2 ( 896/ 290- 945/ 412)  48.4 E27-E19 S38-S29
  

    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