Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1629 1547 (N21W59 L=187)   564    104        0 (1580/1529-1686/1565) -32.0 W54-W63 N17-N22
      2 1259 1492 (N28W28 L=156)   421     80        3 (1220/1461-1291/1531)  23.2 W24-W31 N26-N30
      3 1636 1473 (N17W54 L=182)   373     60        0 (1609/1449-1663/1501) -63.4 W53-W56 N14-N18
      4  200 1050 (N23E62 L= 66)  2071    214        5 ( 150/ 952- 248/1119) -37.6 E70-E56 N17-N28
      5 1853  940 (S20W57 L=185)   664     65        2 (1823/ 914-1880/ 967) -29.7 W53-W61 S22-S18
      6 1640  867 (S17W35 L=163)   692    146        3 (1590/ 806-1694/ 922)  21.8 W30-W40 S20-S14
      7  838  489 (S21E25 L=103)  1795    424        3 ( 711/ 345- 943/ 639)  31.0 E38-E13 S28-S13
      8  952  465 (S25E18 L=110)   458    158        1 ( 893/ 419-1007/ 497) -39.8 E22-E15 S29-S22
      9 1084  440 (S30E11 L=117)   288     77        0 (1046/ 421-1124/ 466)   0.0 E13-E07 S30-S28
     10  618  344 (S27E48 L= 80)  1529    286        2 ( 545/ 300- 716/ 437) -67.4 E53-E39 S31-S18
  

    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