Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1180 1489 (N26W15 L= 50)   988    234        1 (1062/1437-1275/1518) -24.8 W07-W20 N21-N28
      2  497 1346 (N24E32 L=  2)   966    216        3 ( 409/1288- 599/1390) -26.6 E40-E24 N18-N27
      3 1022 1155 (N07W01 L= 36)   272     80        0 ( 989/1120-1056/1195) -68.2 W00-W02 N04-N09
      4  786 1075 (N04E14 L= 21)   251     98        0 ( 772/1022- 801/1119)  78.7 E14-E12 N01-N06
      5 1894  967 (S12W65 L= 99)  1553    260        7 (1853/ 881-1942/1091) -49.8 W59-W74 S17-S04
      6 1695  931 (S13W44 L= 78)   291     70        1 (1660/ 907-1726/ 948)   0.0 W40-W46 S14-S11
      7 1492  724 (S24W28 L= 62)   619    214        1 (1466/ 623-1519/ 808)  85.2 W26-W29 S30-S18
      8 1807  763 (S24W56 L= 91)   376     54        1 (1775/ 742-1837/ 779)  -7.1 W52-W60 S25-S22
      9  629  520 (S28E34 L=  1)  1352    257        1 ( 559/ 412- 695/ 628)  35.8 E44-E26 S34-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