Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  680 1612 (N18E39 L=167)  1129    308        1 ( 640/1469- 710/1753) -41.8 E50-E31 N10-N27
      2  320 1556 (N07E67 L=140)  6265    325        6 ( 226/1435- 401/1680)  82.9 E75-E59 N00-N16
      3  939 1619 (N26E23 L=183)   351    100        1 ( 918/1567- 963/1662) -60.3 E25-E20 N22-N29
      4  278 1133 (S18E52 L=154)  1909    305        1 ( 189/ 996- 353/1234)  52.6 E60-E47 S27-S10
      5  781 1025 (S12E13 L=193)   212     88        1 ( 738/1003- 820/1044)  21.8 E16-E11 S13-S11
      6  396  919 (S28E38 L=169)  2907    568        4 ( 148/ 776- 633/1014)  25.6 E71-E21 S41-S16
      7 1254  880 (S07W17 L=223)   470    224        0 (1205/ 804-1337/ 927)  63.4 W15-W21 S12-S03
      8  981  702 (S25W07 L=213)   286     71        0 ( 937/ 686-1021/ 718)  45.0 W04-W08 S26-S22
      9  423  575 (S50E33 L=173)   815    153        2 ( 377/ 512- 486/ 633) -24.4 E38-E25 S52-S45
     10 1687  484 (S14W64 L=271)  3215    297        4 (1601/ 420-1747/ 544)   8.5 W54-W75 S18-S09
     11  786  403 (S50W06 L=212)   326     69        0 ( 767/ 369- 810/ 439) -33.7 W02-W09 S51-S47
     12  972  329 (S48W24 L=230)  1574    267        4 ( 855/ 268-1082/ 384)  13.4 W14-W36 S51-S44
  

    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