Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1535 1723 (N30W62 L= 73)  1769    191        1 (1478/1672-1635/1752)  -6.8 W53-W78 N27-N32
      2  936 1634 (N41W11 L= 23)   251     47        1 ( 913/1614- 961/1657)  26.6 W09-W13 N40-N42
      3 1893 1209 (S08W66 L= 78)  2725    180        7 (1846/1141-1951/1273) -45.0 W61-W75 S12-S03
      4 1790 1232 (S04W55 L= 66)   382     93        1 (1762/1203-1832/1263) -45.0 W52-W59 S06-S01
      5 1554 1185 (S00W35 L= 46)   326    118        0 (1514/1138-1600/1210)   8.1 W31-W38 S01-N00
      6 1798 1168 (S07W54 L= 65)   328     67        1 (1780/1136-1811/1206)  31.0 W51-W56 S09-S05
      7  193 1140 (N27E58 L=313)   480     99        1 ( 158/1108- 238/1160)   0.0 E63-E52 N25-N28
      8 1293 1006 (S04W15 L= 26)  2029    565        2 (1119/ 827-1443/1137) -56.3 W06-W20 S17-N05
      9 1545 1108 (S04W32 L= 44)   267     42        1 (1519/1087-1575/1126) -14.0 W30-W34 S05-S02
     10  393 1109 (N21E38 L=333)   232     44        0 ( 365/1092- 423/1122) -26.6 E40-E35 N20-N22
     11 1647 1014 (S12W37 L= 49)  1383    266        5 (1566/ 920-1705/1103)  21.0 W30-W43 S16-S07
     12 1582 1047 (S08W33 L= 45)   382     63        0 (1557/1007-1601/1079)  76.0 W32-W35 S10-S06
     13  724 1029 (N10E17 L=355)   560    167        1 ( 634/1003- 796/1058) -36.9 E21-E12 N07-N13
     14 1728  794 (S27W42 L= 53)  1691    281        4 (1635/ 759-1853/ 858) -10.3 W33-W55 S30-S21
     15 1540  603 (S34W24 L= 35)  2703    514        9 (1451/ 486-1609/ 807)  84.8 W16-W30 S42-S20
     16  883  708 (S11E15 L=356)   630    144        2 ( 830/ 677- 943/ 744) -45.0 E17-E11 S13-S08
     17  734  460 (S23E33 L=339)  3719    806        8 ( 501/ 245- 984/ 681)  14.7 E67-E10 S33-S14
     18  473  583 (S10E45 L=326)   372    155        0 ( 397/ 561- 519/ 618)   0.0 E52-E41 S10-S08
  

    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