Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1166 1798 (N52W43 L= 61)   372     52        1 (1143/1775-1196/1819)  12.5 W39-W48 N50-N52
      2  237 1264 (N36E59 L=318)   732     98        1 ( 205/1221- 270/1305) -28.6 E66-E53 N32-N38
      3  787 1066 (N14E13 L=  4)  1052    259        2 ( 758/ 975- 838/1204)  74.1 E15-E10 N08-N22
      4  357 1075 (N22E44 L=333)   433    151        1 ( 296/1046- 408/1131) -24.4 E50-E39 N20-N26
      5  905  949 (N05E08 L=  9)   401    135        1 ( 842/ 928- 972/ 976) -35.5 E11-E04 N03-N08
      6  416  891 (N11E42 L=335)   631    157        5 ( 371/ 819- 460/ 936) -81.9 E45-E38 N06-N13
      7 1526  854 (S14W27 L= 44)   323    102        0 (1503/ 804-1546/ 896)  45.0 W24-W28 S17-S12
      8  201  819 (N09E65 L=312)   816     84        2 ( 166/ 792- 239/ 855) -29.7 E70-E60 N06-N11
      9 1102  582 (S21E05 L= 12)   553    157        1 (1057/ 517-1143/ 651)  41.2 E09-E01 S24-S16
     10  734  326 (S31E40 L=338)   451     91        0 ( 715/ 283- 759/ 353)  56.3 E42-E36 S34-S28
     11  588  273 (S33E58 L=319)  1181    128        2 ( 546/ 218- 616/ 328)  48.0 E65-E53 S37-S27
  

    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