Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1815 (N52E27 L=199)   366     40        1 ( 943/1795- 986/1834) -21.8 E29-E24 N49-N53
      2 1031 1724 (N45E15 L=211)  1097    156        3 ( 975/1688-1095/1761)   5.2 E19-E08 N42-N48
      3 1033 1603 (N36E11 L=216)   522     92        1 ( 990/1580-1083/1638)  23.2 E13-E06 N33-N38
      4  371 1527 (N19E59 L=168)  1190    227        0 ( 302/1446- 443/1615) -13.2 E81-E47 N15-N24
      5  281 1391 (N10E61 L=165)   796     95        3 ( 247/1356- 307/1430)  59.0 E64-E57 N07-N12
      6 1323 1180 (N13W16 L=242)   837    181        3 (1253/1133-1382/1223)  45.0 W11-W19 N09-N16
      7 1251 1112 (N08W12 L=239)   712    213        0 (1195/1047-1311/1215) -41.6 W07-W16 N05-N13
      8  753 1176 (N04E19 L=208)   280     61        0 ( 723/1156- 786/1198) -14.0 E20-E16 N03-N05
      9 1041 1124 (N06E01 L=226)   461    130        0 (1024/1060-1061/1193) -74.1 E02-W01 N02-N09
     10  385 1012 (S11E42 L=185)   862    266        1 ( 369/ 883- 415/1127) -86.2 E43-E39 S18-S03
     11 1544 1016 (N07W33 L=259)   295    128        0 (1494/ 979-1579/1037)  -9.5 W29-W35 N05-N08
     12  297  859 (S22E49 L=177)  1071    167        1 ( 214/ 819- 380/ 881)  14.7 E59-E40 S25-S19
     13 1259  722 (S15W20 L=246)   481    134        0 (1223/ 661-1294/ 773)  81.9 W17-W20 S18-S11
     14 1655  613 (S15W52 L=278)   350     75        0 (1612/ 590-1683/ 634)  -7.1 W47-W55 S15-S13
     15 1484  416 (S30W47 L=274)   446     76        0 (1445/ 395-1525/ 441)  45.0 W44-W50 S32-S27
     16  583  397 (S50E26 L=200)   397     55        1 ( 560/ 374- 605/ 426) -76.0 E28-E23 S52-S47
     17  677  313 (S56E17 L=210)   852     92        1 ( 647/ 277- 718/ 356) -50.2 E20-E11 S58-S51
     18  959  236 (S56W16 L=242)  1107    140        4 ( 904/ 187-1015/ 292) -25.2 W08-W25 S60-S51
  

    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