Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1726 1359 (N20W55 L= 32)   781     93        3 (1682/1315-1775/1409) -20.6 W51-W60 N16-N22
      2  225 1245 (N20E62 L=276)   406     72        0 ( 206/1224- 248/1275) -53.1 E64-E58 N18-N22
      3 1621 1232 (N13W42 L= 19)   317    100        0 (1564/1219-1665/1245) -14.0 W37-W45 N11-N14
      4  599  986 (N04E27 L=310)   801    197        3 ( 565/ 890- 629/1080)  73.3 E29-E24 S01-N09
      5  289  895 (S01E52 L=285)  2940    491        3 ( 234/ 730- 434/1071) -53.1 E57-E41 S11-N09
      6 1269  735 (S15W14 L=351)   424    111        1 (1223/ 700-1309/ 762) -45.0 W10-W16 S17-S12
      7  488  686 (S14E38 L=299)   331     51        1 ( 452/ 665- 516/ 704)   0.0 E40-E35 S15-S13
      8  482  577 (S22E41 L=296)   807    187        0 ( 400/ 538- 562/ 624) -26.6 E47-E35 S24-S18
      9  683  353 (S39E33 L=304)  1147    235        1 ( 613/ 271- 763/ 449)  31.3 E46-E23 S46-S31
  

    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