Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  747 1681 (N25E37 L= 95)  1649    322        3 ( 716/1568- 790/1833) -53.1 E50-E29 N16-N36
      2  372 1570 (N10E61 L= 72)  1155    231        5 ( 336/1478- 412/1635) -61.4 E67-E57 N04-N15
      3 1144 1527 (N25E06 L=126)   768    201        2 (1120/1429-1167/1616) -59.7 E10-E03 N18-N30
      4 1016 1468 (N18E12 L=120)   524    138        2 ( 985/1407-1035/1531) -63.4 E13-E09 N14-N22
      5 1781  968 (N11W50 L=182)   621    164        0 (1767/ 892-1794/1046) -77.5 W47-W52 N06-N15
      6  551  816 (S30E25 L=108)   554    170        0 ( 466/ 797- 629/ 831)  14.0 E31-E19 S32-S28
      7  634  660 (S38E16 L=117)   658    221        2 ( 524/ 638- 736/ 682)  26.6 E25-E08 S42-S34
      8  902  619 (S33W04 L=136)   501    104        1 ( 877/ 560- 921/ 667) -50.2 W01-W06 S36-S30
      9 1696  549 (S13W60 L=193)   655    115        1 (1654/ 509-1743/ 591)  66.8 W57-W62 S16-S09
     10  985  526 (S37W13 L=145)   615    110        2 ( 930/ 505-1043/ 549)  14.0 W08-W16 S38-S34
  

    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