Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  423 1648 (N23E61 L= 89)   947    236        1 ( 405/1597- 441/1701) -31.0 E67-E57 N20-N26
      2  346 1264 (S02E50 L=100)   438    129        0 ( 319/1209- 380/1313) -26.6 E54-E45 S04-N01
      3 1594 1195 (N20W33 L=183)   721    136        4 (1551/1141-1631/1255) -35.5 W29-W36 N17-N23
      4 1837 1118 (N22W55 L=206)   462     68        0 (1811/1098-1861/1138)  18.4 W52-W58 N20-N23
      5  291  928 (S23E48 L=102)   317    108        1 ( 257/ 896- 319/ 961)  45.0 E51-E45 S25-S20
      6  632  899 (S17E21 L=129)   303    132        0 ( 597/ 871- 678/ 923)   9.5 E23-E17 S19-S16
      7  375  869 (S25E40 L=110)  1020    343        1 ( 300/ 780- 489/ 916)  36.3 E46-E31 S31-S20
      8 1159  825 (S11W12 L=162)   256     88        1 (1110/ 811-1202/ 841)  38.7 W09-W14 S12-S08
      9  929  711 (S22W01 L=151)   406    133        0 ( 868/ 678- 975/ 766)  60.3 E02-W03 S25-S18
     10  960  477 (S37W10 L=160)   630    214        0 ( 844/ 462-1056/ 496)  23.2 W01-W16 S40-S33
     11 1430  439 (S28W43 L=193)   998    184        0 (1341/ 410-1494/ 489)  41.2 W36-W46 S31-S23
  

    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