Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1077 1853 (N54W43 L=348)  2743    330        4 ( 941/1777-1207/1924)   0.0 W24-W70 N50-N56
      2 1162 1673 (N38W31 L=336)   317     47        1 (1134/1653-1189/1693) -14.0 W28-W33 N36-N39
      3  545 1417 (N39E22 L=283)   309     74        0 ( 516/1389- 578/1447)  18.4 E25-E19 N38-N41
      4 1013 1380 (N25W09 L=314)  1824    426        3 ( 791/1341-1199/1405) -32.5 E04-W19 N17-N31
      5  685 1275 (N27E14 L=291)   463    159        0 ( 659/1205- 734/1331)  40.6 E16-E09 N22-N29
      6  814 1254 (N22E07 L=299)   291    133        0 ( 746/1226- 871/1281) -40.6 E10-E03 N19-N25
      7 1765 1119 (S10W50 L=355)  1496    264        4 (1716/1001-1804/1247)  75.1 W46-W53 S17-S02
      8  130 1109 (N29E71 L=235)   675    105        1 (  97/1068- 149/1135) -17.1 E80-E66 N26-N30
      9  791 1071 (N13E12 L=293)   332    153        0 ( 736/1015- 847/1114) -63.4 E14-E09 N08-N16
     10  221  845 (N12E60 L=245)   492     83        1 ( 203/ 809- 244/ 885)  38.7 E62-E57 N09-N14
     11  872  641 (S14E18 L=287)   282    116        0 ( 830/ 604- 902/ 686)  23.2 E22-E15 S14-S11
     12 1354  621 (S27W10 L=315)   846    214        1 (1271/ 568-1439/ 659)   4.4 W02-W15 S28-S25
     13 1629  607 (S37W31 L=336)   868    148        0 (1561/ 577-1684/ 649)   0.0 W24-W36 S39-S35
     14 1226  458 (S35E03 L=302)   246     35        1 (1209/ 434-1247/ 481)  33.7 E04-E01 S36-S33
     15  337  416 (S16E73 L=232)   907     86        1 ( 302/ 379- 363/ 448) -63.4 E77-E69 S18-S13
     16  978  342 (S36E24 L=281)   494     67        2 ( 951/ 315-1012/ 377) -90.0 E26-E22 S38-S32
     17  790  256 (S37E43 L=262)  1719    272        1 ( 752/ 160- 833/ 353)  26.6 E58-E32 S45-S29
  

    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