Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1319 1759 (N56W18 L= 85)   815    125        1 (1253/1737-1385/1779)  -3.6 W10-W26 N54-N58
      2  716 1526 (N29E28 L= 39)  1150    249        3 ( 681/1399- 772/1603) -81.9 E31-E24 N20-N34
      3  402 1533 (N25E56 L= 12)   583     97        0 ( 374/1489- 443/1566)  59.0 E57-E52 N21-N27
      4 1708 1152 (N16W46 L=113)   453    123        0 (1670/1101-1732/1206) -50.2 W42-W48 N12-N18
      5  453 1157 (N03E38 L= 29)   752    194        0 ( 380/1088- 502/1203) -16.7 E44-E33 S01-N05
      6 1526  955 (N02W32 L=100)   563     88        1 (1492/ 907-1546/ 998) -76.0 W29-W34 N00-N04
      7  903  921 (S07E06 L= 61)   398    131        0 ( 862/ 879- 953/ 967)  50.2 E08-E03 S09-S03
      8 1773  676 (S13W60 L=128)  1276    146        4 (1738/ 609-1821/ 737) -82.9 W56-W64 S16-S08
      9 1650  590 (S19W51 L=118)   510     98        1 (1611/ 560-1695/ 621)  45.0 W46-W54 S21-S16
     10 1695  495 (S25W62 L=129)   546     65        0 (1665/ 472-1729/ 520)  53.1 W59-W65 S26-S22
     11 1429  318 (S41W47 L=114)  1492    178        1 (1360/ 280-1515/ 369)  38.7 W41-W52 S44-S35
  

    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