Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1117 1888 (N57E34 L=  1)   877     56        1 (1087/1861-1147/1914) -13.0 E41-E28 N53-N60
      2 1483 1712 (N56W13 L= 49)   669     81        1 (1456/1665-1515/1738)  68.2 W10-W17 N51-N61
      3 1692 1574 (N53W41 L= 77)   713     80        1 (1667/1535-1714/1611)  81.9 W36-W45 N49-N56
      4 1053 1427 (N19E09 L= 26)   287    113        1 (1002/1409-1103/1450)   9.5 E12-E06 N17-N20
      5 1628 1252 (N26W32 L= 67)   425    150        0 (1567/1216-1656/1306)  58.0 W27-W34 N21-N29
      6 1640 1115 (N18W35 L= 71)   454    131        0 (1616/1049-1674/1164)  76.0 W34-W37 N13-N21
      7 1769 1043 (N18W48 L= 83)   466     69        3 (1739/1018-1800/1063)   0.0 W45-W51 N16-N19
      8 1890  972 (N18W64 L=100)  1470    121        3 (1862/ 915-1926/1028)  58.0 W60-W71 N14-N22
      9 1729  966 (N12W46 L= 82)   719    139        1 (1670/ 932-1777/1005)  32.0 W41-W50 N09-N15
     10  243  937 (S29E54 L=341)  1058    201        2 ( 205/ 819- 275/1002) -74.7 E59-E50 S35-S24
     11 1830  882 (N11W59 L= 94)   933     83        2 (1792/ 853-1868/ 914)  18.4 W55-W64 N08-N13
     12  367  860 (S31E40 L=355)   464    143        1 ( 307/ 827- 436/ 887)  22.6 E46-E34 S33-S27
     13  192  831 (S37E63 L=333)   664     77        1 ( 168/ 794- 216/ 859) -63.4 E66-E58 S39-S35
     14 1200  408 (S35W33 L= 69)   630     71        1 (1160/ 389-1247/ 438)  31.0 W29-W36 S36-S32
     15  984  415 (S41W19 L= 54)   311     91        0 ( 934/ 402-1035/ 426)  20.6 W14-W22 S43-S39
     16  503  394 (S62E27 L=  9)   467     68        0 ( 471/ 367- 532/ 414) -12.5 E31-E22 S63-S59
  

    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