Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1493 1731 (N37W55 L=236)   539     62        0 (1456/1712-1527/1751)   4.8 W49-W61 N36-N38
      2  893 1689 (N45W03 L=184)   279     44        0 ( 876/1658- 910/1715)  59.0 W00-W04 N42-N47
      3 1764 1541 (N20W70 L=251)   672     76        1 (1744/1517-1785/1580)  76.0 W66-W73 N18-N22
      4 1493 1494 (N21W40 L=220)   275     68        0 (1463/1472-1521/1514) -56.3 W37-W41 N19-N22
      5  498 1476 (N36E32 L=149)   344     98        0 ( 461/1458- 548/1489) -18.4 E35-E26 N34-N37
      6  394 1439 (N36E41 L=140)   716    107        0 ( 360/1393- 432/1487)  45.0 E44-E37 N32-N38
      7 1380  963 (S08W20 L=201)   515    179        1 (1330/ 893-1450/1002)  18.4 W16-W25 S11-S06
      8 1899  914 (S19W66 L=247)   486     56        0 (1878/ 888-1918/ 946)  45.0 W62-W68 S20-S17
      9  965  676 (S19E09 L=172)   542    245        0 ( 876/ 599-1032/ 768)  31.6 E16-E03 S23-S14
     10  233  613 (S12E68 L=113)  2931    207        9 ( 169/ 527- 307/ 691) -42.3 E78-E62 S18-S07
     11  702  540 (S24E29 L=152)   434    154        1 ( 660/ 492- 722/ 611)  41.2 E33-E25 S26-S19
     12 1214  559 (S31W05 L=186)   273     80        0 (1176/ 535-1250/ 580)  11.3 W02-W07 S31-S29
  

    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