Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1077 1847 (N47E28 L=111)   346     68        0 (1054/1822-1111/1859)  26.6 E29-E25 N44-N49
      2  887 1823 (N38E38 L=100)  1168    168        2 ( 823/1789- 979/1845)  52.1 E42-E33 N34-N43
      3  951 1638 (N27E22 L=116)   306     90        0 ( 931/1595- 972/1688) -74.1 E23-E21 N23-N30
      4  513 1633 (N16E52 L= 87)   415     55        0 ( 491/1603- 534/1665) -68.2 E53-E49 N13-N18
      5 1317 1575 (N33W02 L=141)   632     95        1 (1281/1527-1350/1614) -45.0 W00-W05 N30-N35
      6 1234 1472 (N24E00 L=138)   551    206        0 (1184/1390-1299/1532)  84.8 E02-W02 N18-N29
      7 1503 1415 (N29W18 L=157)   609    149        1 (1474/1342-1540/1488) -48.8 W15-W22 N25-N33
      8  755 1422 (N08E26 L=113)   381     86        0 ( 705/1405- 813/1438)  18.4 E28-E22 N06-N10
      9  235 1214 (S14E59 L= 80)   614     83        1 ( 204/1177- 256/1248)  90.0 E62-E56 S15-S11
     10  682 1188 (S06E23 L=115)   403    113        0 ( 649/1140- 714/1237) -33.7 E26-E20 S08-S04
     11  471 1185 (S11E37 L=101)   368     85        0 ( 445/1147- 501/1222)  68.2 E38-E35 S13-S08
     12 1779 1118 (N20W46 L=184)   422    123        0 (1747/1068-1804/1177) -29.7 W41-W48 N17-N22
     13  649 1051 (S15E22 L=116)   381    105        1 ( 600/1020- 684/1072)   0.0 E25-E19 S16-S12
     14 1636 1043 (N11W35 L=174)   334     60        0 (1605/1019-1664/1067)  36.9 W33-W37 N09-N13
     15  731 1004 (S15E16 L=122)   493    136        1 ( 687/ 965- 798/1040)  45.0 E18-E12 S18-S11
     16 1579  584 (S15W48 L=186)  1163    294        1 (1453/ 550-1704/ 649)  40.6 W40-W54 S19-S07
     17 1290  482 (S29W34 L=173)   246     55        1 (1265/ 461-1306/ 512) -26.6 W32-W36 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