Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1554 1470 (N26W45 L=137)   410    126        0 (1501/1439-1606/1503) -45.0 W42-W49 N23-N29
      2 1773 1400 (N18W65 L=156)   704     69        1 (1736/1366-1800/1420) -63.4 W60-W69 N15-N19
      3 1453 1302 (N17W32 L=124)   268    118        1 (1402/1272-1503/1335) -39.8 W28-W34 N15-N20
      4  571 1308 (N27E29 L= 63)   490    150        0 ( 516/1261- 617/1331) -20.6 E33-E25 N24-N29
      5 1277 1169 (N12W17 L=109)   194     24        1 (1257/1148-1295/1190)  45.0 W16-W18 N10-N12
      6  175 1090 (N15E67 L= 25)  1663    234        4 ( 135/1012- 200/1177) -73.3 E75-E62 N10-N20
      7  992  887 (S02E03 L= 88)  1484    384        3 ( 939/ 726-1042/1063)  87.1 E05-E01 S12-N08
      8 1754  855 (S15W50 L=142)   563     27        1 (1729/ 832-1780/ 883) -63.4 W47-W52 S16-S12
      9  195  812 (S02E66 L= 25)   527     68        1 ( 173/ 779- 220/ 841) -31.0 E69-E63 S04-N00
     10  589  635 (S15E33 L= 59)   412    146        1 ( 541/ 609- 647/ 655) -15.9 E37-E29 S16-S13
  

    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