Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  929 1830 (N40E36 L=103)  1660    305        3 ( 818/1789-1100/1866)  45.0 E43-E26 N34-N49
      2  942 1636 (N26E23 L=116)   284     94        0 ( 923/1593- 962/1682) -63.4 E24-E21 N23-N29
      3  507 1633 (N16E52 L= 87)   434     51        2 ( 484/1605- 528/1663)  90.0 E54-E50 N13-N18
      4 1307 1576 (N33W01 L=141)   656     76        1 (1269/1538-1343/1611) -21.8 E01-W04 N31-N35
      5 1228 1480 (N24E01 L=138)   426    178        1 (1175/1420-1290/1536)  78.7 E02-W01 N19-N29
      6 1496 1417 (N29W18 L=157)   711    169        1 (1464/1342-1534/1498) -56.3 W15-W22 N25-N34
      7  749 1425 (N08E26 L=113)   348     93        0 ( 704/1406- 804/1439)   8.1 E29-E22 N06-N10
      8  220 1178 (S16E59 L= 80)  1101    182        3 ( 183/1094- 253/1253)  73.3 E63-E56 S21-S11
      9 1843 1210 (N29W52 L=191)   458     80        0 (1818/1168-1865/1246) -26.6 W48-W54 N26-N30
     10  675 1190 (S06E24 L=115)   367     96        0 ( 642/1146- 708/1227) -26.6 E26-E20 S08-S04
     11  465 1189 (S11E38 L=101)   359     87        0 ( 441/1148- 492/1225)  78.7 E39-E36 S13-S08
     12 1772 1123 (N20W45 L=184)   510    114        0 (1740/1072-1801/1179) -29.7 W41-W48 N17-N22
     13  646 1055 (S15E22 L=117)   401     98        2 ( 601/1024- 679/1077)   0.0 E25-E19 S16-S12
     14 1628 1043 (N11W35 L=174)   359     75        0 (1592/1021-1658/1069)  45.0 W32-W36 N09-N13
     15  729 1010 (S15E16 L=123)   528    152        1 ( 680/ 970- 806/1044)  39.8 E19-E12 S18-S11
     16  896  803 (S23E01 L=138)   197     81        0 ( 872/ 786- 936/ 825)  18.4 E03-W00 S23-S21
     17 1579  661 (S11W44 L=183)   248     30        1 (1560/ 639-1596/ 685)  63.4 W43-W45 S11-S09
     18 1575  587 (S14W47 L=187)  1184    292        2 (1448/ 556-1699/ 654)  38.7 W39-W54 S19-S07
     19 1283  486 (S29W34 L=173)   255     51        1 (1256/ 466-1302/ 516) -33.7 W32-W35 S30-S27
  

    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