Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1443 1599 (N28W41 L=272)   463     97        2 (1400/1579-1499/1619)  -7.1 W37-W45 N26-N30
      2  968 1593 (N37W07 L=238)   289     26        0 ( 944/1570- 987/1612) -76.0 W05-W08 N34-N38
      3  464 1316 (N27E33 L=198)  1689    405        2 ( 335/1216- 668/1396)  10.5 E43-E16 N22-N29
      4  218 1273 (N29E58 L=173)  1586    173        1 ( 142/1231- 289/1309) -18.4 E72-E48 N24-N32
      5 1501  792 (S21W26 L=257)  1227    442        0 (1377/ 648-1627/ 953)  35.0 W16-W36 S28-S12
      6  274  813 (N00E54 L=177)   516    126        0 ( 256/ 743- 296/ 866)  50.2 E57-E51 S03-N03
      7  347  607 (S13E55 L=176)   971    144        2 ( 313/ 532- 375/ 669)  71.6 E58-E52 S17-S08
      8 1122  653 (S23E00 L=231)   255     82        0 (1073/ 635-1157/ 669) -14.0 E03-W01 S24-S21
  

    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