Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1168 1789 (N48E02 L= 11)   426     84        0 (1149/1758-1188/1820)  74.1 E03-W00 N45-N52
      2 1271 1715 (N44W08 L= 21)   290     74        1 (1244/1692-1290/1739)  45.0 W05-W09 N41-N46
      3  878 1703 (N36E21 L=351)   289     62        0 ( 862/1676- 900/1734) -45.0 E23-E19 N33-N38
      4 1460 1596 (N37W23 L= 36)   490    115        0 (1429/1559-1501/1619)  40.6 W20-W27 N33-N39
      5  344 1406 (N10E52 L=321)   472    116        2 ( 306/1387- 391/1426)   7.1 E55-E47 N09-N11
      6  400 1313 (N05E45 L=328)   958    212        4 ( 351/1232- 443/1396)  68.2 E47-E42 N00-N10
      7  732 1313 (N08E21 L=351)   276     55        0 ( 708/1288- 751/1349)  90.0 E22-E20 N06-N10
      8 1862 1181 (N17W58 L= 70)  2360    337        7 (1808/1101-1910/1256) -35.0 W51-W64 N13-N21
      9  692 1112 (S05E21 L=352)   288    104        1 ( 645/1074- 732/1150) -33.7 E24-E18 S06-S02
     10  319 1124 (S07E48 L=324)   277     46        1 ( 296/1106- 339/1143) -14.0 E50-E46 S07-S05
     11 1301  925 (S08W18 L= 31)   627    159        2 (1254/ 841-1336/1011)  79.7 W16-W19 S13-S02
     12  662  948 (S15E21 L=351)   530    157        1 ( 615/ 897- 745/ 991)  41.2 E24-E16 S18-S11
     13  128  897 (S21E67 L=305)  1316    174        2 ( 102/ 853- 166/ 944)  32.0 E70-E62 S23-S18
     14 1820  879 (S02W58 L= 70)   579    158        1 (1777/ 820-1864/ 926)  40.6 W54-W62 S05-N01
     15 1620  816 (S09W41 L= 54)   254     27        0 (1594/ 801-1640/ 833)  33.7 W39-W43 S10-S07
     16  456  705 (S32E36 L=337)   556    120        2 ( 429/ 668- 505/ 746) -21.8 E38-E31 S34-S29
     17 1067  641 (S28W08 L= 21)   332     73        1 (1045/ 602-1085/ 682) -78.7 W07-W09 S30-S25
     18 1751  576 (S20W61 L= 74)  2533    410        5 (1644/ 491-1826/ 649)  38.2 W54-W68 S27-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