Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  737 1678 (N25E38 L= 96)  1598    295        4 ( 701/1572- 776/1835) -53.1 E50-E30 N17-N37
      2  376 1595 (N12E62 L= 71)  1048    141        4 ( 343/1528- 407/1661) -59.0 E66-E59 N07-N17
      3 1008 1510 (N20E14 L=119)   556    243        0 ( 989/1420-1038/1634) -66.8 E17-E10 N14-N29
      4 1134 1531 (N25E07 L=126)   710    212        2 (1107/1439-1157/1628) -57.5 E11-E04 N18-N30
      5 1774  973 (N11W49 L=182)   550    159        0 (1760/ 894-1787/1052) -71.6 W47-W50 N06-N15
      6  543  819 (S30E26 L=108)   571    171        0 ( 453/ 801- 617/ 834)   8.7 E32-E19 S32-S28
      7  638  665 (S38E16 L=118)   664    263        1 ( 514/ 643- 761/ 689)  25.3 E26-E07 S42-S32
      8  895  626 (S33W03 L=137)   396     90        1 ( 879/ 569- 912/ 670) -56.3 W01-W05 S36-S30
      9 1689  550 (S13W60 L=193)   635    137        0 (1643/ 509-1736/ 594)  66.8 W57-W62 S16-S09
     10  984  532 (S36W12 L=146)   635    144        2 ( 922/ 512-1044/ 561)  32.0 W07-W15 S38-S32
  

    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