Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1450 1527 (N21W42 L=  9)   275     76        0 (1433/1496-1468/1570)  45.0 W40-W44 N18-N22
      2  895 1192 (N18E03 L=325)   303    115        0 ( 843/1175- 960/1215) -26.6 E05-W01 N16-N20
      3  130 1012 (N25E65 L=263)   619    131        1 ( 105/ 975- 152/1050) -53.1 E67-E62 N22-N27
      4  661  878 (N07E24 L=303)   722    205        1 ( 565/ 838- 756/ 917) -32.5 E30-E18 N03-N11
      5  268  834 (N13E55 L=273)  1494    448        6 ( 162/ 757- 355/ 910)   2.7 E66-E45 N08-N15
      6 1332  719 (S20W09 L=337)   282    106        2 (1315/ 666-1343/ 774)  56.3 W07-W11 S22-S16
      7  841  418 (S24E29 L=298)  1388    370        1 ( 812/ 287- 877/ 553)  45.0 E38-E22 S32-S16
      8  310  447 (S11E72 L=256)  1801    156        2 ( 279/ 402- 349/ 498) -90.0 E75-E67 S14-S07
      9  382  428 (S13E63 L=264)   842    182        2 ( 353/ 391- 411/ 465) -80.5 E66-E60 S16-S10
     10 1181  318 (S42E14 L=314)   289     91        1 (1154/ 294-1212/ 354)  80.5 E14-E11 S44-S38
  

    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