Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  840 1550 (N28E24 L=159)  3371    745       11 ( 643/1358-1115/1763) -19.4 E51-W00 N20-N39
      2  388 1552 (N19E60 L=122)   940    121        3 ( 366/1493- 414/1604) -32.5 E66-E53 N15-N22
      3  395 1359 (N08E49 L=133)   338     40        1 ( 378/1335- 419/1383)  45.0 E50-E47 N06-N08
      4 1786 1244 (N27W53 L=235)  1030    194        4 (1727/1173-1844/1317) -30.3 W46-W59 N23-N30
      5  449 1266 (N03E42 L=141)   632     97        1 ( 404/1236- 489/1302)  38.7 E45-E39 N01-N05
      6 1691 1193 (N22W43 L=225)   587    121        1 (1655/1142-1729/1242) -39.8 W39-W46 N18-N24
      7 1243 1128 (N09W11 L=193)   700    286        2 (1131/1091-1358/1200)  -3.8 W03-W18 N07-N11
      8 1059 1050 (N01W02 L=184)   699    224        2 ( 976/1004-1161/1071)  33.7 E03-W07 S03-N03
      9  427  761 (S28E37 L=145)   356    114        3 ( 403/ 702- 443/ 808)  69.4 E39-E35 S32-S24
     10  931  742 (S20E00 L=182)   544    169        0 ( 891/ 658- 965/ 806)  84.3 E01-W00 S25-S15
     11 1007  639 (S25W06 L=189)   482    145        2 ( 945/ 610-1076/ 663)   6.3 W02-W11 S25-S23
     12  322  620 (S40E51 L=132)   720     87        2 ( 280/ 598- 356/ 651)  18.4 E57-E45 S42-S37
     13 1135  535 (S29W18 L=200)   437    132        3 (1114/ 467-1155/ 593) -71.6 W15-W19 S33-S24
     14 1388  524 (S24W35 L=218)   356    101        2 (1347/ 486-1418/ 565)  90.0 W33-W37 S27-S21
     15 1616  500 (S21W56 L=238)   391    146        1 (1553/ 459-1660/ 549)  60.3 W51-W57 S24-S17
     16  859  364 (S48W03 L=186)   353     72        2 ( 831/ 331- 886/ 393) -38.7 W00-W06 S49-S45
  

    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