Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1450 1741 (N32W56 L=101)   895    206        2 (1394/1700-1509/1765) -15.3 W49-W61 N29-N34
      2 1053 1253 (N18W08 L= 53)   697    270        1 ( 939/1231-1213/1277) -26.6 W01-W17 N13-N21
      3 1534 1225 (N02W35 L= 80)   401    107        3 (1498/1194-1562/1248) -45.0 W33-W37 N00-N04
      4 1891 1140 (S14W63 L=108)  3149    447        8 (1852/1014-1949/1204)  34.7 W55-W73 S20-S09
      5 1697  818 (S25W36 L= 82)  2316    498        2 (1571/ 684-1762/1070)  43.5 W25-W46 S32-S12
      6 1177  852 (S07W04 L= 49)  1064    391        0 (1049/ 782-1332/ 963) -32.7 E02-W13 S11-N01
      7 1574  818 (S21W27 L= 72)   289     79        1 (1554/ 761-1594/ 853)  51.3 W24-W28 S24-S19
      8 1276  651 (S22W04 L= 49)   947    398        3 (1240/ 465-1340/ 821)  53.6 E03-W12 S33-S13
      9  293  668 (N03E59 L=346)  1322    308        4 ( 210/ 627- 373/ 717) -10.0 E68-E50 N00-N06
     10  903  498 (S21E22 L= 23)   238     75        0 ( 884/ 460- 921/ 521)  78.7 E24-E20 S23-S18
     11  944  337 (S32E27 L= 18)   318     75        1 ( 912/ 318- 967/ 362)   9.5 E30-E24 S34-S31
  

    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