Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1355 1592 (N31W31 L= 58)   416    137        0 (1290/1554-1413/1618) -32.0 W26-W34 N27-N32
      2 1709 1451 (N18W55 L= 82)   321     39        1 (1684/1433-1727/1471)  11.3 W52-W57 N16-N19
      3  875 1403 (N24E05 L= 22)   424    163        0 ( 810/1347- 946/1432) -36.9 E09-E01 N19-N26
      4 1784  818 (S21W53 L= 80)   353     74        0 (1769/ 776-1799/ 860)  90.0 W51-W53 S23-S18
      5 1717  650 (S31W49 L= 76)   396     22        1 (1695/ 627-1737/ 674) -45.0 W46-W51 S32-S29
      6  287  501 (S24E69 L=319)   994    129        3 ( 243/ 452- 309/ 558)  37.9 E76-E63 S27-S20
      7 1619  406 (S49W49 L= 76)   396     26        0 (1595/ 389-1641/ 422) -20.6 W44-W52 S50-S46
      8 1565  375 (S50W43 L= 71)   591     69        0 (1529/ 353-1596/ 392) -17.1 W37-W50 S52-S48
      9 1372  328 (S52W22 L= 49)   876    103        1 (1316/ 292-1417/ 379) -21.0 W15-W28 S55-S47
     10 1464  332 (S53W33 L= 60)   516     82        1 (1434/ 298-1493/ 358) -33.7 W28-W37 S56-S49
  

    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