Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  323 1591 (N10E69 L= 24)  1506    300        4 ( 285/1524- 345/1640) -38.7 E77-E64 N05-N13
      2  352 1515 (N05E60 L= 33)  1351    236        4 ( 288/1435- 397/1591)  90.0 E64-E57 N00-N10
      3 1037 1463 (N18E12 L= 82)   682    217        1 (1014/1366-1057/1567) -65.2 E14-E08 N12-N25
      4 1404 1476 (N31W10 L=103)   530     94        1 (1372/1425-1432/1508) -31.0 W07-W12 N27-N33
      5 1230 1326 (N16W03 L= 97)   501    157        1 (1200/1249-1252/1406) -58.0 W00-W05 N11-N19
      6  236 1292 (S10E61 L= 32)   810    163        5 ( 212/1239- 263/1345) -81.9 E64-E58 S13-S06
      7 1753 1104 (N20W43 L=137)   458    102        1 (1719/1081-1791/1130)  33.7 W40-W46 N17-N22
      8 1028 1039 (S06E00 L= 93)   358     71        1 (1007/ 992-1055/1075)  90.0 E02-W00 S08-S03
      9  385  869 (S32E38 L= 56)  2931    649        9 ( 202/ 752- 532/1026)  43.8 E55-E28 S43-S19
     10  748  719 (S31E08 L= 86)   927    285        0 ( 669/ 624- 850/ 832)  67.6 E11-E04 S39-S22
     11 1787  771 (N03W57 L=151)  1120    327        2 (1677/ 737-1867/ 832)  32.0 W47-W64 S01-N09
     12 1340  724 (S14W27 L=120)   452    192        1 (1287/ 671-1414/ 799) -18.4 W21-W33 S15-S10
     13 1080  720 (S22W12 L=106)   804    235        2 ( 956/ 672-1171/ 783)  53.1 W06-W15 S27-S15
     14 1527  474 (S21W50 L=144)  3627    760       14 (1271/ 333-1784/ 687)  58.3 W39-W61 S35-S01
     15 1235  619 (S23W25 L=118)   781    295        1 (1080/ 551-1374/ 680)  34.6 W18-W30 S30-S15
     16 1561  354 (S25W61 L=155)   565     84        1 (1534/ 323-1598/ 377)  78.7 W58-W62 S27-S22
  

    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