Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1202 1707 (N40W30 L=111)   639    119        0 (1163/1668-1246/1741)  24.0 W25-W34 N37-N41
      2  819 1654 (N45E01 L= 80)   259     51        0 ( 801/1636- 843/1678) -71.6 E02-W00 N43-N47
      3 1641 1581 (N22W56 L=137)   761    133        1 (1597/1544-1675/1610)   6.3 W51-W60 N19-N23
      4  922 1570 (N36W05 L= 86)   276     48        0 ( 894/1553- 951/1587) -33.7 W02-W07 N35-N37
      5  754 1366 (N26E10 L= 70)   492    197        0 ( 658/1329- 844/1406) -30.3 E16-E04 N23-N30
      6 1674 1019 (S11W41 L=122)  1944    472        0 (1579/ 815-1825/1166)  31.0 W31-W57 S21-S04
      7 1878 1082 (S11W61 L=142)   599     87        0 (1846/1041-1905/1114) -51.3 W58-W64 S14-S09
      8  799  290 (S40E35 L= 45)   443     98        0 ( 780/ 255- 819/ 327)  50.2 E37-E31 S43-S37
  

    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