Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1222 1579 (N35W28 L=170)   294     94        0 (1176/1565-1270/1593) -15.9 W24-W31 N33-N36
      2  583 1389 (N38E23 L=119)   764    161        1 ( 529/1336- 648/1436) -45.0 E27-E18 N33-N42
      3  348 1148 (N26E44 L= 98)  1737    457        2 ( 292/ 954- 414/1310) -70.8 E49-E40 N14-N37
      4 1790 1184 (S02W56 L=197)   634     93        2 (1756/1155-1818/1212) -26.6 W52-W58 S04-N00
      5  823 1183 (N20E09 L=133)   233     66        0 ( 793/1159- 854/1206) -53.1 E10-E07 N18-N22
      6  901 1139 (N16E05 L=137)   289    102        1 ( 857/1120- 964/1155) -18.4 E07-E01 N14-N17
      7  163  978 (N18E65 L= 77)   661     90        1 ( 146/ 926- 182/1018)  80.5 E66-E63 N14-N20
      8 1125  934 (S00W04 L=146)   904    212        2 (1055/ 828-1175/1015) -84.8 W01-W06 S06-N05
      9 1727  941 (S14W44 L=186)   721    108        1 (1708/ 876-1751/ 996)  60.3 W41-W46 S18-S11
     10 1086  651 (S16E04 L=138)   810    140        2 (1014/ 590-1125/ 696) -63.4 E07-E01 S20-S12
     11  868  315 (S35E29 L=112)  1065    217        1 ( 802/ 247- 906/ 377)  32.7 E38-E24 S39-S30
  

    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