Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  823 1812 (N38E45 L=154)   987    179        3 ( 747/1778- 913/1842)   3.2 E54-E36 N35-N39
      2  944 1726 (N35E30 L=170)   310     49        2 ( 918/1703- 968/1750) -31.0 E31-E26 N33-N36
      3  670 1722 (N27E48 L=151)   595     48        1 ( 643/1695- 710/1745)   0.0 E52-E45 N25-N28
      4  512 1667 (N19E57 L=143)   580    133        0 ( 468/1634- 571/1706)  66.8 E59-E55 N16-N23
      5  666 1520 (N14E36 L=163)   623    126        2 ( 615/1467- 711/1570)  66.8 E38-E33 N10-N17
      6 1631 1367 (N33W30 L=230)   693    195        1 (1581/1289-1658/1454) -52.1 W24-W33 N28-N37
      7  547 1287 (S03E35 L=164)  1213    255        2 ( 472/1190- 600/1403)  85.9 E38-E33 S09-N05
      8 1486 1308 (N24W20 L=219)   624    185        3 (1437/1218-1512/1381) -45.0 W15-W23 N19-N29
      9  386 1293 (S06E48 L=152)   932    182        3 ( 367/1190- 412/1367) -48.0 E51-E42 S11-S01
     10  233 1280 (S10E63 L=137)   853    148        0 ( 197/1216- 272/1343)  69.4 E65-E59 S13-S05
     11 1208 1228 (N11W05 L=204)   335     89        1 (1180/1197-1245/1264) -11.3 W02-W07 N09-N12
     12 1450 1173 (N15W20 L=220)   237     93        0 (1416/1142-1490/1202)  -9.5 W17-W23 N14-N16
     13  601  969 (S19E24 L=176)   873    302        1 ( 505/ 875- 703/1078)  62.1 E28-E19 S27-S10
     14  956  937 (S12E01 L=198)   193     82        0 ( 932/ 909- 981/ 962)  63.4 E02-W00 S13-S09
     15 1302  793 (S10W22 L=222)   660    208        2 (1238/ 757-1375/ 831)  12.5 W18-W27 S12-S07
     16 1168  776 (S15W15 L=215)   355    113        2 (1135/ 741-1213/ 815)  68.2 W13-W17 S17-S11
     17  935  646 (S29W06 L=206)  1155    348        0 ( 805/ 610-1103/ 711)  15.5 E03-W15 S31-S23
     18  782  375 (S52W07 L=207)   835    201        0 ( 693/ 351- 878/ 416)  10.0 E02-W15 S54-S48
  

    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