Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  796 1729 (N59E03 L=  1)  1849    271        5 ( 704/1652- 911/1787)   4.2 E14-W13 N53-N62
      2 1252 1724 (N46W39 L= 44)  1097    178        1 (1158/1697-1333/1749) -10.0 W30-W47 N43-N48
      3 1468 1697 (N38W56 L= 60)   890    116        2 (1422/1665-1526/1722) -24.0 W50-W62 N35-N40
      4 1638 1463 (N20W54 L= 58)   398     96        2 (1604/1428-1681/1493) -45.0 W51-W57 N17-N22
      5  656 1347 (N32E20 L=345)   242     83        0 ( 630/1321- 696/1365)  11.3 E21-E16 N31-N33
      6 1387 1266 (N14W27 L= 31)  1111    257        3 (1336/1160-1422/1360)  80.5 W23-W29 N08-N20
      7  642 1094 (N17E23 L=341)   414    133        0 ( 584/1066- 709/1123)   6.3 E27-E18 N15-N17
      8  520  650 (S08E39 L=325)   260     58        1 ( 484/ 629- 553/ 669) -21.8 E42-E37 S09-S06
      9 1411  587 (S28W18 L= 22)   295     70        1 (1395/ 556-1434/ 620) -78.7 W16-W19 S30-S25
     10 1009  242 (S47E21 L=343)   873    143        1 ( 953/ 206-1078/ 289)  19.7 E27-E13 S49-S43
  

    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