Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  760 1657 (N37E22 L=257)  1616    316        2 ( 632/1583- 897/1726) -23.7 E35-E10 N31-N42
      2 1564 1577 (N34W41 L=320)   611    164        2 (1494/1523-1628/1617) -33.7 W35-W45 N29-N36
      3 1910 1095 (N04W67 L=346)   554     74        0 (1890/1055-1925/1120)  59.0 W64-W70 N01-N06
      4 1728  884 (S10W48 L=327)   331     60        0 (1708/ 855-1749/ 925) -36.9 W46-W50 S11-S07
      5 1556  807 (S15W35 L=314)   655    177        0 (1523/ 729-1579/ 877) -83.7 W33-W37 S19-S10
      6  320  811 (S17E49 L=230)   689      0       -1 ( 279/ 779- 351/ 846)  -8.1 E53-E46 S18-S14
      7  218  733 (S21E63 L=216)   607     69        2 ( 189/ 708- 254/ 766)   5.7 E67-E57 S22-S18
      8 1709  645 (S24W52 L=332)   627    138        0 (1661/ 591-1742/ 702)  74.1 W49-W55 S27-S20
      9  636  601 (S30E27 L=252)   633    156        2 ( 592/ 550- 672/ 662)  66.8 E31-E23 S33-S26
     10  414  566 (S32E47 L=232)   817    178        2 ( 337/ 512- 488/ 601)  14.0 E59-E39 S35-S30
     11 1331  215 (S59W42 L=321)   480     44        1 (1305/ 199-1352/ 234)  26.6 W37-W47 S61-S57
     12 1185  190 (S63W25 L=304)   764     69        1 (1146/ 161-1223/ 224)  18.4 W21-W31 S66-S59
  

    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