Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1199 1786 (N48W01 L= 11)   439     81        1 (1178/1755-1227/1812) -53.1 E01-W03 N45-N50
      2 1495 1593 (N38W26 L= 36)   491     85        0 (1468/1560-1528/1617)  39.8 W23-W29 N34-N39
      3  571 1535 (N20E37 L=333)   262     63        2 ( 547/1514- 601/1563)  18.4 E38-E35 N18-N22
      4  425 1422 (N12E45 L=325)   659    222        2 ( 333/1390- 515/1459)  12.1 E52-E38 N09-N14
      5  442 1341 (N07E42 L=328)   408    114        2 ( 407/1305- 475/1395)  90.0 E44-E39 N04-N10
      6  769 1301 (N07E19 L=351)   212     59        0 ( 738/1281- 791/1329)  63.4 E20-E17 N05-N09
      7 1876 1202 (N19W59 L= 69)  1133    211        3 (1836/1166-1917/1257) -26.6 W54-W64 N16-N21
      8  345 1141 (S06E46 L=324)   398     51        1 ( 323/1110- 370/1169) -45.0 E47-E44 S07-S04
      9 1890 1127 (N14W61 L= 71)   705    103        2 (1870/1088-1911/1158) -36.9 W58-W64 N12-N16
     10  727 1108 (S05E19 L=351)   256     76        0 ( 685/1079- 763/1143) -31.0 E21-E16 S06-S02
     11  697  938 (S15E19 L=351)   593    159        1 ( 647/ 882- 774/ 985)  45.0 E21-E14 S18-S11
     12  142  895 (S21E65 L=305)  1313    192        1 ( 108/ 852- 191/ 950)  36.9 E69-E59 S23-S17
     13 1325  879 (S10W20 L= 30)   407    107        1 (1280/ 830-1361/ 932)  74.1 W17-W21 S13-S06
     14  311  879 (S22E47 L=323)   270     74        1 ( 293/ 857- 343/ 907) -26.6 E48-E44 S23-S20
     15 1652  809 (S09W44 L= 54)   242     36        0 (1627/ 795-1674/ 824)  14.0 W41-W45 S09-S07
     16  485  702 (S32E33 L=337)   498    124        3 ( 449/ 668- 532/ 742) -26.6 E36-E29 S33-S29
  

    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