Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1286 1683 (N35W42 L=148)  1395    218        3 (1182/1644-1389/1710) -33.7 W34-W48 N29-N38
      2  373 1396 (N45E39 L= 67)  1489    194        3 ( 286/1332- 456/1441) -17.7 E52-E30 N39-N49
      3 1854 1226 (S08W63 L=169)   840    139        6 (1812/1168-1886/1266) -49.4 W58-W67 S12-S05
      4 1527 1138 (S02W32 L=138)   374     80        1 (1510/1091-1549/1184)  78.7 W30-W33 S05-N00
      5 1802 1026 (S17W50 L=157)  1353    332        0 (1712/ 919-1888/1146)  16.3 W40-W64 S20-S13
      6  149 1057 (N28E67 L= 39)  1811    149        2 (  86/1002- 202/1104) -16.9 E84-E59 N23-N31
      7 1581  662 (S32W25 L=132)  1398    291        3 (1543/ 545-1614/ 800)  77.5 W20-W27 S39-S21
      8  405  744 (N05E46 L= 61)   489    113        0 ( 362/ 713- 452/ 787) -56.3 E48-E43 N02-N08
      9  231  726 (N07E64 L= 43)   470     23        1 ( 208/ 707- 252/ 747)   9.5 E66-E60 N06-N08
     10  896  583 (S16E20 L= 87)   256     63        0 ( 883/ 547- 909/ 618)  63.4 E20-E18 S17-S13
     11  734  425 (S21E36 L= 70)   870    237        1 ( 653/ 357- 800/ 497)  16.4 E45-E28 S24-S18
     12 1374  413 (S42W04 L=111)   366     93        0 (1331/ 393-1418/ 436)  -8.1 W00-W07 S43-S40
  

    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