Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  765 1722 (N33E39 L= 93)   508     91        1 ( 735/1689- 789/1765) -71.6 E41-E37 N30-N36
      2 1218 1670 (N42E06 L=127)  1097    178        1 (1165/1590-1273/1738) -42.0 E10-W00 N37-N47
      3  522 1537 (N15E46 L= 86)   822    177        2 ( 485/1494- 564/1613) -33.7 E51-E41 N13-N19
      4  857 1540 (N23E24 L=108)   262     49        1 ( 836/1518- 873/1569)  90.0 E24-E22 N21-N25
      5  594 1495 (N14E39 L= 93)   265     64        1 ( 579/1463- 621/1516)  90.0 E40-E37 N12-N16
      6  781 1432 (N15E25 L=108)   310     71        1 ( 767/1393- 800/1469) -45.0 E26-E22 N12-N16
      7 1094 1310 (N15E03 L=129)   366     95        2 (1063/1280-1129/1351) -26.6 E04-W00 N12-N17
      8  182 1250 (S08E68 L= 64)  1448    300        6 ( 139/1199- 223/1297)  54.5 E76-E62 S11-S04
      9  986  839 (S15W02 L=134)   224     49        1 ( 955/ 816-1006/ 858)   0.0 W00-W03 S15-S13
     10 1719  791 (N01W51 L=184)   467     95        1 (1681/ 746-1754/ 829)  59.0 W49-W53 S02-N03
     11  650  638 (S35E15 L=117)   256     45        0 ( 631/ 614- 675/ 666) -36.9 E16-E12 S36-S33
     12 1264  311 (S40W41 L=173)  1265    293        1 (1157/ 282-1402/ 340)  14.0 W32-W52 S42-S36
  

    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