Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  813 1749 (N32E35 L=314)   605    101        0 ( 778/1701- 842/1785) -24.0 E39-E30 N28-N34
      2  966 1593 (N25E18 L=331)   323    126        0 ( 947/1535- 982/1651) -54.5 E20-E15 N21-N28
      3  834 1546 (N19E24 L=325)   424    140        1 ( 762/1529- 908/1565)  12.5 E28-E19 N17-N19
      4 1561 1361 (N26W25 L= 14)   550    172        1 (1525/1292-1612/1406)  77.5 W23-W28 N21-N30
      5 1655 1316 (N26W33 L= 22)   382    116        0 (1612/1263-1683/1354) -26.6 W29-W35 N23-N27
      6 1283 1305 (N15W08 L=357)   292    124        2 (1257/1266-1331/1339) -11.3 W06-W11 N11-N16
      7  954 1252 (N04E09 L=340)   775    339        0 ( 803/1212-1101/1310)  29.4 E17-E01 N00-N10
      8 1591 1212 (N18W29 L= 19)   601    159        1 (1543/1175-1666/1248)  11.3 W25-W35 N15-N19
      9 1830 1012 (N13W53 L= 42)   818    153        0 (1816/ 941-1851/1090) -83.7 W51-W55 N09-N18
     10  194  982 (S25E60 L=289)  1720    187        6 ( 122/ 929- 261/1016)  11.9 E71-E51 S27-S23
     11 1819  805 (N02W58 L= 47)   425     88        1 (1780/ 771-1851/ 837)  53.1 W55-W62 N00-N04
     12  196  671 (S45E68 L=281)  1369    114        2 ( 135/ 640- 253/ 708)  -6.6 E87-E58 S46-S42
     13 1140  657 (S25W16 L=  6)   327     44        0 (1112/ 633-1172/ 683) -14.0 W14-W18 S25-S23
     14 1448  524 (S24W41 L= 31)  2004    482        4 (1250/ 470-1678/ 589)  35.8 W30-W56 S31-S13
     15 1761  490 (S15W72 L= 62)   899     74        1 (1722/ 472-1798/ 514)  21.8 W66-W79 S17-S12
     16 1365  222 (S43W62 L= 51)   506     58        2 (1341/ 189-1382/ 252)  -9.5 W57-W69 S44-S41
  

    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