Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1048 1272 (N19W08 L=325)   481    113        0 ( 992/1251-1116/1299) -29.7 W05-W12 N16-N21
      2  164 1104 (N30E60 L=257)  2111    362        2 ( 106/1009- 216/1193) -52.4 E65-E55 N23-N36
      3 1862  985 (S22W55 L= 12)   394     70        0 (1842/ 955-1884/1017)  21.8 W52-W58 S23-S20
      4  801  953 (N07E14 L=302)   857    270        2 ( 672/ 910- 930/ 987) -29.7 E21-E06 N03-N12
      5  406  896 (N13E41 L=275)  1086    377        4 ( 257/ 832- 517/ 957)  -2.5 E55-E32 N11-N15
      6 1274  866 (S10W10 L=326)   347    101        0 (1223/ 836-1313/ 907)   9.5 W06-W12 S10-S08
      7 1476  797 (S20W20 L=337)   388     93        1 (1449/ 730-1501/ 842)  51.3 W16-W22 S22-S17
      8  187  779 (N11E66 L=250)   496     68        0 ( 168/ 751- 211/ 818) -68.2 E67-E64 N08-N13
      9  984  493 (S24E18 L=299)  1402    370        3 ( 947/ 350-1035/ 638)  45.0 E26-E10 S32-S16
     10  360  455 (S11E64 L=252)  7151    709       15 ( 281/ 350- 466/ 551) -60.3 E74-E56 S19-S04
     11  479  445 (S14E53 L=264)   487     76        1 ( 454/ 419- 503/ 468)  18.4 E55-E49 S16-S12
  

    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