Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  350 1446 (N14E56 L=282)  2110    444        9 ( 248/1285- 454/1615)  73.7 E62-E51 N02-N26
      2  418 1426 (N14E49 L=289)   311     69        2 ( 388/1396- 449/1447)  53.1 E51-E46 N11-N15
      3 1667 1196 (N21W41 L= 18)   293    130        0 (1638/1118-1685/1248) -60.3 W37-W42 N16-N23
      4  173 1097 (S10E65 L=273)  2314    399        6 ( 128/ 981- 204/1218) -35.8 E78-E59 S16-S03
      5  776 1147 (N03E17 L=321)   235     41        2 ( 747/1125- 801/1173)  18.4 E18-E15 N01-N04
      6 1334 1092 (N09W18 L=355)   593    145        2 (1284/1039-1394/1161) -35.5 W13-W21 N05-N11
      7  352  942 (S16E44 L=294)   366    154        3 ( 296/ 922- 404/ 971)  24.0 E48-E39 S18-S13
      8 1618  557 (S19W51 L= 29)   485    161        1 (1592/ 509-1651/ 599) -24.0 W46-W55 S21-S15
  

    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