Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  959 1641 (N32E13 L= 49)   299     73        0 ( 919/1618- 995/1659)  -9.5 E16-E10 N30-N32
      2  311 1488 (N16E58 L=  4)  6774    484       15 ( 207/1342- 442/1608)  45.0 E75-E47 N06-N24
      3  179 1405 (N11E72 L=350)   962     87        2 ( 160/1356- 199/1444) -78.7 E76-E68 N07-N13
      4  326 1386 (N09E52 L=  9)   339     67        0 ( 295/1371- 355/1406)   9.5 E55-E49 N08-N10
      5 1705 1153 (N11W42 L=104)   807    211        2 (1670/1066-1746/1259)  72.9 W39-W45 N05-N18
      6  423 1229 (N00E40 L= 21)   434    107        2 ( 378/1199- 467/1252)   0.0 E43-E37 S01-N01
      7 1849 1139 (N13W56 L=118)  1180    153        2 (1802/1069-1880/1205) -54.5 W51-W60 N09-N16
      8  484  822 (S23E33 L= 28)  2437    709        4 ( 302/ 642- 683/1146)  46.8 E53-E21 S35-S02
      9  166 1058 (S11E63 L=359)  2119    264        6 ( 106/ 956- 256/1117)  27.9 E72-E53 S16-S07
     10  390 1003 (S13E41 L= 21)   392     92        2 ( 347/ 970- 423/1032) -18.4 E44-E38 S14-S11
     11  242  964 (S16E54 L=  8)   332     25        0 ( 222/ 942- 265/ 981)   0.0 E56-E51 S17-S15
     12 1389  904 (S08W23 L= 85)   397    138        0 (1354/ 837-1413/ 969)  82.9 W22-W24 S12-S04
     13  757  919 (S15E15 L= 47)   304    109        0 ( 724/ 872- 785/ 965) -45.0 E17-E12 S17-S12
     14  206  872 (S22E59 L=  3)  2190    342        3 (  86/ 824- 374/ 904)   5.4 E83-E41 S24-S19
     15 1887  839 (S03W66 L=127)   547     72        1 (1867/ 803-1907/ 869) -31.0 W62-W69 S04-S01
     16 1434  776 (S15W29 L= 90)   238     67        0 (1411/ 747-1454/ 805)  63.4 W27-W29 S16-S12
     17 1323  713 (S20W23 L= 85)   648    235        2 (1249/ 676-1408/ 779)  -5.2 W17-W28 S22-S17
     18  977  442 (S42W06 L= 67)  2612    487        7 ( 753/ 402-1195/ 487)   8.6 E12-W22 S45-S37
  

    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