Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1521 (N28W21 L=112)   359     45        0 (1312/1495-1378/1548) -45.0 W18-W23 N26-N30
      2  333 1447 (N22E51 L= 40)   444     89        1 ( 313/1412- 358/1480) -53.1 E53-E48 N20-N24
      3 1389 1448 (N23W23 L=115)   323     58        1 (1354/1426-1416/1468)  21.8 W20-W25 N22-N24
      4  521 1439 (N21E35 L= 56)   323     96        0 ( 472/1425- 562/1451)  -8.1 E38-E31 N20-N22
      5  283 1193 (N06E52 L= 40)   511    142        1 ( 264/1140- 311/1273) -63.4 E53-E48 N03-N11
      6  148 1162 (N05E67 L= 24)  1408    225        3 ( 116/1094- 183/1222)  45.0 E72-E62 N01-N08
      7  960 1125 (N02E04 L= 87)   514    141        2 ( 923/1082- 986/1196)  71.6 E06-E02 N00-N06
      8  985 1033 (S03E02 L= 89)   174    114        0 ( 961/ 980-1035/1062)  59.0 E03-W00 S06-S01
      9 1208 1016 (S04W11 L=102)   281     74        3 (1155/ 994-1249/1035)   0.0 W07-W13 S05-S02
     10 1292  877 (S12W17 L=108)   192     92        0 (1246/ 866-1342/ 893)  18.4 W13-W20 S12-S10
     11 1770  860 (S11W52 L=144)   414     93        2 (1740/ 822-1794/ 886) -31.0 W49-W54 S12-S09
     12 1364  773 (S18W22 L=114)   694    193        1 (1272/ 734-1426/ 852)  33.7 W16-W25 S20-S12
     13  918  737 (S21E06 L= 85)   718    178        1 ( 887/ 640- 946/ 821)  84.8 E08-E04 S27-S16
     14  314  769 (S19E50 L= 41)   732    211        5 ( 252/ 719- 365/ 810)  18.4 E56-E44 S22-S16
     15  619  743 (S21E26 L= 65)   361     75        0 ( 588/ 712- 661/ 775)  38.7 E28-E23 S23-S19
     16 1295  657 (S26W19 L=110)  1167    441        3 (1169/ 612-1496/ 743)   9.9 W10-W33 S28-S20
     17  256  610 (S29E60 L= 32)  1487    387        3 ( 226/ 511- 286/ 685)  39.8 E66-E54 S35-S24
     18  557  539 (S35E34 L= 57)   322    118        2 ( 542/ 488- 576/ 583)  60.3 E36-E31 S38-S31
     19 1046  526 (S35W03 L= 94)   334     90        0 (1001/ 492-1084/ 569)  45.0 W00-W05 S37-S32
  

    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