Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  624 1743 (N27E51 L=284)   915    113        1 ( 595/1702- 662/1789) -22.6 E56-E44 N23-N30
      2  433 1708 (N21E67 L=268)  1656    133        3 ( 375/1671- 481/1753)  54.5 E73-E62 N17-N24
      3  996 1452 (N16E13 L=322)   441     83        0 ( 966/1411-1020/1492) -78.7 E14-E10 N13-N18
      4 1697 1223 (N23W37 L= 12)  1323    324        1 (1641/1074-1754/1368) -47.1 W31-W44 N15-N30
      5  886  975 (S13E07 L=329)   256     46        1 ( 862/ 951- 918/ 992)  14.0 E08-E04 S14-S11
      6  824  923 (S18E09 L=326)   199     78        0 ( 796/ 893- 866/ 947)  63.4 E10-E06 S19-S15
      7  184  726 (S43E69 L=266)   876    102        0 ( 157/ 678- 219/ 765)  15.9 E81-E60 S46-S40
      8 1704  645 (S08W55 L= 31)   367     39        1 (1682/ 621-1725/ 671) -26.6 W53-W58 S09-S06
      9  948  615 (S32W07 L=342)   349    107        0 ( 920/ 593- 999/ 653)  -9.5 W04-W10 S33-S30
     10 1640  528 (S16W57 L= 32)   357     25        0 (1621/ 506-1664/ 543)  56.3 W54-W58 S17-S14
     11  536  397 (S61E24 L=311)   479     54        1 ( 506/ 380- 565/ 430)   5.7 E29-E19 S63-S57
     12  694  332 (S61E02 L=333)   514     88        0 ( 661/ 302- 734/ 359) -12.5 E07-W03 S62-S58
  

    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