Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1513 1741 (N35W60 L=345)  1145    126        1 (1467/1698-1575/1777)   5.2 W51-W73 N33-N36
      2 1566 1616 (N26W53 L=338)   479    123        0 (1510/1595-1624/1628) -24.0 W48-W57 N23-N27
      3  989 1421 (N25W05 L=290)  1502    405        4 ( 822/1328-1168/1500)   9.5 E06-W18 N22-N27
      4  309 1387 (N36E48 L=237)   617     90        1 ( 265/1364- 359/1406) -14.0 E54-E42 N34-N38
      5  604 1378 (N30E21 L=264)   282     55        1 ( 580/1347- 629/1401)  33.7 E22-E19 N28-N31
      6  711 1370 (N27E14 L=272)   252     79        1 ( 685/1341- 752/1391)  11.3 E15-E10 N25-N28
      7  159 1218 (N28E66 L=220)  1186    134        0 (  87/1210- 230/1225)  -3.5 E88-E55 N26-N29
      8  152 1167 (N25E66 L=219)  1518    147        0 (  78/1157- 226/1177)  -3.5 E89-E55 N23-N26
      9 1954 1094 (S12W76 L=  1)   860     32        0 (1928/1080-1978/1108)   3.6 W00-W88 S13-S10
     10 1634  916 (S16W37 L=322)  1540    375        0 (1574/ 766-1665/1089)  60.9 W31-W41 S24-S06
     11  540 1050 (N11E29 L=256)   202     45        0 ( 518/1034- 565/1070) -33.7 E30-E27 N10-N12
     12 1511  939 (S12W28 L=313)   243     63        0 (1486/ 910-1540/ 966)  26.6 W26-W30 S13-S11
     13 1714  796 (S25W43 L=328)   328     83        0 (1688/ 756-1740/ 835)  31.0 W40-W45 S27-S22
     14  157  790 (N02E70 L=215)   746    120        0 (  97/ 785- 225/ 796)  -4.2 E88-E60 N00-N03
     15  735  639 (S16E25 L=260)   240    122        1 ( 711/ 584- 770/ 685)  39.8 E27-E21 S19-S13
     16  269  469 (S19E78 L=207)  1589    115        1 ( 230/ 430- 309/ 510) -45.0 E82-E72 S21-S15
  

    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