Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  440 1323 (N14E44 L=244)   406     74        1 ( 415/1288- 465/1361) -51.3 E46-E40 N11-N16
      2  721 1257 (N13E22 L=266)   460    142        1 ( 663/1228- 782/1276)   0.0 E25-E17 N10-N14
      3  850 1031 (N00E11 L=277)   257     49        1 ( 830/1004- 873/1055) -63.4 E12-E09 S01-N01
      4  879  862 (S10E07 L=280)   389     86        2 ( 851/ 830- 915/ 885)  45.0 E09-E05 S11-S08
      5  627  789 (S17E24 L=264)   244     53        1 ( 605/ 766- 647/ 821) -76.0 E25-E22 S18-S14
      6  911  748 (S17E04 L=283)   613    185        2 ( 831/ 725- 981/ 796)  -6.3 E09-W00 S18-S14
      7 1572  715 (S13W40 L=328)   908    187        0 (1541/ 640-1620/ 788)  78.7 W38-W42 S17-S07
      8  234  741 (S25E61 L=227)   759    105        1 ( 187/ 719- 279/ 775)  14.9 E69-E54 S26-S21
      9 1333  680 (S17W24 L=311)   389    120        0 (1279/ 654-1376/ 724)  45.0 W20-W25 S18-S13
     10 1716  641 (S16W56 L=344)   456     57        2 (1693/ 613-1734/ 671) -76.0 W53-W58 S18-S14
     11  763  537 (S33E14 L=274)   286     80        0 ( 720/ 527- 805/ 550)   0.0 E16-E10 S33-S31
     12 1511  360 (S38W53 L=340)  1100    118        3 (1476/ 319-1549/ 408) -69.4 W47-W56 S41-S33
  

    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