Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  446 1582 (N26E56 L=241)  1089    170        1 ( 413/1501- 486/1642) -68.2 E59-E53 N20-N30
      2 1414 1587 (N43W21 L=319)   480     61        2 (1384/1554-1443/1612) -78.7 W18-W24 N40-N45
      3  662 1488 (N23E33 L=265)   358    121        1 ( 610/1460- 703/1506)   8.1 E35-E28 N21-N24
      4 1699 1301 (N27W46 L=343)   389     82        0 (1676/1258-1726/1338) -68.2 W43-W48 N24-N29
      5 1419 1194 (N16W23 L=320)   521    178        0 (1363/1140-1475/1262) -36.9 W18-W27 N13-N19
      6  143 1035 (S12E71 L=226)  2930    176        4 ( 116/ 957- 183/1114)  68.2 E77-E64 S17-S07
      7  308  969 (S14E49 L=248)   498     84        0 ( 289/ 923- 326/1010) -78.7 E51-E47 S16-S11
      8  224  793 (S26E60 L=237)  1368    253        3 ( 161/ 720- 276/ 877) -36.0 E69-E53 S31-S21
      9 1554  775 (S07W38 L=336)   250     89        0 (1513/ 752-1595/ 797)  -7.1 W34-W42 S08-S06
     10  630  752 (S22E22 L=276)   269     74        0 ( 609/ 727- 663/ 782) -26.6 E23-E19 S23-S20
     11 1548  654 (S15W41 L=339)  1267    294        4 (1479/ 572-1596/ 731) -45.0 W34-W46 S20-S10
     12 1613  587 (S18W50 L=347)  1736    452        3 (1496/ 477-1683/ 730)  80.0 W42-W55 S26-S08
     13 1720  609 (S15W59 L=357)   774    142        1 (1661/ 573-1768/ 659)  36.9 W54-W63 S18-S11
     14  925  551 (S31W01 L=298)  1573    440        3 ( 744/ 502-1130/ 614)  28.4 E11-W14 S37-S24
     15 1238  303 (S44W34 L=332)  1655    251        1 (1137/ 273-1360/ 343)  26.6 W25-W43 S48-S38
     16  895  314 (S51W04 L=301)   432     86        1 ( 846/ 298- 942/ 331)   0.0 W00-W08 S52-S48
  

    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